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ABSTRACT 


This report summarizes the econometric theory surrounding estimatijon of 
residential demand functions and describes the methodology used to develop 
residential electricity demand functions for 136 utility service areas. A 
Single equation function is described; its functional form is log-linear. 
Special treatment is given cost-of-living differences and climatic variables. 
The demand equations were estimated in two stages. The selected model 
contained humidity as a second weather variable, is adjusted for 
heteroskedasticity, and has all variables correct in Sign and approximate 
Magnitude. Elasticities are reported for income and prices and compared to 
previous studies. Price and income projections are then used to predict 1985 
electricity sales. 


Keywords: electricity, demand econometrics, coal, households, model 


This paper was prepared for limited distribution to the research community 
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Authors 


Dr. John R. McKean is a Professor in the Economics Department at Colorado 
State University, Fort Collins and was Principal Investigator on the contract 
which developed this work. 

Dr. John W. Green is a Regional Economist and Project Leader, Natural 
Resource Economics Division, Economic Research service, U.S. Department of 
Agriculture located in the Economics Department, Colorado State University, 
Fort Collins, Colorado and contract monitor for this work. 


11 


FOREWORD 


The research activities and results described in this report are part of a 
project entitled "Integrated Assessment: Economic and Social Consequences of 
Coal and Oil Shale Development." This project is supported jointly by the 
U.S. Department of Agriculture (USDA) and the Environmental Protection Agency 
(EPA) as a component of the Federal Interagency Energy/Environment Research 
and Development Program. The overall project focuses on identifying and 
analyzing certain interrelationships between society's needs for energy, 
environmental quality, and agricultural production and the impacts of 
alternative public policy strategies for dealing with these needs. 


The econometric model developed in this contractor's report is both 
utility service area and national in scope. It illustrates the logic and 
structure of an econometric household electricity demand model useful in 
deriving the demand for coal to produce electricity. The model development 
work described in this report has been halted. Our purpose of this report is 
to describe completed work, to simulate discussion of the basic concepts, and 
to solicit comments and suggestions. Colleagues in the energy demand modeling 
field, as well as individuals actively involved in related policy analysis, 
are encouraged to send their comments to the authors or to us. 


Paul Schwengels, Project Officer Joseph R. Barse 
Office of Environmental Engineering Natural Resource Economics 
and Technology Division 
Office of Research and Development U.S. Department of Agriculture 


U.S. Environmental Protection Agency 
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SUMMARY 


The average residence in the U.S. consumes about 660 kWh of electricity 
per month, more than two-thirds of which is produced using coal as a fuel. 
Power plants require about 65 gallons of oil, 660 pounds of coal, or 9,000 
cubic feet of natural gas to generate that 660 kWh of electricity each month. 
Income and other factors play a role in residential electricity consumption 
but the relationship is not as direct as the relationship between production 
and electricity use in the industrial sector. 


The goal of the project reported in this study was to provide a model that 
could be used to project electricity consumption by households in specific 
service areas. The projections can be used in the Interregional Coal Analysis 
Model to evaluate impacts of coal development on rural resources and 
communities. 


A single equation least-squares demand function was selected to project 
household electricity consumption in 1985. The functional form of the demand 
equation is log-linear and special treatment is given variables describing 
service areas and climate. The primary data sources include Federal Energy 
Regulatory Commission, Bureau of Census, County and City Data Book, American 
Gas Association, and Department of Commerce-Department of Agriculture OBERS 
projections. 


Existing studies of residential electricity demand are extremely 
heterogeneous in the structuring of their models and in their research 
findings. All studies agree that prices are important but different 
methodologies are employed to describe buyer and seller behavior. Several 
controversies concerning model specification exist among researchers. These 
include marginal versus average price, simultaneous supply and demand 
modeling, and impacts of historical fuel prices. 


Over three-fifths of previous studies use either time series or pooled 
cross section-time series methods. The pooled and pure time series models 
invariably use the distributed lag technique to model the demand function. 
Adoption of this technique allows the use of short run time series data to 
estimate both short run and longer run price elasticities. Both time series 
and pooled time series-cross section analyses also incorporate distributed lag 
relationships in order to remove serial correlation. Unfortunately, the use 
of time series data presents several problems, including trends in the 
variables. It is impossible to say to what degree the relationship among 
these variables is spurious rather than reflecting a true interdependence. 
Thus, time series data present some insoluble problems. Improper 
specification of the model can result in bias in the elasticity estimates. 


Time series data also suffer from collinearity since nominal "independent" 
variables contain trends which make them correlate highly to each other. 
Cross section analyses avoids the problem of contemporaneous trends although 
some collinearity may still exist. 


Having rejected the time series and dynamic modeling approaches implied by 
time series data, only the pure cross section alternative remains. The 
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approach chosen utilizes cities making up utility service areas as sample 
observations. Thus, a cross section regression equation is estimated across 
selected cities in the U.S. The statistical model is consistent with the 
comparative static framework of analysis whose properties are more completely 
understood than dynamic economic models. 


The model incorporates both monetary and physical variables, i.e., 
electricity consumption is affected by the prices of electricity and closely 
related goods, consumer income, and by nonmonetary variables such as climate 
and weather. All prices and income are expressed in real terms to remove 
unwanted variation caused by changes in the cost of living. 


The estimating procedure was undertaken in two stages. The first stage 
involved the examination of several alternative measures of the hypothesized 
causal variables for 102 cities. Demand equations were specified and a 
preliminary recursive model approximated for further testing uSing better 
American Gas Association consumer price data. The second stage consisted of 
reestimating the demand equations specified in the first stage using a 
surrogate price index to adjust the monetary variables in an expanded 138 city 
sample necessitated by the need for forecasts for utility service areas. The 
use of the surrogate price deflator index allowed the expansion of the 
deflated demand estimates from a sample of 57 to the total of 138 cities. The 
use of deflated cross section data greatly improved the accuracy of the 
estimated equations. Average real spending by city per family for household 
electricity in 1975 varied from $151 to $545. The mean deflated electricity 
price of the sample data for 1975 was $38.64 per thousand kWh, the mean 
deflated household income was $11,854. The estimated model was very robust. 


Six equations are presented in the second stage of the empirical 
investigation. All the final equations may be described as log-linear. 
Separate and pooled regressions are presented to allow adjustment for bias 
resulting from inclusion of surrogate deflators. Elasticities are shown for 
electricity price, income, and natural gas price by city. The demand 
equations are shown in both log-linear and product form on page 89 and will 
not be repeated here. Tables 10 through 15 present 6 demand equations. 
Coefficients of determination range from 0.842 to 0.978. 


Model I is the demand equation for the subsample of 57 cities for which 
actual price deflators were available. Model II is similar except it was 
fitted to the sample of 81 cities for which the surrogate inflator was 
required. Models I and II incorporate dummy variables which shift the demand 
equation up or down to adjust for variation among states in the rate of 
household electricity consumption. 


Model III combines the sample data included in Models I and II fora 
sample size of 138 cities. The sign on the natural gas variables is “wrong” 
if natural gas is to be considered a substitute for electricity; a result 
which has been reported in other studies. Model VI does not differentiate 
among real or surrogate price deflators and contains humidity as a second 
weather variable. Model VI is also adjusted for heteroskedasticity. 


Models IV and V are similar to Models I and II except that the stepwise 
regression program was calibrated to allow state price elasticity dummy 


variables to enter. This allows both the slope and intercept coefficients to 
vary by state. Inclusion of both effects leads to high correlation among 
independent variables and possible estimation errors. The net addition to 
explanation resulting from addition of independent variables in models IV and 
V is quite small. 


Model VI, which pools the two data sets contained in Models I and II, has 
been selected as the prediction equation best suited to provide consumption 
forecasts for 136 utility service areas. Model VI is superior to Model III 
which estimated the wrong sign for natural gas prices deflated by the "real" 
price index. The use of weighted least squares to adjust for 
heteroskedasticity increased the efficiency of Model IV and a second weather 
variable was also significant. All signs were theoretically correct and 
reasonable in magnitude. The historical data, projected inputs, econometric 
model, and consumption forecasts are summarized in Figure 5. 


Our price elasticity estimate of -0.72 is slightly lower than those in 
previously reported studies (Tables 2 and 3). Our cross elasticity to natural 
gas price (0.10) is also lower than previous estimates. Our income elasticity 
of 0.35 is also below the average of other pure cross section models. 

Although our elasticity estimates tend to be smaller in absolute value than 
previously reported estimates, few of the estimates differ markedly. Several 
factors could account for our smaller elasticity estimates; they are discussed 
beginning on page 101. 


Department of Energy electricity price change projections were applied to 
our electricity price data in order to utilize our estimated demand equation 
to predict future electricity sales. DOE historical and projected natural gas 
prices were used in the same manner. Family income projections were derived 
using Bureau of Economic Analysis data. Long run average heating degree days 
was used as the predictor of heating degree days in 1985. Thus, the separate 
influences of changes in a particular city's energy price structure and income 
level on electricity consumption are indicated. The effects of changes in 
prices and incomes depend both on the elasticity estimated by the regression 
model and on the amount of price and income changes projected by DOE and BEA 
for 1985. 7 


The primary goal of the research presented here was to project household 
electricity demand in 1985 for specific utility service areas. For projection 
purposes, we have utilized a simplified model which combines changes in 
appliance stocks and changes in utilization rates into the price effect. But 
we have also determined that long run electricity demand projections should 
consider changes in shares of appliances powered by electricity. Further 
research support would enable us to (1) continue the development of our data 
set, (2) more fully utilize the superior AGA substitute fuel price data, (3) 
continue the development and testing of several alternative models of the 
market for residential electricity, and (4) estimate the coal-fired 
electricity demand of those utilities not represented by SMA data. A two 
equation model should be developed in which the first equation relates 
quantity of electricity demanded per customer to the price of electricity, 
weather conditions, and the level of saturation of ownership of residential 
electrical appliance stocks. The second equation would relate the saturation 
levels of ownership to the price of electricity, the prices of competing 
fuels, income, climate, and other regional socioeconomic characteristics. 

Some exploratory modeling has been performed along these lines and is reported 
in Appendix II. 
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COAL FROM THE NATION: HOUSEHOLD DEMAND FOR ELECTRICITY 


John R. McKean* and John W. Green** 


CHAPTER I. INTRODUCTION 


Household Electricity Demand and Coal 


Residential electricity consumption is already the largest single source 
of demand for coal in the United States. More than two-thirds of U.S. coal 
production is used to generate electricity. Potentially, much larger demands 
for coal may occur as alternative fuels grow more expensive to extract or to 
import and if nuclear generation continues to lose appeal. Household use of 
electricity provides a "derived" demandl/ for coal. The average residence 
in our nation consumes about 660 kWh per month of electricity. In order for 
an electric utility to supply a household's monthly demand, approximately 65 
gallons of 011, 660 pounds of coal, or 9,000 cubic feet of natural gas must be 
consumed at the generating plant. Additional amounts of energy are used 
indirectly to mine and transport the fuels used to power the generators and to 
provide the human and capital inputs necessary to the electric power industry 
infrastructure. Electricity provided about 14 percent of the residential 
energy requirements in 1978. Electric utilities accounted for 30 percent of 
total energy use in the U.S. in 1978. 


The nature of consumer spending is such that forecasts are necessarily 
more complex than for industrial spending. Coal or electricity is a factor of 
production purchased by manufacturers in proportion to the value of their 
output. Industrial energy use is closely related to production. Although 
consumer spending is ultimately constrained by the value of production, i.e., 
by their income, consumers in the U.S. normally have a large amount of 
discretionary income to spend as they desire. Income plays a role in 
residential electricity consumption but there is no direct tie between 
production and electricity use as exists in the industrial sector. 

Electricity as a final demand sale to households is subject to a large variety 
of influences. 


Electricity Demand Growth 


An energy study (Sponsored in part by the Ford Foundation Energy Policy 
Project) by E. A. Hudson and D. W. Jorgenson (67) utilized a pioneering 
approach combining econometric models with input-output (Leontief) models to 


T/ Coal is consumed indirectly by households who purchase electricity which 
is produced from coal. 

*Dr. John R. McKean is a Professor in the Economics Department at Colorado 
State University, Fort Collins and was Principal Investigator on the contract 
which developed this work. 

**Dr. John W. Green is a Regional Economist and Project Leader, Natural 
Resource Economics Division, Economic Research Service, U.S. Department of 
Agriculture located in the Economics Department, Colorado State University, 
Fort Collins, Colorado and contract monitor for this work. 
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relate patterns of economic growth to the demand and supply of energy .2/ 

The study simulates the workings of the U.S. economy with exogenous impacts by 
foreign sectors. These two authors have used the model to study the effects 
of tax policies on future patterns of energy consumption. Their comment on 
their electricity consumption "base case" under existing conditions is 
revealing. 


"Electricity consumption increases at a rapid rate--at 5.2 percent a 
year between 1975 and 2000. ... Household electricity consumption 
increases at an average annual rate of 6.9 percent between 1975 and 
2000, with household and industrial demand increasing at an average of 
5.9 percent and industrial demand at 3.8 percent. The increase in 
household use of electricity is the most rapid of any use of any fuel 
and, in this, it continues the trend that has been observed in the 
past. However, the projected rate of increase is slower than those 
observed in the past due to the more rapid increase in electricity 
prices, the approach to saturation of heavy electricity using 
appliances such as home air-conditioning, laundries, dishwashers and 
so on, and to the decline in the rate of new dwelling construction." 


Table 1 depicts the Hudson-Jorgenson "base case" U.S. energy use 
projections. 


Purpose of This Study 


As opposed to many other empirical studies of residential electricity 
demand in the United States, this study has a very limited goal. Our task is 
to provide a model that can be used to project electricity consumption by 
households in certain specified service areas (e.g., cities) from 1975 to 1985 
(see Table 2). These projections can be used in an extensive linear 
programming model concerned with coal utilization in the u.S.3/ The nature 
of the econometric model which has been selected to perform this projection 
task may differ considerably from that which would be desired if the task were 
to test a hypothesis or to examine the general characteristics of markets. 

Two conflicting qualities are desired in the projections model. First, it 
should provide a good fit or explanation of historical data which are used to 
estimate the parameters of the model. Second, it should be simple and contain 
a minimum of independent variables. This second requirement is imposed 


2/ The nine-sector I-O model has technical coefficients which are 
endogenously determined as a function of prices of outputs and inputs. The 
second component of the model incorporates a dynamic macroeconometric system 
to project components of national income. 

3/ The national Interregional Coal Analysis Model is a large-scale linear 
programming model intended to assist in coal development planning and 
analysis. Coal development has expanded considerably in the mountain states 
and has caused severe impacts on agriculture and rural communities. The 
interregional competition model can be used to systematically assess the 
impacts of alternative policy decisions on coal development in the United 
States. 

The interregional competition model is described in a report by John W. 
Green, Western Energy: The Interregional Coal Analysis Model August, 1980, 
Natural Resource Economics Division, Economic Research Service, U.S. 
Department of Agriculture, Technical Bulletin No. 1627. 
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Table 1--Composition of U.S. energy use; 1970-2000 (energy flows in trillion Btu) 


Consuming Sector 
Household/ Indus- Trans- Electricity Synthetic Total 


Fuel Commercial trial portation generation gas Input 
1970 Coal 427 5,004 8 7,483 0 E2922 
Petroleum 6,453 5. 207) 15. 592 2,087 0 29,614 
Natural Gas 7,108 10,162 745 4,015 0 22,029 
Hydro, Nuclear 0 0 0 2,879 0 2,879 
Total Input 13,988 20,433 16,345 16,464 0 67,444 
Electricity 3,000 2,210 16 0 0 5,206 
Synthetic Gas 0 0 0 0 0 0 
Total 16, 988 22.045 2) 16,3561 16, 464 0 67,444 
1975 Coal 200 4,393 0 8,625 0 13, 219 
Petroleum 6,819 62627" 2 a17..9357 2,400 0 355/65 
Natural Gas 8,415 1a NGZ 945 4,290 0 ZATSS 
Hydro, Nuclear 0 0 0 5.475 0 5475 
Total Input 15,4355 22,122, PabS, B82 20,790 0 T1520 
Electricity 3,982 2,709 22 0 0 6,713 
Synthetic Gas 0 ‘V0 0 0 0 0 
Total 19, 418 24,831 18,904 20, 790 0 77,230 
1980 Coal 180 3,897 0 LOK227 0 14,304 
Petroleum 7,996 7,665 20,967 2,567 0 39,195 
Natural Gas 8,515 12,778 1,003 4,505 0 26, 800 
Hydro, Nuclear 0 0 0 TOSEES 0 LOSLES 
Total Input 16,691 24,340 21,970 27,413 0 90,414 
Electricity 5,427 3,396 30 0 0 8,853 
Synthetic Gas 0 0 0 0 0 0 
Total 22,118 21,7350 | 22,000 2 ALS 0 90,414 
1985 Coal 175 S607 0 11,779 Dad 16,741 
Petroleum 9,149 8,938 23,948 2,907 89 45,031 
Natural Gas 7,499 14,690 993 4,810 0 28,022 
Hydro, Nuclear 0 0 0 16,010 0 16,010 
Total Input 16,822 27,304 24,942 35,536 1,199 105,804 
Electricity 7,490 3,951 37 0 0 11,478 
Synthetic Gas 310 379 0 0 0 690 
Total 24,623 31,635 24,978 35,536 1,199 105, 804 
1990 Coal SS S002 0 13,545 2,594 19, 887 
Petroleum 9,964 10,655 25,902 4,053 208 50, 783 
Natural Gas 6,462 15,978 1,000 SS 0 28,612 
Hydro, Nuclear Die: 0 0 23,616 0 23,616 
Total Input 16,561 30,245 26,902 46,387 2,802 122,897 
Electricity 10,064 4,876 43 0 0 14,984 
Synthetic Gas 808 988 0 0 0 1,798 
Total 27,434 36,109 26,945 46, 387 2,802 122,897 
2000 Coal 174 5,048 0 15,878 5,061 26,160 
Petroleum 2727 14,405 29,810 4,085 405 61,432 
Natural Gas 4,742 16,632 954 6,310 0 28,639 
Hydro, Nuclear 0 0 0 AT ASZ 0 A7 152 
Total Input 17,643 36,085 30,764 73,425 5,466 163,383 
Electricity 16,737 6,927 56 0 0 23,720 
Synthetic Gas 1,638 2,002 0 0 0 3,640 
Total 36,017 45,014 30,820 73,425 5, 466 163,383 


ource: Econometric Studies of U.S. Ener olicy, Dale W. Jorgenson, Editor, 
North-Holland Publishing Company, Amsterdam, 1976. 
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Table 2--Standard Metropolitan Areas appearing in the 


ALABAMA 
Birmingham 

ALASKA 

ARIZONA 


Phoenix 
Tucson 


ARKANSAS 
CALIFORNIA 


Anaheim 

Los Angeles 
San Diego 
San Francisco 
San Jose 


COLORADO 


Colorado Springs 


Denver 
CONNECTICUT 
Bridgeport 
Hartford 
Waterbury 
DELAWARE 


Wilmington 


DISTRICT OF COLUMBIA 


Washington 
FLORIDA 


Fort Lauderdale 
Orlando 


GEORGIA 


Atlanta 


residential demand analysis 


Alabama Power Company 


None 


Salt River Power District 
Tucson Gas and Electric Co. 


None 


City of Anaheim Light and Power 
L.A. Dept. of Water and Power 
San Diego Gas and Electric Co. 
Pacific Gas and Electric Co. 
Pacific Gas and Electric Co. 


Colo. Springs Dept. of Public Utilities 
Public Service of Colorado 


The United Illuminating Co. 
The Hartford Electric Light Co. 
The Conn. Light and Power Co. 


Delmarva Power and Light Co. 


Potomac Electric Power Co. 


Florida Power and Light Co. 
Orlando Utilities Commission 


Georgia Power Co. 


Continued 


Table 2--Continued 


IDAHO 
Boise City 
ILLINOIS 
Chicago 
Decatur 
Peoria 
INDIANA 
Evansville 
Fort Wayne 
Indianapolis 
Lafayette 
Muncie 
IOWA 


Des. Moines 
Sioux City 


KANSAS 
Topeka 
KENTUCKY 
Lexington 
Louisville 
Owensboro 
LOUISIANA 
Lafayette 
Lake Charles 
New Orleans 
MAINE 
Portland 
MARYLAND 


Baltimore 


Idaho Power Co. 


Commonwealth Edison Co. 
Illinois Power Co. 
Central Illinois Light Co. 


Southern Indiana Gas and Electric Co. 
Indiana and Michigan Electric Co. 
Indianapolis Power and Light Co. 
Public Service Co. of Indiana, Inc. 
Indiana and Michigan Electric Co. 


Iowa Power and Light Co. 
Iowa Public Service Co. 


The Kansas Power and Light Co. 


Kentucky Utilities Co. 
Louisville Gas and Electric Co. 
Owensboro Municipal Utilities 


Lafayette Utilities System 
Gulf States Utilities Co. 
New Orleans Public Service, Inc. 


Central Main Power Co. 


Baltimore Gas and Electric Co. 


Continued 


Table 2--—Continued 


MASSACHUSETTS 
Boston 
Fall River 
Pittsfield 
Springfield 
MICHIGAN 
Flint 
MINNESOTA 
Duluth 
Minneapolis 
Rochester 
MISSISSIPPI 
Jackson 
MISSOURI 
Columbia 
Kansas City 
St. Louis 
MONTANA 


Billings 
Great Falls 


NEBRASKA 


Lincoln 
Omaha 


NEVADA 
NEW HAMPSHIRE 


Manchester 
Nashua 


NEW JERSEY 


Atlantic City 


Newark 
NEW MEXICO 


Albuquerque 


Boston Edison Co. 

Fall River Electric Light Co. 
Western Massachusetts Electric Co. 
Western Massachusetts Electric Co. 


Consumers Power Co. 


Minnesota Power and Light Co. 
Northern States Power Co. 
Rochester Dept. of Public Utilities 


Mississippi Power and Light Co. 


Columbia Water and Light Dept. 
Kansas City Power and Light Co. 
Union Electric Co. 


The Montana Power Co. 
The Montana Power Co. 


Lincoln Electric System 
Omaha Public Power District 


None 


Public Service Co. of New Hampshire 
Public Service Co. of New Hampshire 


Atlantic City Electric Co. 
Public Service Electric and Gas Co. 


Public Service Co. of New Mexico 


Continued 


Table 2--Continued 


NEW YORK 


Binghamton 
Buffalo 
New York City 
Rochester 
NORTH CAROLINA 
Charlotte 
Fayetteville 
Raleigh 
NORTH DAKOTA 
Fargo 
OHIO 


Cincinnati 
Cleveland 
Columbus 
Dayton 
OKLAHOMA 
Oklahoma City 
OREGON 
Portland 
PENNSYLVANIA 
Erie 
Harrisburg 
Philadelphia 
Pittsburg 
Scranton 
RHODE ISLAND 
Providence 


SOUTH CAROLINA 


Columbia 
Greenville 


SOUTH DAKOTA 


Sioux Falls 


New York State Electric & Gas Corp. 
Niagara Mohawk Power Corp. 
Consolidated Edison Co. of N.Y., Inc. 
Rochester Gas and Electric Corp. 


Duke Power Co. 
Fayetteville Pub. Works Commission 
Carolina Power and Light Co. 


Northern States Power Co. 


The Cincinnati Gas and Electric Co. 
The Cleveland Electric Illuminating Co. 
Columbus and Southern Ohio Electric Co. 
The Dayton Power and Light Co. 


Oklahoma Gas and Electric Co. 


Portland General Electric Co. 


Pennsylvania Electric Co. 
Pennsylvania Power and Light Co. 
Philadelphia Electric Co. 
Duquesne Light Co. 

Pennsylvania Power and Light Co. 


The Narragansett Electric Co. 


Southern Carolina Electric and Gas Co. 
Duke Power Co. 


Northern States Power Co. 


Continued 


Table 2--Continued 


TENNESSEE 
Knoxville 
Memphis 
Nashville 


TEXAS 


Dallas 

El Paso 
Fort Worth 
Galveston 
Houston 
Lubbock 

San Antonio 


UTAH 


Salt Lake City 


VERMONT 
VIRGINIA 
Norfolk 
Richmond 
Roanoke 
WASHINGTON 
Seattle 
WEST VIRGINIA 
Charleston 
WISCONSIN 
Appleton 
Green Bay 
Milwaukee 


Racine 


WYOMING 


Knoxville Utilities Board 


Memphis Light Gas and Water Division 


Nashville Electric Service 


Dallas Power and Light Co. 

El Paso Electric Co. 

Texas Electric Service Co. 
Houston Lighting and Power Co. 
Houston Lighting and Power Co. 
Southwestern Public Service Co. 
San Antonio Public Service Board 


Utah Power and Light Co. 


None 


Virginia Electric and Power Co. 
Virginia Electric and Power Co. 
Appalachian Power Co. 


Seattle Dept. of Lighting 


Appalachian Power Co. 


Wisconsin Michigan Power Co. 
Wisconsin Public Service Corp. 
Wisconsin Electric Power Co. 
Wisconsin Electric Power Co. 


None 
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because to use a model to make projections requires substituting estimated 
1985 values for each independent variable into the projection equation. Each 
of these estimated independent variable values is subject to sampling error. 
Ideally, household electricity demand would only depend on a few rather stable 
variables such as climate or population. Unfortunately, recent history and 
empirical studies reveal that incomes and prices have considerable influence 
on electricity consumption. Forecasts cannot be expected to succeed if these 
important and volatile influences are ignored. Nevertheless, this study 
examines several types of econometric models both conceptually and empirically 
and when an equally good fit can be obtained with a less "realistic" model 
(generally a simpler model) the simpler model is selected. Reducing the 
number of independent variables which must be projected to 1985 should reduce 
the amount of sampling error in the forecasts. 


Methodology and Data 


In keeping with the requirements imposed by the projections goal, a single 
equation least-squares demand function has been selected to project household 
electricity consumption for 1985. The functional form of the demand equation 
is log-linear, that is, each of the continuous variables is subjected toa 
logarithmic transformation and a linear equation is estimated relating the 
logarithm of quantity demanded per customer to the logarithms of several 
continuous independent variables. Special treatment is given certain dummy 
variables which designate regions or states which encompass service areas. A 
different functional form is also applied to the climate variable, heating 
degree days. Details of the equation form are discussed in Chapter III. 


Two separate data inputs are necessary to construct projections of 
household electricity consumption by city in 1985. The first data requirement 
is to provide a sample to statistically estimate a forecasting equation. In 
this case we desire to estimate a demand equation relating quantity per 
household to causal variables such as prices, income, climate and the like. 
The second data requirement is to provide estimates of the values assumed by 
the causal variables in 1985. In some cases, the latter requirement is very 
Simple. For example, the appropriate measure of 1985 heating degree days in a 
given city is likely to be best predicted by the historical average heating 
degree days for that city. Projection of other independent variables, such as 
prices, is not as obvious. Further discussion of this problem occurs in 
Chapter V. 


The primary sources of data used to build the econometric model include 
Federal Energy Regulatory Commission (FERC) reports, the Census, County and 
City Data Book and uncirculated data provided by the American Gas 
Association.4/ Numerous other government publications were required and 
these are documented in the bibliography. Sources of projections for the 
independent variables include price projections by the U.S. Department of 
Energy and population and income projections by the Department of Commerce and 
U.S. Department of Agriculture (OBERS). Since some of these projections are 
subject to considerable variation because of world events, a range of price 
scenarios are included to reveal the sensitivity of the forecast to world 
energy availability. 


4/ The American Gas Association data are confidential. 
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Elasticity 


Household demand for electricity in cities throughout the United States 
depends on the interaction of many factors. Climate, life-style, population 
and prices are all important determinants of demand. Simple models, formerly 
popular with utilities, such as the one that projects trends in the ratio of 
energy to population, cannot hope to track reality in recent years as prices 
have varied. No simple constant proportion exists by which one can multiply 
population to obtain aggregate future energy demand. 


The accumulated evidence supports the assertion that prices are an 
important, though not the only, determinant of energy demand. This body of 
evidence is already substantial and is documented in Chapter II. Although 
this report asserts the importance of prices other factors cannot be 
neglected. The fact that buyers tend to purchase less electricity as prices 
rise, other things held constant, does not deny the effect of climatic 
differences among regions on relative consumption rates. Expectations for 
further increases in the volatility of prices, however, requires careful 
analysis of price effects. Several general conclusions can be drawn from the 
already accumulated evidence and these can be most succinctly stated in terms 
of price and income elasticities of demand. 


The price elasticity of demand, denoted Eiy) refers to the percentage 
change in the quantity demanded of good i in response to a one percent change 
in the price of good j. If i and j are the same good, for example, we might 
want to know how electricity sales change as electricity prices change. Then 
we have what is known as an "own-price elasticity." If 1 and j are different 
goods, for example, we might wish to know what effect an increase in natural 
gas price would have on electricity sales. Then we have what 1s known as a 
"cross elasticity." The income elasticity of demand, N;, denotes the 
percentage change in demand accompanying a one percent increase in income. 


Note that since the concept of elasticity is stated in terms of percentage 
changes in quantity demanded and percentage changes in price, the measure is 
independent of any units of measure. Thus, it is possible to directly compare 
energy studies whose various units such as Btu, kWh, gallons, tons, etc. are 
used without tedious conversion calculations. This convenient attribute makes 
elasticities the most useful available measure of demand responsiveness to 
price changes. 


Own-price elasticities are generally the most useful magnitudes and 
several benchmarks exist against which to judge price responsiveness. If the 
elasticity coefficient is zero then demand is said to be perfectly price 
inelastic, that is, price has no effect on sales; if the elasticity is -1l then 
a 1 percent increase in price will cause quantity demanded to fall by l 
percent. In this case, demand is termed "unitary elastic." If the price 
elasticity of demand is between zero and -l1 then demand is said to be “price 
inelastic." An elasticity larger than 1 in absolute value is termed "price 
elastic." Elasticities refer to demand curves and some of these demand curves 
and matching elasticities are depicted in Figure 1. The correspondence 
between price elasticities and economic terminology is shown below. 


In the case of zero price elasticity, consumers ignore price changes and 
spend whatever is required to obtain desired quantities. The opposite extreme 
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Figure 1--Own-price elasticities and the corresponding demand curves. 


iad) 
= 
— 
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Cross-elasticities and terms 
Cross-elasticity | Term 
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Own-Elasticities, Terms and Effects 
Change in Spending 
Elasticity Term When Price Increases 
0 Perfectly Inelastic Spending Rises 
Proportionately with Price 
O> Byy > -1 Inelastic Spending Rises 
=~] Unitary Elastic Spending is Constant 
~l> Ejj > -» Elastic Spending Declines 


- © Perfectly Elastic Spending Approaches Zero 
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is shown by perfect price elasticity where a very small rise in price is 
sufficient to cause consumers to substitute other fuels in place of 
electricity. Neither extreme is appropriate to describe household electricity 
demand. 


Only negative own-elasticities have been considered because economic 
theory suggests that, as prices rise, quantity sold falls (other things held 
constant). If Ej; > 0, the response can be said to be perverse. 


While own-elasticities should always be negative, cross-elasticities can 
take on any value. If a cross-elasticity is positive the two goods can be 
said to be substitutes. For example, electricity and natural gas are 
substitutes if, as the price of electricity increases, the demand for natural 
gas increases.5/ If two goods are complements it follows that Ejj < 0. 


Energy inputs differ from many other consumer goods because demand 
elasticities differ between the short run and the long run. Energy is a 
product which is not consumed directly. For example, space heat is the 
product the consumer enjoys rather than electricity. Because electricity is 
an intermediate good, it is consumed in conjunction with other inputs 
including capital stocks. For the consumer, this means refrigerators and 
stoves and other appliances. As long as a given capital stock remains in 
place, the possibility of substitution is limited due to the technical 
features of the refrigerators, stoves, or machines. Each piece of capital 
requires a given fuel type and has fixed efficiency. Therefore, in the short 
run, substitution of cheaper fuel types is difficult though not always 
impossible. Thus, household energy demands are generally inelastic in the 
short run and more elastic in the long run. 


In his comprehensive survey of demand price elasticities for energy 
Edmonds (30) concludes 


"While agreement is unanimous that elasticities are nonzero, and 
correct (negative) in sign, no agreement exists with regard to their 
actual numerical values. Estimates sometimes diverge by an entire 
order of magnitude. ... At the writing of this paper, economists 
still search for the perfect model, data base, and R2 and the 
funding to make it all possible." 


Chapter II contains a comprehensive survey of price and income elasticity 
estimates for household electricity consumption. The truth of Edmond's 
assertion is borne out in that review. 


5/This statement refers only to the consumer demand market where electricity 
and natural gas compete to produce an output like space heat. The statement 
is not likely to be true in a supply market where natural gas is a fuel used 
in the production of electricity. In the latter market the demand for natural 
gas may fall as the price of electricity increases. 
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CHAPTER II. SURVEY OF PREVIOUS WORK 


Conceptual Issues 


Existing studies of residential electricity demand are extremely 
heterogeneous in the structuring of their models and in their research 
findings. All studies agree that prices are important in projecting future 
sales of electricity. Differing methodologies employed, reflecting differing 
assumptions concerning the appropriate way to describe the behavior of 
purchasers and sellers of electricity, result in quite different estimates of 
the effects of price changes on quantity sold.6 


In the process of replicating some of the existing studies uSing our 
improved data set, new insights suggest certain modifications in the 
structuring of model specifications. Also, combining the ideas of previous 
researchers has resulted in new models to be investigated. 


We have determined empirically that long-run electricity demand 
projections should consider changes in household share of appliances powered 
by electricity. However, our test of this hypothesis depends on a single year 
(1960) of historical substitute fuel price data used in conjunction with 1975 
data. A complete and accurate statistical analysis of the lagged and 
cumulative effects of historical fuel prices on current and future electric 
appliance stock should be possible using the retail energy price data for 
alternative fuels by city which have been supplied to us by the American Gas 
Association. Unfortunately, the confidential data were received too late to 
undertake a dynamic modeling effort in this study. 


Several controversies concerning the specification of the model exist 
among researchers doing research on electrical demand. Some of the 
mostimportant issues are described below. 


Marginal versus average price. Early researchers (including several 
well-known econometricians) argued that consumers are unaware of block rate 
prices and, therefore, it was appropriate to use average price as an indicator 
of the total costs associated with electricity use. More recently, 
prestigious research organizations, such as the Electric Power Research 
Institute (EPRI) (95), have published models using marginal prices. Many 
electric utility officials remain skeptical of marginal prices as a 
determinant of consumer behavior. 


Simultaneity. Some researchers have argued for a simultaneous supply and 
demand model while others have either ignored or minimized the importance of a 
two-equation model. Most recently, EPRI contractors// have developed a 
model which attempts to take account of the simultaneous nature of the 
electricity market by picking the prices associated with various quantities 


6/ Contrary to past history, electricity prices are expected to be among 


the more volatile determinants of electricity sales in the future. 


7/ Taylor, et. al., see tables II-l and II-2 and the references following 
the tables. 
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sold off the utility rate schedules. Unfortunately, this technique fails to 
recognize the fact that rate schedules vary systematically among regions as a 
function of output and that this variation is an important component in the 
total explanation of the electricity market, i.e., to take rate schedules as 
being given and fixed is to ignore the supply and cost factors underlying 
them. This will be discussed in detail later in this report. 


Historical fuel prices. A third difference among researchers results from 
the fact that while the discussion of appliance stocks abounds in the 
literature, little has been done to assess the effects of past fuel prices on 
the present and future shares of electric appliances. The recent EPRI study 
(96) is one of the few to use this approach but they are forced to rely on 
distributed-lag-model techniques due to their limited historical data set and 
the very weak price data used for alternative fuels. Their research and ours, 
even with the limited data processed so far from the AGA, shows the importance 
of analyzing the influence of past prices of electricity and substitute fuels 
on the share of electrical appliances. Shifts in the appliance fuel mix is 
the primary way in which changes in relative fuel prices is translated into 
changes in quantity demanded. 


Simultaneity 


Several investigators into the demand for residential electricity have 
noted that the use of ex-post average price (total receipts divided by the 
quantity sold) not only ignores the marginal pricing structure which the 
consumer faces but also neglects to account for the simultaneity of supply and 
demand. Discussions of simultaneity, however, show a lack of agreement among 
researchers as to the causes of, and solution for, the simultaneity problem. 
Several investigators utilize ex-ante rate schedules in some manner to pick 
off the prices charged for various amounts of electricity thus claiming to 
have eliminated the simultaneity problem. A few investigators attempt , 8. 
with little success, to include a supply function based upon factors which 
presumably would affect costs of production. Neither approach seems to result 
in markedly different price elasticity estimates. This may be attributable to 
a deficiency in the structural models estimated and inadequate data. A more 
accurate depiction of the market for residential electricity may include a 
number of ingredients not found in existing studies.2/ The simultaneity 
which exists in the electricity market is made more complicated because of the 
existence of utility pricing schemes which usually amount to some form of 
decreasing rate structure (this is historically true). Most researchers have 
assumed that rate schedules are autonomously imposed by utilities and 
regulatory agencies without regard to demand or cost.10/ If this assumption 
were true it would be possible to treat prices shown by actual rate schedules 
as exogenous. Unfortunately, this technique is not sufficient. In order to 
identify demand, the residential electricity demand model must take account of 
factors which affect cost and efficiency if (as seems extremely likely) the 


8/ See Table 3. Halvorsen, Levy, Asbury, and Griffin attempted two stage 
least squares estimation of simultaneous equation models. 

9/ One ingredient concerns the use of appliance stocks in a manner 
appropriate to measure the total sensitivity of quantity sold to electricity 
price. This case will be discussed later. 

10/ For example, Chern states, "Simultaneity exists in an average price 
model because, with declining block rates, the average price depends on the 
quantity consumed even though rates are determined exogenously" (20). 
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cost-efficiency factors are related to output (65). By using only the given 
rate schedules to find average and/or marginal prices researchers have 
ignoredthe influence of the cost-output relation on the setting of rates. 
Utilities are not subSidized and thus are constrained to set prices so that 
receipts cover costs. If costs are a function of output then simply picking 
prices from given rate schedules fails to account for the influence of output 
on cost efficiency and therefore on the average price required to cover 
costs. Omitting the output-cost effects on the selection of a rate schedule 
may tend to falsely attribute to price elasticity of demand those effects 
caused by cost-efficiency changes due to changing output levels. 


To assume that rate schedules are set autonomously without regard to 
production costs is to disregard the basic cost-output relationship 
conventionally used in the supply component of supply-demand analysis. While 
it is true that some characteristics of the rate schedule could be set 
autonomously, it is also true that utility districts need to cover costs of 
operation from receipts. This implies that, while the nature of the rate 
schedule (such as declining block rates) can be autonomously set, the average 
price must be a function of output if cost-efficiency is a function of 
output. For example, a utility and its regulators could choose to make a rate 
schedule decline more rapidly but they might have to simultaneously shift the 
schedule upward in order to continue to cover costs of production (depending 
on the elasticity of demand and on whether steepening the rate schedule 
increased the discriminatory pricing extraction of consumer surplus) 

(Figure 2). 


An estimation technigue which allows for the discretionary power of 
utilities to set the height and slope of rate schedules, but which takes 
account of supply (cost efficiency as a function of output) factors 
determining the required average price (to cover costs for a given output 
level), has not been achieved. One possibility might be to use the technigue 
of picking marginal prices from ex-ante rate schedules and also estimating, 
via a simultaneous model, the average price as required by supply 
considerations. Both the exogenously determined ex-ante marginal price and 
the simultaneously determined average price would be included when estimating 
the simultaneous model. As far aS we are aware, no Such technique has ever 
been attempted or even suggested. 


Long-Run Demand 


The second ingredient which we feel is necessary to accurately describe 
the residential electricity market concerns appliance stocks. Several authors 
(Taylor, for example) have argued that the demand for appliance stocks is the 
long-run demand for electricity. As noted by Cooper(27), this argument is 
incorrect. To say that the factors which determine appliance stocks totally 
measure the sensitivity of electricity sales to price changes is to assume 
implicitly that appliance stocks must be used at a constant rate. Historical 
data refute this assumption. 


Electricity sales are determined both by saturation levels of electric 
appliance ownership and by the rate of utilization of these appliances. A 
pure appliance stock model which explains only electric appliance ownership 
saturation levels neglects the effect of prices on rates of utilization. 
Surely, as prices of electricity rise, a potential consumer may desire to 
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Figure 2--The identification problem for residential 
electricity demand 


(Assume the observational unit is the electric utility) 


Price Case 1: Demand is traced 
when the rate 
schedule fluctuates 
more than demand. 

Quantity Per Time Period 
Price Case 2: Rate Schedule is 


traced when demand 
fluctuates more than 
the rate schedule. 


Mongrel Equation is 
traced when both the 
rate schedule and 
demand fluctuate. 


Quantity Per Time Period 
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maintain ownership of certain electric appliances and use them less. Electric 
heat might be supplemented by more insulation and greater use of a wood 
burning fireplace or stove. Cooler water might be used in the washer and for 
other household purposes thus reducing the utilization rate of electric water 
heaters. Many other examples could be made to demonstrate this point. 


The Almon Lag Technique 


One of the earlier investigators of residential electricity demand (35) 
reject a stock adjustment model (distributed lag) for use in explaining the 
demand for appliances. 


",...our analysis is concerned with the stock of certain physical units; 
units moreover, which have the property that a given household generally 
owns one of them at the most. Net changes in that stock, therefore, come 
about overwhelmingly by the purchase of new units by households not 
previously owning one. While it is true...that the demonstration effects 
of other household's possessions may influence such purchases, it is 
unreasonable to use a model which makes such purchases proportional to the 
difference between desired and actual stock. The actual stock that counts 
here is zero and the fact that other households in the economic aggregate 
being considered own the good already is not relevant in this place." 


The Fisher-Kaysen argument against a stock adjustment model (Koyck model) 
to predict appliances certainly seems convincing. However, the stock 
adjustment rationale is only one of several models leading to a flow 
adjustment, distributed lag formulation. We suggest that appliance stocks may 
be purchased on the basis of an expectations model leading to a distributed 
lag effect. Assume that (among other things) the stock of electric appliances 
is based upon the expected price for electricity vis-a-vis other fuels (for 
appliances where other fuels are feasible). Expected prices are not 
observable directly but a second equation which formulates how price 
expectations are formed in relation to actual historical prices can also be 
hypothesized. Either the Nerlove model of expectations (described elsewhere) 
or a simple assumption that expected price is based upon a geometrically 
declining weighted average of past pricesll/ will result in a distributed 
lag model. 


One of our exploratory investigations on electric appliance demand 
utilized the “short-cut” distributed lag model in which lagged values of 
appliance stocks are used to explain current stocks. At the time of this 
study, the complete historical data set had not been received from AGA and 
only 1960 appliance stocks data were available with which to form a 
Gistributed lag model to predict 1970 electric appliance saturation. Future 
modeling could include weighted averages of alternative fuel prices and 
electric fuel prices from 1960 to 1975 as predictors of current electric 
appliance saturation, i.e., the Almon distributed-lag method which reduces 
problems associated with serial correlation. Expected prices are incorporated 
directly into the model in place of the “short-cut” technique which utilized 
previous saturation rates. 


ony, This weighting scheme is required in order to achieve a reduced form 
equation containing observable variables if the "short-cut" technique (where 
the dependent variable is regressed on its previous value) is used. The 
shortcomings of this technique are discussed at length in a later chapter. 
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Previous Major Findings 


The salient features of 20 studies conducted between 1962 and 1977 are 
presented in Table 3 and Table 4. Table 3 shows the estimated price 
elasticity, data type, and estimation method while Table 4 shows the cross-— 


elasticity to natural gas, the income elasticity, and describes the variables 
used in the model. 


The difference between the price elasticity estimates arises partly 
because electrical prices are defined differently between studies. Studies l, 
2, 4, 6, 7, 8, 10, 14, 16, 17, 18, and 19 use an average price of electricity 
defined as revenue per kWh. Studies 3, 5, and 13 use marginal prices defined 
as the FERC's typical bill for a certain number of kilowatt hours per period. 
Study 20 uses actual rate schedules to determine marginal prices. Study 9 
uses the implicit price deflator for electricity estimated by the U.S. 
Department of Commerce. 


There are conceptual difficulties associated with using either an average 
Or a marginal price in the demand function. The problem with uSing average 
price in the demand function is that average price includes demand in its 
definition. Studies 6, 7, 10, 12, and 16 attempt to get around this problem 
by simultaneously estimating price and quantity demanded. In studies 6 and 16 
two average specifications for the price equations are tried, one with average 
price determined by cost variables (fuels, etc.) and the other with price 
explained by typical electric bills at three levels of usage. It is 
interesting to note that similar price elasticity estimates are obtained with 
both methods. The problem with including only the marginal price in the 
demand function is that the price elasticity will change if customers jump to 
new rate schedules over the period in which the price elasticities are | 
estimated. It can be shown that using only a marginal tariff rate can produce 
an upward bias in the price elasticity estimate. Study 13 attempts to 
circumvent this problem by estimating three different demand functions, each 
using a marginal price for a different rate class. Taylor argues that both 
average and marginal prices should be included in the demand function but that 
actual rate schedules rather than expost schedules should be used (95) .12/ 
It does not appear that the theory of price determination in electricity 


12/ If Taylor is correct, the analysis of underspecification in Chapter III 
will show the bias in the elasticity estimate. Let P, stand for marginal 


price and P, stand for average price. The primes on each variable indicate 
a logarithmic transformation. If the true relation is 


, ah P f 
Q' = 8, +8 Pins + 83 hace On 


and P’,, is excluded, then the estimate for 82 will be 
Bou eons 3 m a 
where Bo < 0, 83 < 0, Var (Pm) > 0 and, with declining rate schedules, 


Covar (P Pay > og. Thus B82 is more strongly negative than is the true 
elasticity with respect to marginal price (82). 


10. 
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Table 3--Selected estimates of long-run price elasticity 
of residential electricity demand* 


Authors and Date** 


Fisher-Kaysen, 1962 


Houthakker-Taylor, 1970 


Wilson, 1971 


Anderson, 1972 


Halverson, 1973 


Halverson, 1973 


Levy, 1973 


Mount, 1973, 
with Chapman 
and Tyrrell 


Verleger, 1973 


Asbury, 1974 


Price Elasticity 


-0.22 to -0.55 


=1i,09 


3 


-0.91 and -0.85 


solver ee 


-1.00 


alee 


-1.20 


=N62)5) 


19597) 0. 99 
965), 0.97, 
108-~)—1.03 


Data Type and Estimation Technique 


1946-49 and 1951-57 pooled 
cross section and time series of 
states. 


1947-64 time Series, a state ad- 
justment model, ordinary least 
squares technique 


1966 cross section of states, 
ordinary least squares technique 


1969 cross section of 77 cities, 
ordinary least squares technique 


1969 and 1970 cross section of 

states, ordinary least squares 

technique (also tries appliance 
stock demand elasticity -.84) 


1969 cross section of states, 
ordinary least squares, similar 
results with 2SLS technique 
(supply variables generally not 
significant ) 


1970 cross section of 67 utilities 
in New England, 2SLS, similar re 
sults from ordinary least squares 
technique 


1947-70 pooled cross section and 
time series of states, variable 
elasticity (by state) geometric 
lag model, ordinary least squares, 
similar results with instrumental 
variables technique 


Time-series, ordinary least 
squares, includes commercial 
customers 


1959, 1965, and 1970 cross-sec- 
ton of states, ordinary least 
squares technique and 2SLS 
(elasticity -0.88 to -1.14) 


Continued 


Table 3--Continued. 
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Authors and Date 


Price Elasticity 


Data Type and Estimation Technique 


ale 


12. 


13. 


14. 


ie 


16. 


Ii. 


18. 


Gill, 1974, with 
Chern, Ellison, and 
Tyrrell 


Griffin, 1974 


Houthakker, 1974, 
with Verleger and 
Sheehan 


Randall, 1974, 
with Ives and Ryan 


Gill, 1975, with 
Tyrrell, Ellison, 
Chern, Haynes, and 
Kaplan 


Halvorsen, 1975 


Lyman, 1975 


Uri, 1975 


-0.85 


-0.52 


100-250kWh, -1.2 
200-500kWh, -1.02 
250-500kWh, -0.44 


-1.06 


eS 


calle dis) 


West coast -1.02 
Northwest -1.15 
N. Midwest-0.13 


Midwest -1.05 
MidAtlantic-0.82 
S. West -1.17 


N. Texas ~1.13 
S. Texas -0.58 


South -1.36 
Florida & 
-] ® $2 


1962-1972 pooled cross section 
and time series of distributors in 
TVA, geometric lag model, general- 
ized least squares technique 


1970 U.S. time series data, Almon 
polynomial distributed lag model, 
2SLS technique, multiequation 
model 


1961-71 pooled cross section and 
time series of states, flow ad- 
justment model, generalized least 
squares technique 


1970 cross section, Southwest 
U.S., ordinary least squares tech- 
nique 


See Gill above 


1961-69 pooled cross section and 
time series of states, 2SLS 
technique 


1959-68, pooled cross section and 
time series separated into 10 
regions, data on 67 utilities, 
author says price elasticity esti- 
mates lie between the short run 
and the long run 


1947-70 pooled cross section and 
time series of states, variable 
elasticity geometric lag model, 
ordinary least squares technique 


Continued 
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Table 3--Continued. 


Authors and Date Price Elasticity Data Type and Estimation Technique 

19. Hassan, 1977, Uncompensated 1947-72 Canadian time series data, 
with Johnson and Model demand for rent, fuel, and elec- 
Green Rotterdam-- -0.74 tric power 


Linear-- -0.86 
Log-Log-— -0.91 


Compensated Model 
Rotterdam-- -0.54 


Linear-- —0.62 
Log—Log-- -0.79 


State Adjustment Model 
76 


Dynamic Linear Expenditure Model 


Uncompensated -0.62 
Compensated -0.37 


20. Taylor, 1977, Log-Log Flow 1956-72 pooled cross section and 
with Blattenberger and Adjustment Model time series of states, the author 
Verleger -0.82 _ declares the larger elasticity 

Linear Flow estimates “Implausible," however, 
Adjustment Model the R2 is .99 for all distri- 
—6.21 buted lag models, estimated by 


Log-Log Koyck Dis- by variance components technique 


tributed Lag Model 
-0. 81 


Linear Koyck Dis- 
tributed Lag Model 


-3.68 
Appliance Stock 
Model 
-0.46 to -.90 
Zia Yang, 1976 -.57 to -1.119 1962-75 pooled cross section and 
time series of states, comparisons 
of cross sections for different 
periods are made. 
227 UGE, LOTT -.35 Log-linear flow adjustment model 


money illusion is assumed so that 
deflation is not attempted, 
Zzellner's seemingly unrelated re- 
gression method is used on ag- 
gregate U.S. time series data from 
1971-1976 


*Nominally Presented by the Authors as Long Run 


**One of the earliest comparable studies is by Houthakker in 1951. He finds a price 


elasticity with fixed stocks of applicances of -.89. 
for 42 British towns in 1937 and 1938. 


His data are cross section 
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Table 4--Selected estimates of long-run elasticity of natural gas price 
and income* | 


Cross Elasticity Income Dependent Variable and Other 
Authors with Natural Gas Elasticity Variables in the Model 


1. Fisher-Kaysen Appliance stocks;income 
population, number of wired 
households (uses average 
revenues for price) 


2. Houthakker-Taylor 0.00 1.93 Per capita consumption (use 
average revenues for price) 


3. Wilson 0.31 —0.46 KWh/household, house size, 
degree days (uses TEB 500 
kWh for price) 


demographic and climatic 
variables (hard to inter- 
pret Since not tied to elec 
tricity sales) 


4. Anderson O2t3re: S20 1.13 & 0.94 Appliance fuel choice; 


| 


5. Anderson 0.30 0.80 KWh/household;: (alternate 
fuel prices and climatic 
variables play minor roles) 
(uses TEB 500 to 1,000 kWh 
for price) 


6. Halverson 0.16 0.70 Avg. use and average 
revenue (simultaneous); incq 
climate, demographic variab 
alternate fuel prices. 


7. Levy Ook 0.44 Similar to Halverson's stud 
Uses utility service areas J 
Place of state data, price, | 
income, utility ownership, 
family size 


8. Mount, Chapman, 0.19 0.20 Total kWh consumption; 
and Tyrrell : per capita income, gas 
and appliance prices, popu- ® 
lation, regional dummies 
(Average revenues for price} 
9. Verleger Not estimated KWh/capita; price 
data is the implicit price 
deflator for electricity 


Continued 


Table 4--Continued 


Cross Elasticity Income Dependent Variable and Other 

Authors with Natural Gas Elasticity Variables in the Model 

10. Asbury 195977 O19 0.43 Population density (Average 

1965, 0.14 0.40 revenues for price) 
1970 70.17 0.18 
11. Gill, Chern, 0.66 0.125 Total kWh consumption, number 
Ellison, & of customers, appliance 
Tyrrell prices 
12. "Griffin Not estimated 0.88 KWh/capita, stocks of 
electricity consuming equip- 
ment; (part of a 25-equation 
macroeconomic model) (appliance 
stocks not related to elec- 
tricity price) 

13. Houthakker, Not estimated 100-250kWh, KWh/capita (marginal 
Verleger, & 1.60 prices using TEB data) 
Sheehan 200-500kWh, 

1.64 
250-500kWh, 
2.20 

14. Randall, Ives, 0.25 0.27 Total kWh consumption, 

& Ryan number of customers (uses aver- 
age price) 

15. Gill, Tyrrell, 0.53 0.23 See Gill above. 

Ellison, Haynes, 
& Kaplan 
16. Halvorsen 0.04 0.52. KWh/household; % rural 
population, household 
size (cost variables in supply 
curve) 

17. Lyman West coast 0.33 0.17 KWh/customer; index of 

Northwest -0.28 0.55 other prices, after tax 

N. Midwest 0.36 0.24 income/number of households, 
Midwest 0.03 0.35 mean household size, median 

MidAtlantic0.04 0.30 household income and propor- 
S. West -0.02 0.08 tion of households by income 
N. Texas 0.29 £-2.07 range, average temperature, 

S. Texas -1.27 -2.14 average degree days, average 
South -0.80 -0.54 relative humidity, variation 
Florida & in climate variables, number 
Gulf 0.49 Oo of customers in utility 


service area, population 
density, cities/land area 
interaction terms (uses 
average revenue for price) 


Continued 
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Table 4--Continued 


Income 
Elasticity 


Cross Elasticity 
Authors with Natural Gas 


IWS}. Wiest 0.42 


19. Rotterdam 
model 
1.04 

Linear model 
1.30 

Log-Log model 
0.65 


Hassan, Not estimated 
Johnson, & 


Green 


20. Taylor 0.00 Log-Log Flow 


Adjustment Model 
1.08 


0.00 Linear Flow 


Adjustment Model 
8.56 


0.00 Log-Log Koyck 
Distributed 


Lag Model 
1.05 


0.00 Linear Koyck 
Distributed 


Lag Model 
3.84 


0.00 Appliance Stock 
Model 
EO; Eorr.36 


21. Yang -05 to .158 =—/20) tOr.28 


Electric Bills, average gas 


Dependent Vari1labie and Other 
Variables in the Model] 


Total kWh consumption 
(uses average revenue for 
price) 


Per capita, spending J 
for rent, fuel, and electric 1 


power (quantity per capita i 
some models) 


KWh/customer, (Marginal 
prices using ex-ante rate 
schedules. ) 


Average per capita quantity 
of electricity, marginal ) 
electricity price for 250-5 
250-500 kWh from Typical 1 


prices from Gas Facts, 
average per capita income. 
All monetary data are 
deflated. Average electri-g 


city price is also included f 


Continued 
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Table 4——Continued 


Cross Elasticity Income Dependent Variable and Other 
Authors with Natural Gas Elasticity Variables in the Model 
POE) al Not estimated 2.00 Monthly sales per customer; 


average electricity price, 
monthly average of average 
weekly earnings, heating 
degree days, oil-embargo 
dummy variable. 


*Nominally Presented by the Authors as Long Run 
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demand has been adequately worked out. It may be that these theoretical 
problems are not very important when applied to a broad aggregate of electric 
utility customers. 

Methodological differences among the studies also produce variations in 
the short run (SR) and long run (LR) price elasticity estimates. Differences 
between SR and LR elasticity estimates are explained in part by the factors 
discussed above but also by the fact that different distributed lags are 
used. For those studies that utilize either time series or time series and 
pooled cross section data studies (1, 2, 8, 9, 11, 12, 13, 15, 16, 17, 18, and 
19) LR price elasticities are derived functionally from the SR elasticities 
and the distributed lag-parameter estimates. The most common method for 
estimating the distributed lag coefficients is to use the Koyck transformation 
which assumes a geometrically declining set of weights in the lag 
distribution. An example of the Koyck distributed lag is given by 


BE’. = 6a 8, re ene g 
ie POP ee ee ee a (1) 
where DE}; = electric power demand per customer in the it economic 

sector in period t, 

Kit = all variables other than the price of electricity that 
explain the demand for electric power (e.g., per capita 
income, price of gas, etc.), 

PE}t = price of electricity, and 

B9,81,82,83 = parameter estimates. 


The primes in equation (1) indicate that all variables are measured in 
logarithms so the parameter estimate 82 is a direct measure of short run 
price elasticity. To obtain the long-run price elasticity, we set 
DE~;,_17DE., 13/ and resolve for DE~., thus obtaining the LR demand 


equation 
DE’, = Bo + By Mae Fea PE at Ee (2) 


The LR price elasticity estimate is then given as 


13/ Quantity demanded is assumed constant when market equilibrium is 
reached, thus s 
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Where 83 is the distributed lag coefficient estimate. Studies 8, 9, 13, and 


16 estimate the distributed lag using the "short cut" technique described 
above. 


Although the Koyck method of estimating the distributed lag is the most 
common method used (largely because of its small data requirements) it is well 
known in the econometric literature that there are severe methodological 
problems associated with the technique. First, the geometric lag imposes 
exactly the same unrealistic distributed lag on all the independent variables 
in the equation. Second, if the random disturbances are autocorrelated bias 
results in the estimated distributed lag coefficient. Study 12 attempts to 
get around these problems by using the Almon polynomial distributed lag 
procedure. The basic idea of this method is that the lag distribution can be 
expressed as a fairly low order polynomial function of the lag itself. Thus, 
if the distributed lag model is written as 


t S 
DE’. Sheg a eu Ai ag 7 Bo x W, PE Sata it (3) 


and 


= 2 s-1 
We SO o> a, t+ 0 eee le ah OL t (4) 


then the distributed lag in equation (4) can be rewritten as weighted averages 
of PEjt-n. The Almon lag places fewer constraints on the lag structure than 
the Koyck method and eliminates the autocorrelation bias problem. Study 12 
uses the Almon lag in estimating power demand and has considerably lower SR 
and LR price elasticity estimates than study 8 which uses the Koyck lag 
distribution. On the other hand study 9, which also uses the Koyck lag 
distribution, produces SR and LR estimates similar in magnitude to those found 
in study 12. 


Time Series Versus Cross Section Estimates 


Over three-fifths of previous studies are either time series or pooled 
cross section and time series. The pooled and pure time series models invari- 
ably use the distributed lag technique to model the demand function. The 
adoption of a distributed lag technique allows the use of short run time series 
data to estimate both short run and longer run price elasticities. Long run 
elasticities are appropriate if long term predictions of electricity usage are 
desired. The time series and combination cross section and time series data 
sets (pooled) also incorporate distributed lag relationships in order to 
remove serial correlation (which can reduce efficiency and which violates the 
assumptions that the observations are independent). Unfortunately, the use of 
time series data presents several problems including contemporaneous trends in 
several of the variables (notably, prices and quantities as well as income and 
population). It is impossible to say to what degree the relation among these 
variables might be spurious (each one increasing due to exogenous influences) 
rather than reflecting a true interdependence. The distributed lag technique 
is an attempt to adjust for this but, unfortunately, tests used to show that 
serial correlation has been removed are not effective when the "short cut" 
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technique is used to generate the estimated equation. Thus, time series data 
present some insoluble problems. Improper specification of the model can 
result in bias in the elasticity estimates. These problems severely restrict 
the usefulness of increasing the sample size through the incorporation of time 
series in combination with cross section data. 


Time series data also suffers from collinearity since most all nominal 
"independent" variables contain trends which makes them correlate highly to 
each other. This means that, statistically, any one variable explains for 
all, but, in fact, many causal variables exist and it is simply the existence 
of trend which prevents the proper inclusion of the appropriate independent 
variables. Cross section analysis avoids the problem of contemporaneous 
trends although some collinearity may still exist. 


It is interesting to note in Table 5 that the price elasticity estimates 
found in the studies using cross section data all lie very closely around 
-1.0. In contrast, the studies which utilize time series data show a very 
wide variation (from nearly zero to -6.21) in price elasticity. Time series 
studies also have a wide variation among their income elasticity and cross 
elasticity estimates. 
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Table 5--Comparison of cross section and time series price 


elasticity estimates 


Cross Section Estimates 


Source* Price Elasticity 
(3) 135 
(4) -0.91 and -0.85 
(5) -1.12 
(6) —1 00 
W) -1.12 
(10) -0.93, -0.97, and -1.03 
(14) ~1.06 
(21) =Lh9 
(22) SHSS 
Time Series Estimates 
(1) 0.22 to -0.55 
(2) ala 
(8) -1.20 
(9) -0.55 
Cit} -0.85 
(12) -0.52 
(13) -0.44 to -1.2 
(15) =1 5 
(16) -1.15 
(17) =0.58" to =I. 36 
(18) -1.92 
(19) ~0.37 to -0.91 
(20) -0.46 to -6.21 


*Authors are referenced in Tables 3 and 4. 


CHAPTER III. RESEARCH METHODOLOGY 


In this section we discuss the selection of data, variables, functional 
form and mode] structure used in this research effort. 


The Single Equation Model 


Three alternatives exist to define the data included in the sample used to 
estimate the forecasting equation for household electricity consumption. 
These are the pure time series sample, the pure cross section sample, and the 
pooled cross section and time series sample. Many of the studies described in 
Chapter Two were of the pooled time series and cross section variety. For the 
purposes of this study it was decided to avoid the use of time series data 
altogether. In principle, time series data afford the opportunity to estimate 
dynamic relationships which measure the speed of reactions such as might be 
shown by the time path of household electricity sales determined by rapid 
changes in prices. Unfortunately, the ability of statistical models to reveal 
true relationships depends partly on prior knowledge of relationships. Very 
little is known about dynamic adjustments in most economic or market 
Situations and less is known about the dynamics of residential electricity 
consumption.14/ Estimation of a hypothesized dynamic model from time series 
data could serve a useful purpose in enlarging our understanding of the market 
but it would be deceptive to declare that such a model was more than 
"exploratory" testing. We have indeed engaged in some "exploratory" testing 
and it is reported in the appendix. For the purposes of a forecast it was 
decided to utilize the relatively simple cross section sample which is more 
consistent with the usual comparative static microeconomic framework of 
analysis, i.e., weekly or monthly changes are not predicted but rather "long 
run" adjustments in 1975 are estimated. 


14/ Several other comments on the use of time series data are necessary. 
Most economic data contain time trends, and time series analysis is always 
plagued with the question of whether true causation is shown by time series 
regressions or whether simple correlation of trends is involved. Spurious 
correlation of trends may account for the fact that high correlations are 
almost assured when time series data is used. Adjustment to remove Serial 
correlation, when done mechanically, may "throw out the baby with the bath 
water” while adjustments based on prior knowledge of the dynamic properties of 
the market are impossible given our limited knowledge. Serial correlation, 
common in market data also results in overstatements of the goodness of fit 
calculated by the usual multiple regression techniques which assume that no 
serial correlation exists. 

15/ Under the assumption that both appliance stocks and rates of use are at 
equilibrium levels in each regional observation rationing or abrupt price 
changes would invalidate more recent cross section data. 

16/ Many of the previous studies have selected the distributed lag 
technique (using the data saving "Short cut” technique where the dependent 
variable in time period t is regressed on the value of the dependent variable 
in time period t-l) and thus it is necessary to explain briefly why this 
approach was discarded. 
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Several alternative explanations of consumer behavior are consistent with 
the ubiquitous "short cut" version of the distributed lag model. Thus the 
model is not powerful in discriminating among alternative hypotheses but very 
often provides a "good" statistical fit to the data. The "short cut" version 
is popular among economists because it allows estimation of the lagged effects 
of causal variables on electricity consumption without requiring measurement 
of the lagged observations on the causal variables. Unfortunately, the 
requirement for geometrically declining weights on past values of causal 
variables (implied by the structure of the model) may not be appropriate to 
explain household electricity demand. Also, the lack of power to distinguish 
among alternative explanations of consumer behavior makes the technique more 
one of pragmatic estimation than one of hypothesis testing. Despite its 
acknowledged statistical weaknesses 17/ the technique is currently used by 
many who seek to explain the behavior of household consumers of electricity. 


The following discussion shows three rationalizations which are consistent 
with the geometrically declining weights in the distributed lag model. It 
should be kept in mind that these comments and later discussion pertain only 
to the "short cut" data saving model which replaces actual lagged observations 
on causal variables with the lagged observation on the dependent variable. A 
model in which past values of the causal variables are included does not 
suffer any restriction concerning weights applied to past values since the 
weights can be determined by the regression itself. A paucity of historical 
data may partly account for the lack of use of the less restrictive model 
incorporating actual lagged values of the independent variables. 


A special form of distributed lag with geometrically declining weights is 
often used to estimate long run demand for electricity. This form can be used 
to measure the lagged effects of prices or income on current electricity 
consumption without actual lagged data for the causal variables. Let EC+z 
stand for current electricity consumption. The distributed lag reduced form 
is 


eS a TRE. (5) 


where X can be measures of electricity price, prices of closely related 
goods, income or other causal variables. This type of distributed lag model 
can be interpreted in several ways. Koyck (69) has shown that, given a 
geometric lag distribution (simplified by including only one exogenous 
variable summed over 7 time periods), i.e. 


zi os 
EC, = 0 : (I-v)vX, + Uy (6) 


1/ Inclusion of a lagged value of the dependent variable as an explanatory 
variable causes added statistical difficulties; this will be discussed later. 
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that an equivalent algebraic expression to equation 6 is 


EC =e + a (1-V)X, + vu, - Us 4 (7) 


iB il 


Since equation 7 is operationally equivalent to equation 5 the Koyck 
interpretation of 5 can be expressed by equation 6. That is, electricity 
consumption would be based on past and current values of the causal variables 
with past values declining in importance by geometric progression. 


A second interpretation is provided by Cagan (16) in which consumption 
would be based on expected values of the independent variables. Since 
expectations cannot be measured directly an adaptive adjustment relation is 
assumed. Current expectations for prices or income are revised in each period 
by some fraction, k’, of the difference between previous expected values and 
actual values. The coefficient,k’, is often required to fall in the range 0 
to ] since this implies that the adjustment in expectations will not exceed 
the error in the previous forecast. The Cagan model can be summarized as 


at 3 U 
NORE Baie es (8) 


and 


oh al Ae wll + 
Ds, Xx K (X, Xx ) 


t#i t (9) 


where * denotes expected values. Solving equations 8 and 9 simultaneously 
also results in a function that is operationally equivalent to equation 5. 


Yet a third model, which is equivalent to equation 5, is due to Nerlove 
(79). Current values of the exogenous variables are assumed to determine the 
desired values of electricity consumption. However, adjustment of actual 
consumption to the desired level is not accomplished in a Single period 
(year). Hence, the adjustment is some fraction of the difference between 
actual and desired consumption. The Nerlove version can be expressed as 


EC,** = 0 Sa tee (10) 


and 


WC = BC, - BC = D(EC *4-FC 
AEC = EC, EC. en BC, 4? (11) 


where ** denotes desired values. 
Thus the distributed lag formulation might be justified by three different 
arguments: (1) consumption is determined on the basis of a weighted average 
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of past prices and income, (2) consumers purchase electricity on the basis of 
expected prices and income with an adjustment mechanism for expectations, or 
(3) households purchase electricity on the basis of current prices and income 
with a particular type of adjustment lag in setting their rate of 
consumption. Only the latter rationalization of the model has the 
characteristic that the reduced or operational form of the model does not 
exhibit serial correlation if there was none to start with. Conversely, if 
prior residuals should contain serial correlation, the latter model would 
allow no chance of improving the situation. 


Least Squares Estimates 


We can carry out the least squares calculations without making assumptions 
about the underlying situation. However, the method will only have optimal 
properties under certain assumptions. Let us begin by assuming a linear model 
such as 


Y=a +gZ + v 


uN 
S 


where E(v2) 


iT] 
© 


and a) 


for all t (residual variance constant over time). To get an unbiased estimate 
of the variance one should, in general, divide by T - K, whgre K is the number 
of estimated coefficients. Hence an unbiased estimate of w°” is given by 


The first consideration about the least squares estimates is whether the 
coefficients bj are unbiased estimates. Let us consider bj in particular 
and, within the hypothetical framework of repeated sampling, ask "If we £1x 
the values 2],...,z7, observe the y's and calculate the coefficient bj 
which results from a given sample of v's and do this many times for the same 
Z's, in the long run is the average value of bj equal to 8 (E(b) = 8)?" It 
is clear that we do need some assumptions about the error term to answer this 
question, because the statistical properties of bj; depend on those of V. 


It is convenient to work in mean deviation form. By taking the original 
equation 


Yt = OL BZ, 2s Ve 
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and summing over T observations and dividing by T, we have 


SOR nan eta) 


which, by substracting, yields 


yo lee eae) 


where y» is the sample mean of the errors and must be distinguished from the 
theoretical mean E(v)). 


Now express bj in terms of v in order to study its statistical 
properties by substituting from the previous equation to yield 


Us nec CE 
fa, 2 2 MEIC ann” 


Multiplying out the terms in the square brackets and simplifying yields 


i Be Ct - Zz) i e a as 
cue liz = 22 prt a tz, =a) 


Now we are in a position to see whether E(b) =8. The expected value of a sum 
of terms is equal to the sum of the expected values of each term considered 
independently so 


E(b) = ; 


2 ane ae 
Bivens E[(v,-v) ee 


Now if z is fixed the expected value of a constant multiplied by a random 
variable is equal to the constant multiplied by the expected value of the 
random variable so that | 


Dies 2 
E(b) Si Bo ats (ts) E(v, - we 
t 


But E(v.)= 0 so E(v) = 0 hence every element in the summation is zero and E(b) 


= 8, If the z is a random variable itself but exogenous and therefore 
independent of v we need one extra step. Independence of two random variables 
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implies that the expected value of their product is equal to the product of 
thei: separate expected values so that 


‘ ale Z 
E(b) = 8 +s Et, =) ( ‘: 


Again E (v,)= E(v) = 0 hence E(b) = 8 and b is an unbiased estimate of Bie rin 


any Single sample b will not equal 8 and there will be a sampling error of 
b- 8 but these will average out to zero in repeated samples. 


No assumption has been made about the presence or absence of autocor- 
relation in the vy; series nor have we used the assumption about the variance 
of v; being constant, i.e., b is unbiased either way. Trouble arises, 
however, when Z is a lagged endogenous variable since 


: x Zia ja 
a ee fe a= Zz) ] 


can no longer be divided into a product of two expected values because the two 
factors inSide the square brackets are not independent. For a lagged 
endogenous variable to be predetermined we require it to be independent of 
current and future disturbances but we cannot require independence of past 
errors. (If 2 is the lagged dependent variable z+ = yr-1, then 2 

depends on v~-1). So v_ in the first factor is related to some 2t¢44 

which appears in the denominator of the second factor and independence no 
longer holds and we cannot prove E(b) =8. In general, there will be a bias 
E(b) - 8 #0. Therefore, bias is to be expected in the "short cut" 
distributed lag model. 


Asymptotic properties. In a number of situations it is only possible to 
derive properties of estimators as the sample size T tends to infinity. If no 
results concerning bias, variance, etc. are obtainable for finite samples then 
the limiting case is the best we can do. The coefficient b is a consistent 
estimator of 8 if the sampling distribution of b approaches a degenerate 
distribution with all the probability concentrated at 8 as the sample size T 
becomes infinite. The probability limit of b is 8, i.e. 


p limb = 8 


fe 


Consider the possibility of repeating the regression experiment for a 
given set of z values by drawing a number of samples and calculating b for 
each sample of y that results. We shall eventually obtain a smooth 
distribution of b values, the mean and variance of this sampling distribution 
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have already been found for some cases. Now repeat the process for a new, 
larger value of T and obtain another distribution of b. Consistency requires 
that, as we do this for larger and larger sample sizes, the resulting 
distribution of b collapses to the true value 8 which we can illustrate as 
follows. 


J(,b) A J( bys 


The second diagram illustrates a situation in which the estimate is 
biased in small (finite) samples so, for it to be consistent, two things must 
happen as T increases (1) the distribution must move toward 8 and (2) become 
Narrower as in the first diagram. 


This heuristic discussion rests upon two conditions which are actually 
sufficient conditions for consistency, asymptotic unbiasedness and a variance 
tending toward zero i.e., if lim E(b) = 8 and lim var(b) = 0, then plim(b) 

T+ © T> Tenge 
=I. In the simplest case of an exogenous z and independent y's b is a 


consistent estimate of 8 for all sample sizes and var (b) = w /E(zt-Z)? >oas To... 


Plims have a very useful property. The plim of a function of two random 
variables is equal to the value of the function when the separate plims of the 
two variables are substituted in, specifically ee ORT eT ee TN) 

We also have p Tim Me Ds lamb tye) 
2 


: ] B 
Zia B Lamy Ev Cz adil 
ey avi 12 Vp (ope ——_j—1_# ) 


ue UE Wo i a 
De ae z(Z, z) p lim E za (Z, = 4Z)4 


ie 


enabling us to consider the question of bias in large samples. (We divide the 
numerator and denominator by T to ensure that the denominator remains finite 
as T>~i.e., it is now the variance of z). If Zz is predetermined 
(independent of current and future disturbances) we obtain a numerator of zero 
yielding the result that b is an asymptotically unbiased estimate of 8. It is 


also possible to prove consistency. However, even these techniques do not 
help if z is lagged endogenous and v is autocorrelated, a likely condition in 
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most "short cut" distributed lag models. Thus, the "Short cut" distributed 
lag ‘odel will yield biased coefficients and it is likely that the bias will 
not. disappear even when the sample size approaches infinity. 


Having rejected the time series sample data and the dynamic (as opposed to 
static) modelling effort implied by time series data, only the pure cross 
section alternative remains. For the purpose of projecting household 
electricity consumption in certain specified utility service areas throughout 
the United States it seemed appropriate to utilize the cities which make up 
those service areas aS sample observations. Thus a cross section regression 
equation is estimated across selected cities in the United States. The 
statistical model is consistent with the comparative static framework of 
analysis whose properties are more completely understood than dynamic economic 
jnodels. 


The economic model incorporates both monetary and physical variables. We 
stipulate that electricity consumption is affected both by prices of 
electricity and closely related goods, consumer income, and also by 
nonmonetary variables such as climate or weather. The neoclassical theory of 
consumer demand, which is the basis for our model, requires that "rational" 
consumers base their decisions to allocate their incomes among competing uses 
on the basis of relative worth to the consumer at the margin and on relative 
cost. Unless there is reason to suspect that consumers are fooled by 
inflation of incomes or prices into making irrational decisions, or unless 
there is reason to think that consumers are acting speculatively (buying to 
resell at a higher price, for example), it is appropriate to assume that 
consumers examine the real values and real costs associated with their 
purchase of electricity. Thus, all prices and income should be expressed in 
real terms to remove the unwanted variation due to changes in the cost of 
living (exchange value of money for goods). Deflation of monetary variables 
has been attempted in some of the previous studies but the accuracy of the 
deflation is suspect since comparative cost of living indices are only 
available for cities and most of the pooled cross section and time series 
models or pure cross section models utilized state rather than city data. 
Deflation also proved to be a problem for some of the cities in our 
sample. 18/ The methods applied to the deflation problem in this study are 
discussed in Chapter Four. 


A final difficulty occurs with regard to the model structure. 
Conventional theory describes a market in terms of two relationships, the 
behavior of the consumer or demand and the behavior of the producers and 
sellers or supply. Some previous studies have included a supply equation in 
their model of household electricity consumption but with little apparent 
effect. Nonetheless, correct specification of the structure of a market 
should include price quantity relationships on the supply as well as the 
demand side of the market. There is some question as to the meaning of supply 
in a regulated industry whose price may not always reflect marginal production 
costs as they would in competitive markets.19/ In any event, the result of 


18/ Nonetheless, experimental regression analysis utilizing non-deflated 
monetary variables consistently resulted in much less accurate equations. 
Typically, the coefficient of determination was lowered by 20-30 percent when 
non-deflated variables were used. 

19/ Unless utilities are subsidized by general taxes, revenues must cover 
costs in the long run. 
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neglecting the supply equation is to estimate a "mongrel" demand equation 
whose price parameters may be biased through the inclusion of price as a 
causal variable when it is truly simultaneously determined along with quantity 
demanded.20/ Bias in the estimated price coefficients is due to correlaton 
between the price variable and the error term in the regression equation. The 
least squares regression technique assumes that no correlation exists and thus 
calculates a biased price coefficient. The bias is minimized, however, if the 
error term in the regression is small. That is, if the regression fit (R2) 

is very high, then the residual variance is small and the extent of the 
correlation of the price variable and error term is reduced. Bias should be 
minimized if the fit of the regression equation to the sample data is very 
good. This factor may partly explain the lack of large differences among 
studies which incorporated two equations as opposed to studies which only 
estimated a demand relationship. Every effort has been expended here to 
optimize the fit of the estimated demand equation. Efforts toward this goal 
seem to bring more immediate reward than do the efforts toward estimation of a 
simultaneous model. The latter is most important, however, if a better 
understanding of the market is the primary goal (as opposed to developing a 
forecast). 


Specification Error 


Once the correct specification of the model is assumed, model estimation 
and model testing become relatively straightforward. In reality, however, one 
can never be sure that a given model is correctly specified. In fact, applied 
researchers usually examine more than one possible model specification, 
attempting to find the specification which best describes the physical process 
under study. 


Underspecification 
Consider first the case in which a variable is omitted from a "true" or 
correct model specification. Assume that the true model is given by 


= Bos ae Seve aves 
Va Be ee See lee vue 


while the regression model is given by 


= * ° * 
ey Bo voF oF U, ' 


If all the assumptions of the classical linear model hold, the estimated slope 
parameter is 


Bo = 


Bs 


vines 
L574 


20/ Note the graph shown in Figure 2. 
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Substituting y;, as defined in the true model, into the estimator equation 
for £2 yields 
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Since there is no guarantee that the second term will be 0, the least-squares 
Slope estimate yields a biased estimate of the true slope parameter 82 . This 
bias will not disappear as the sample size grows large so that the omission of 
a variable from the true model yields inconsistent parameter estimates as 
well. The only case in which the bias (and inconsistency) will completely 
disappear occurs when COV ( x5, %,) = 0, that is, whenxz and x3 are 
uncorrelated in the sample. this result generalizes if there are numerous 
independent variables. Only when the omitted variable is uncorrelated with 
all the included independent variables does the bias disappear, which is 
extremely unlikely with data on the demand for electricity. 


The formula above is useful because it tells that the direction of any 
bias which might occur depends upon the correlation between the omitted 
variable and all included variables as well as on the sign of the true slope 
coefficient 83 . To the extent that x2 and xz are highly correlated the 
coefficient of x» will include the effect of the xz variable and will be 
biased. When zx» and x3 are uncorrelated xg picks up none of the effect 
of xz and no bias occurs. As a practical matter it is the extent of the 
specification bias which is important. This suggests that a careful 
researcher will consider not only the question of missing variables but their 
possible correlation with included model variables as well. 


To be complete, we should pause to consider the effect of variable 
omission on the error variance and its estimate. First consider the case in 
Which x5 and x3 are uncorrelated. Then 83 will be an unbiased estimator 
of 82 and will have an identical variance with 82 . The only difficulty 
with model misspecification arises because the usual estimate of the variance 
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of g5 will be biased. However, in the more general case when Zo» and x 
are correlated the two estimators will not have identical variances. In the 
two variable model, the actual variance of 85 will be less than the actual 
variance of 62 , even though the model is misspecified. If one is willing to 
give up lack of bias as an important objective then the omitted variable 
specification has some merit. In certain cases, specification bias can prove 
useful as a tool of analysis. 


The theory of specification error is germane to our study because 
variables which measure electric appliance stocks have been intentionally 
excluded from our model. It is reasonable to hypothesize that households with 
larger numbers of electric appliances will tend to purchase more electricity 
than would households with fewer appliances and/or with nonelectric 
appliances. The results of previous studies, and our own experiments 
(reported in a later chapter) verify this conjecture. Discussion in Chapter 
VI argues for the need to perform added research making use of recently 
acquired data on alternative fuel prices by city and the further development 
of some rather complex econometric models of the household electricity market 
in conjunction with appliance stocks. For the purposes of this part of the 
research, data with which to project appliance stocks into the future were 
sketchy and the models which were built and tested using the available data 
did not fit more accurately than the less complicated model presented here. 
Since the use of econometric models for making forecasts requires projecting 
future values of the causal variables it is preferable to select the model 
which requires the fewest independent variables, all other things equal. (It 
is our belief, however, that continued development of the more realistic 
models of the household electricity market will be worthwhile in the light of 
the recently acquired fuel price data.) 


The exclusion of electric appliance stocks from the electricity demand 
equation has the effect of biasing the coefficients of other variables 
included in the estimated equation if they are correlated with the excluded 
appliance stocks variables. The following examples will demonstrate the 
nature of the bias. It is shown that the exclusion of electric appliance 
stocks from the model tends to increase the direct price effects estimated in 
the model. It is argued that this is useful since electricity price would 
have influenced appliance stocks which then would have affected electricity 
sales. The exclusion of appliance stocks simply short cuts the process and 
attributes changes in electricity consumption directly to its price rather 
than distributing the effect of price changes between the change in the 
utilization rate of electric appliances and the changes in electric appliance 
stocks. 21/ 


For example consider a true demand model given by 


Q. = DOP. + DAC. + b, SH, + U. 


where Q is kWh per customer per year, P is electricity price, AC is the 
percent of households which have electric air conditioning and SH is the 


21/ We believe that rising electricity prices would cause reductions in the 
rate of utilization of electric appliances and encourage the substitution of 
nonelectric appliances. 
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percent of households with electric space heating. A regression is estimated 
in which AC and SH are intentionally excluded, i.e., 


Q. = boP. ot = 


It is a delish. shown from the previous discussion that the expected value of 
the price coefficient in the underspecified regression equation will be 


E(b,) = By + b3 Cac * Pa “st 


ZAC, as 
where Cac aS in the regression 


as in the regression 


SH. = CoP. tle 


If the signs of the Cac and CSH terms are known then the direction of the 
bias may be deduced. Cypc and Coy measure the effects of electricity price 
on the ownership of electric eonitna and heating appliances. A prtort | 
reasoning and previous research have both shown that Cac < 0 and Cou < 

It has previously been shown that b3 > 0 and bq > 0 since households with 
a larger share of electric appliances tend to purchase more electricity. Thus 
the net effect of the bias in bos introduced by exclusion of AC and SH 
from the demand model, is negative, i.e., by Cae + Dy CoO. 


Since the effect of electricity price on electricity consumption is also 
negative (b, < 0) the bias tends to increase the absolute value of the price 
coefficients. Price is attributed more direct influence over consumption than 
truly exists. The bias is innocuous, however, because the biased price 
coefficient combines both the direct effect of price on appliance utilization 
rates and the indirect effect of price on electric appliance stocks. 
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Measurement of both phenomenon is thus incorporated in the forecasting 
equation. 


The Log-Linear Functional Form 


Any equation form selected to represent a demand function must be 
considered an approximation because of limitations imposed by the regression 
technique. The basic premise, supplied by neoclassical economic theory, is 
that a negative or inverse price-quantity relationship holds. In addition to 
the criteria of the best statistical fit to the data, a particular functional 
form may be selected for ease of use and other attributes. The log-linear 
functional form is extremely popular in demand studies. In this section, we 
point out a particular advantage of this equation form which allows the use of 
easily obtained average price data (which some would argue is also the most 
appropriate) to estimate a model which is also consistent with marginal price 
data (the "correct" price measure according to economic theory). The second 
topic of this section points out an often ignored statistical problem 


associated with the use of log-transformations applied to variables used in 
regression analysis. 


Use _in Avoiding the Average-Marginal Price Controversy 


A marginal-average price discrepancy exists because electricity is not 
sold at a flat rate. The average price a consumer pays depends on the amount 
of energy purchased. Additional, or "marginal," units of the product become 
less expensive, usually at set intervals of "blocks" purchased per month. 
Consumers choose the amount to be consumed or the basis of an entire rate 
schedule of prices. This creates a kinked budget constraint or total cost 
curve for the consumer (Figure 3). It appears impossible to define a 
traditional demand schedule for an individual because there are as many prices 
as there are rates in the declining block rate structure. If one constructs a 
demand schedule for each of these rates, assuming that the others are held 
constant, the kinked budget constraint results in a discontinuous “marginal 
price” demand schedule (Figure 4). The discontinuity occurs when the consumer 
moves from one rate to the next. This problem is resolved when data are 
aggregated over many individual demand schedules. Blattenberger (13) is able 
to show that a discontinuity does not appear in the aggregated demand function 
if the individuals' utility functions are randomly distributed and if the 
population is large. The discontinuities that appear in the aggregation 
process become so small relative to the total aggregate purchase that they can 
be ignored. 


A second problem remains, however, with the use of average price in the 
elasticity demand model. Traditional economic theory assumes that consumers 
react to changes in the marginal rates, not the average price of fuel. It is 
argued that elasticities which are calculated using average price data may 
yield a distorted view of reality. It is unclear whether many consumers 
know what the marginal price is for the rate block of their consumption 
level. While average price is a result of previous consumption, as well as 
current consumption, consumers may be better able to predict average price 
than marginal price. 


The controversy is moot if one adopts the log-linear functional form. The 
log-linear form has the property that the price elasticities (and other 
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Figure 3--Budget constraint kinked by a three-rate schedule 
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Figure 4--Individual consumer's demand schedule for one rate 


with all other rates, prices and income constant 
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elasticities) are shown directly by the estimated regression slope 
coefficients. Since these coefficients are constants, it follows that the 
log-linear model assumes constant elasticities. Halvorsen (51) first showed 
that if one approximated the actual rate schedule with a smooth function, and 
if the elasticity of that function were constant, then the marginal price of 
energy would be the average price times a scalar. Halvorsen is then able to 
show that the marginal price elasticity of demand will equal the price 
elasticities calculated from average price data. The proof follows. 


Let P, stand for the average price of electricity, Pm for the marginal 
price, Q for the kWh per customer per year and 2, for the other causal 
variables in the demand function. The log-linear demand model could be stated 
as 


eae de 
Pa = ¢ Q 1,4,7. 


Total revenue is defined as 


PQ = = {9 P da 


where i = f(Q) is the rate schedule or 
(ca? Ml, za ee Pad. 


Taking the derivative with respect to Q we have 


(14b) (c) (Q?n,2%%) = P 
i m 
Or 
1 = 
(1 + b) a ne 
If the log-linear model is estimated using average prices, i.e., 
B Y 
Q. = aP TA gare en 
ys Pp 

the average price elasticity will be shown by 8 and Since RA = i Mm the 
regression equation will be aw 


In Q = antilog 5 + gn Ee pee i. 
1 
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which can be written as 


In Q = antilog me g (In Ee - In (1 +b) + ¥ i rd 
, 1 | 


where 1ln(1 + b) is a constant. 


As is well known, the addition of a constant to a variable has no effect 
on the estimation of the partial regression coefficient. Therefore, 
estimated on [lnP - 1n(l1 + b)] is identical to’ on . Use of the average 
Price does not affect the estimation of although the constant term in the 
demand function will be biased. In this manner the use of the log-linear 
functional form sidesteps the marginal-average price controversy. The 
elasticity estimates apply either to average or marginal prices. 


An Unresolved Problem 


The doublelog function form 


n 
TO Vou ey Moo Bint Ly Ri (12) 


is obtained by taking the logarithms of a multiplicative, or "Cobb-Douglas" 
type model, such as 


It is the multiplicative model (13) which holds theoretical significance but, 
Since it is nonlinear, its estimation cannot be obtained by the standard 
techniques for linear models. Thus, model (13) is transformed into 


n 
AMO ctig yb Sie bs lig a. 


for the specific purpose of easier estimation (Johnston (64)). 


Since an exact one-to-one transformation exists between these two forms, 
it may appear that the properties of the estimated coefficients in the 
logarithmically transformed model (12) will remain the same when they are 
converted back to the theoretical model (13), which specifies the 
relationships between the variables which are actually observed. It must, 
however, be pointed out that the mathematical models given above, and the 
statistical model used for estimation, are distinct. When statistical models 
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are transformed by nonlinear operations the distribution of the random 
variables involved must also be transformed. Statistical properties which 
hold in one distribution need not hold in a different distribution. | 


The statistical form of the logarithmically transformed model (12) is gen- 
erally written as 


= eis 5 aes ie 5 
IG: octets Gy ay a8 ap aate tl (14) 


This statistical specification of model (14) automatically requires that the 
original multiplicative model (13) be statistically specified as 


* 
where 25 =e 


Thus, the additive error term in model (14) becomes multiplicative in model 
(15). The appropriateness of a multiplicative error must be decided on an in- 
dividual case basis (Sprent (89) since the variance of Q; increases as the 
values of the explanatory variables increase. 


It is also well known that any desirable properties (such as unbiasedness) 
which may hold for the estimates of the parameters in the logarithmically 
transformed model (14) need not hold for the corresponding parameters in the 
model (15) of theoretic interest (Goldberger (41)). 


In order to obtain desirable properties (such as unbiasedness) in the 
transformed model (14) the expected or mean value of the error term €4 must 
be zero. It is also desirable that the expected or mean value of the stochas- 
tic theoretic model (15) equal the function of the explanatory variables 
specified in the original multiplicative model (13), i.e., that Ee; = 0. By 
utilizating the Jensen Inequality it is easily shown that these two conditions 
cannot be simultaneously met in any practical problem (see Appendix I). If 
the expected value of €+ in model (14) is zero, so that the standard estima- 
tion techniques can be ttilized, the average value of Q in model (15) will be 
greater than 


Since it will also be a function of the variance of e4 in model (14). 
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This functional form also has the inherent requirement that each of the 
explanatory variables used in the model has constant elasticity (56). If this 
functional form is used to estimate a relationship where the explanatory vari- 
ables do not have constant elasticities, the estimates of these elasticities 
will be biased and the variation of the dependent variable will not be ade- 
quately explained because the model does not adequately reflect reality. 
Therefore, even if all the relevant explanétory variables were contained in 
the model, the estimates of their elasticities could be severely biased due to 
the constraint of constant elasticity. Functional forms which require 
constant elasticity (Such as double-log) should therefore be avoided whenever 
the elasticities of the explanatory variables are known to vary considerably 
as quantity levels change. 


Since every functional form always possesses certain advantages and dis- 
advantages, it is highly advisable that a group of functional forms, each of 
which is consistent with the theory involved, be investigated. Other func- 
tional forms which are consistent with demand theory and widely used are the 
loghyperbola, semilog, linear, and hyperbola (83). 


In order to use standard estimation techniques on model (14), E(e) must be 
zero. By the Jensen Inequality22/, E(e°) in model (15) must be greater than 


or equal to e E(©)with the equality holding only in the trivial case where 
var(c) = 0. Therefore, 


since 


where eO = 1,23/ 


Selecting the "Best" Regression Equation 


Suppose we wish to establish a linear regression equation to project elec- 
tricity consumption in terms of "independent" or predictor variables Xj, 


22/ Jensen Inequality: Let x be a random variable with mean E(x) and let 


g(-) be a convex function. Then E[g(x)] > g[E(x)]. 

23/ Example:Let x be a uniformly distributed discrete random variable which 
can have values of -l1, 0, or 1, each with a probability of 
E [e*] = 1/3 aula, pygmee 1/3e! ~ 1.362053 and E(x) = 1/3(-1)+1/3(0)+1/3(1)=0. 


Therefore 


E [e”] > gb (x) 
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X9,---, Xk. We assume this is the complete set of variables from which 

the equation is to be chosen and includes any transformation, such as loga- 
rithms, squares and cross products, thought to be desirable and necessary. 

Two opposing criteria for selecting an equation are usually involved. They are 


1. to make the equation useful for predictive purposes, we should want our 
model to include as many X's as possible so that reliable fitted values 
can be determined and |. 


2. because of the costs and possible errors involved in obtaining accurate 
projections on a large number of X's we should like the equation to 
include as few X's as possible. 


The compromise between these extremes is what is usually called selecting the 
best regression eguation. There is no unique statistical procedure for doing 
this and personal judgment is a necessary part of any of the statistical 
methods discussed. In this section we describe several procedures which ap- 
pear to be in current use. To add to the confusion, they may not all 
necessarily lead to the same solution when applied to the same problem. 24/ 
The procedures are discussed (1) all possible regressions, (2) backward 
elimination, (3) forward selection, and (4) stepwise regression. 292 


All Possible Regressions 


This procedure is extremely cumbersome and is very expensive, even with 
access to a high speed computer. The procedure first requires the fitting of 
every possible regression equation which involves a constant Xg plus any 
number of the variables Xj,...,X,. Since each X; can either be, or not 
be, in the equation (two possibilities) and this is true for every Xj;,i= 
1,2,...,k (kX's), there are 2K equations. (The Xg term is always in the 
equation.) If k = 10 (a not unusually excessive number) 2k = 1024 equations 
must be examined for the second step in the procedure: Here the regressions 
are divided into sets of runs which involve p variables, p = 1,2,...,k and 
each set is ordered according to some criterion. Usually the criterion is the 
value of R* achieved by the least squares fit. The leaders in this ordering 
within each set are then selected for further examination and a decision is 
made, after viewing the R2 values, on which equation is best to use. As an 
example, consider a model with a possibility of as many as four independent 
variables. The procedure is 


naa Divide the runs into five sets 
Set A consists of the run with only the mean value 
Set B consists of the four l-variable runs 
Set C consists of all the 2-variable runs 
Set D consists of all the 3-variable runs. 
Set E consists of the run with 4 variables. 


24/7 Indeed, such is the case in some of the experimental models tested in 
this study. 

25/ Future analysis of our cross city data should include the application 
of flexible functional forms (Box-Cox curves). A recent application is shown 
in "Analysis of Flexible Engel Functions" by Blaylock and Green (14). 
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2. Order the runs within each set by the value of the square of the 
multiple correlation coefficient R2. 


3. Examine the leaders and see if there is any consistent pattern of 
variables in the leading equations in each set. For example, suppose 


Variables 
Set in equation R2 
Set B Y = £(Xq) 67.5% 
Set C tue Ee kp 97.9% 
Y = £(X1, X4) 97.2% 
Set D ¥ = £(X, Xp, X4) 98.2348 
Set E (RG Ke KY) 98.2378 


(Notice that in set C there are two leaders with practically the same R2). 
If we view these results, we see that after two variables have been introduced 
further gain in R2 is minor. Examination of the correlation matrix for the 
data reveals that (X; and X3) and (Xz and Xq4) are highly corre- 
lated. Thus, the addition of further variables, when X) and X2 or Xj 
and X4 are already in the regression equation, will remove very little of 
the unexplained variation in the response. This is clearly shown by the 
slight increase in R2 from set C to set D. The gain in R¢ from set D to 
set E is extremely small. 

What equation should be selected for further attention? One of the equa- 
tions in set C is clearly indicated but which one? If £(X], Xz) is cho- 
sen there is some inconsistency because the best single variable equation 
involves X4. For this reason, many researchers would prefer to use f£(X1, 
Xa). The examination of all possible regressions does not provide a clear- 
cut answer to the problem. Other information, such as knowledge of the char- 
acteristics of the process studied and the physical role of the X-variables 
must aS always be added to enable a decision to be made. 


In general, the analysis of all regressions is quite unwarranted. While 
it means that the statistician has "looked at all possibilities” it also means 
he has examined a number of regression equations that intelligent thought 
would reject. The amount of computer time used is wasteful and the sheer 
physical effort of examining all the computer printouts is enormous when more 
than a few variables are being examined. A selection procedure which shortens 
this task is preferable. 


The Backward Elimination Procedure 


The backward elimination method is an improvement on the "all regressions" 
method in that it attempts to permit the examination, not of all regressions 
but of only the "best" regression containing a certain number of variables. 
The basic steps in the procedure are 
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1. A regression equation containing all variables is computed. 


fh) 


The partial F-test value is calculated for every variable treated as 
though zt were the last vartable to enter the regresston equation, 


3. The lowest partial F-test value, F,, is compared with a prese- 
lected significance level Fo. 


(a) If Fy, < Fg, remove the variable X; which gave rise to Fy, 
from consideration and recompute the regression equation in the 
remaining variables. Then repeat stage 2. 


(b) If Fy, > Fo, adopt the regression equation as calculated. 


First perform the regression on all independent variables. In the 
previous example find the least squares equation Y = £(X), X9, X3, 
X4). In order to eliminate variables one must determine the contribution of 
each of the variables X,, X2, X3, and Xq4 to the regression sum of 
Squares as if each were in the Tiee position. The partial F-values provide 
measures of these contributions. Using the partial F-test, choose the 
smallest value and compare it to some critical value of F based ona 
predetermined risk. If the F value for X3 for example, is smaller than 
the critical value, X3 is removed from the equation. Next, find the least 
squares equation Y = f£(X), X2, X4). Examining this equation one sees 
that Xq4 should be removed. The procedure for this elimination is similar to 
the preceding elimination. 


We now find the least squares equation Y = £(X), Xj). Both variables 
X, and Xz are significant, regardless of position, as indicated by the 
Significant partial F's. Thus, the backward elimination selection procedure 
is terminated and yields 


Y= bg + bjX) + bX. 


This is, in general, a satisfactory procedure, especially for 
statisticians who like to see all the variables in the equation so as "not to 
miss anything." It is much more economical of computer time and manpower than 
the "all regressions” method. However, if the total input data set yields an 
X'X matrix which is ill-conditioned, i.e. singular, then this procedure may 
yield nonsense because of rounding errors. With more refined computer 
programs, this is not usually a serious problem. 


The Forward Selection Procedure. 


The backward elimination method begins with the largest regression using 
all variables and subsequently reduces the number of variables until a 
decision is reached on the equation to use. The forward selection procedure 
is an attempt to achieve a similar conclusion working in the other direction, 
i.e., to insert variables until the regression equation is satisfactory. The 
order of insertion is determined by using the partial correlation coefficient 
aS a measure of the importance of variables not yet in the equation. The 
basic procedure is to first select the X most correlated with Y (suppose it is 
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X,) and find the first-order, linear regression equation Y = £(X}). Next 
find the partial F values for each of the remaining variables to enter the 
equations. As each variable is entered into the regression the following 
values are examined. 


1. R42, the multiple correlation coefficient. 


2. The partial F-test26/ value for the variable most recently entered 
which shows whether the variable has taken up a significant amount of 
variation over that removed by variables previously in the regression. 


As soon as the partial F value related to the most recently entered variable 
becomes nonsignificant the process is terminated. 


The forward selection procedure is more economical of computer facilities 
than the methods previously discussed and it avoids working with more X's than 
are necessary while improving the equation at every stage. One of its 
disadvantages is that it makes no effort to explore the effect that the 
introduction of a new variable may have on the role played by a variable which 
entered at an earlier stage. This deficiency is overcome by the stepwise 
procedure. 


The Stepwise Regression Procedure. 


In spite of its entirely different name, this procedure is an improved 
version of the forward selection procedure discussed in the previous section. 
The improvements involve the reexamination, at every stage, of the regression 
of the variables incorporated into the model in previous stages. A variable 
which may have been the best single variable to enter at an early stage may, 
at a later stage, be superfluous because of the relationships between it and 
other variables in the regression. To check on this the partial F criterion 
for each variable in the regression at any stage of calculation is evaluated 
and compared with a preselected percentage point of the appropriate F 
distribution. This provides an estimate of the contribution made by each 
variable irrespective of its actual point of entry into the model. Any 
variable which provides a nonsignificant contribution is removed. This 
process is continued until no more variables are admitted or rejected. 


We believe this to be the best of the variable selection procedures and 
recommend its use. However, stepwise regression can easily be abused. As 
with all the procedures discussed, sensible judgment is still required in the 
initial selection of variables and in the critical examination of the model 
via examination of residuals. It is easy to rely too heavily on the automatic 
selection made by the computer. All variables tested should be consistent 
with "a priori"™ constructs. It is methodologically improper to use the 
iterative regression techniques to sort variables selected at random and then 
to claim that a hypothesis test has been performed. 


In a theoretical sense, the "all regressions" procedure is best because it 
enables one to "look at everything." Much would depend on the rejection 
levels for the F-tests, however, and also on the statisticians’ feelings about 


26/ The t test is equivalent in this instance. 
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the increase desired in R2 when all regressions are viewed. Inconsistent 
choices might lead to entirely different equations being reached. 


Our own preference, for a practical regression method, is the stepwise 
procedure. Our second choice is the backward elimination. In cases of doubt, 
all regressions could be viewed if the effort is reasonable. The other 
methods given are, in our opinion, of lesser value but each has its use in 
special cases. 


Conditional Forecasting 


Our previous discussion has implicitly assumed that the explanatory 
variables are known without error. This can be an unrealistic assumption 
since some explanatory variables may have to be predicted before the single 
equation regression model can be used. One may intuitively expect that the 
stochastic nature of the predicted value of the X's will lead to forecasts of 
Y which are less reliable. We shall show that the 95 percent confidence 
intervals for the error of forecast are indeed increased in size when the X's 
themselves must be predicted. However, it is quite difficult to derive 
analytic results for the error of forecast in a general setting. We deal, 
therefore, with a special case which, though somewhat restrictive, should be 
instructive. 


Consider the following model 


SO BX, ae = sh -Cinucnepsasy lh 


where 


et Sal N(O, o7), Ub N(O, ou2), a, 


E [tky, oo) Xo) Wapeon et = Bol (aoe Gale ae) n= eenand 


and U. are uncorrelated, 


Brel 


a and 8 are the OLS estimates of ao ands. 


The model presumes that Xp4,] is projected in such a way that its error of 
forecast has 0 mean and constant variance. In addition, the error process 
associated with the forecast of Xm4, is assumed to be independent of the 
error process associated with each of the Y's in the model. Even though they 
are stochastic in nature, the X's are still presumed to be exogenous in the 
sense that they are uncorrelated with the equation error term. The 
restrictiveness of this model becomes clear when we consider the means by 
Which forecasted values of Xp41 might be obtained. One frequently used 
procedure is to extrapolate from the sample values of X. While numerous 
extrapolation routines are available, the likelihood that the X variable is 
autocorrelated in a time-series model suggests that the error of forecast 
associated with the extrapolation procedure is itself likely to be serially 
correlated. 
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The forecasted value of Y in the time period T+1 is defined by 


The error of forecast is then 
en ee cen ee part et CF aK epyyees: (16) 


It is easy to see that this error has 0 mean by computing 


E(AT 1) a E(a-o) a E (aX TES Pe kT ia) i, E(e7,,) 
= EL B(X,) + Ya,)] - Xa) re eka, = 7) 

since : and Up4,] are uncorrelated. 

The variance of the error of forecast is 

eee Nee WER e ra)?) t EPC BR - BX. 7)7] 
2 
HCE (ery) ] + 2E Le = a) (4 -BX741)]. (18) 

But since 

OXze1 ~ OXtay = BUR, - Xtay) + Xtag(@ ~ 8) (19) 


then, 


Fe Drip Bogie B 2 
E [ (BX, aa BX 144) ] i ELa(X+,, = X44) 7 Xray (8 8) ] 


E[2(x 


fo aikaaee YOKE iE o08)? 


[e* + Var(8)o7 - Mead Var(8) (20) 


, in arriving at Equation (20), we took advantage of the fact that Un, = 


Amy] — Xp41 ‘that upy) and 8 are uncorrelated, and, finally, that 


g? = g2+ Var (8). 
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Next, we can use Equation (19) to simplify the last term in Equation (18), 
1.e., 


E[(a - a) (8X Bogq)] = EL(G = a) By, - Xray) 


fT 


+ Xoyq EL(G - 0) (8 - 8)] 


“aA “aw 


Sa Cov(a, 8). 


Now, by combining terms, we find that 


nw 


Die zi 2 és 2 fi Z 
o. = Var(a) + [8° + Var g)] Cite No Var(g) + 2X44 Cov(a, 8) + 07 


= Var (a) + Var(a) [x2 Hie ee Ed ava ee DP 
T+ * %] + 2k, Cov(a, 8) +o" + Bo. (22) 


Putting this in terms of the least-squares estimators the formula for the 
variance of the error of forecast becomes 


(Ncw =p n i Aa Meiers 
op=o fl +¢+—Hl Et] + ao (23) 
: (X4 = X) 
(x RA 
oe (X known) = o°[1 + Cer tea lit 5 Ay: (24) 
=(X, = 0) 


A comparison of Equation (23) with Equation (24) makes it clear that the 
stochastic nature of the forecasted value of X increases the error of 
forecast. There are two additional nonnegative terms involved, both of which 
are minimized only when the forecast of Xp4}] is made with 0 error 

variance. This result makes it quite clear that the 95 percent confidence 
intervals associated with the forecasted value of Yq4}) will be larger than 
the intervals described in the unconditional forecasting derivation. 
Unfortunately, it is quite difficult to describe the confidence intervals for 
the conditional error of forecast. The difficulty arises because the forecast 
variable Yq,; is not normally distributed since it involves the sum 

of products of normally distributed variables. While confidence intervals 
cannot be derived analytically they can be approximated using computer 
simulation techniques. A rather crude estimate of the confidence interval 
might be obtained as follows: 


l. Calculate the 95 percent confidence intervals associated with the 
forecast that would be obtained were we to select X41 to be two 
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standard deviations higher or lower, i.e., the confidence intervals 
associated with 


2. The final interval prediction is taken to be the union of the two 
confidence intervals, i.e., it contains all the values of Yp4j 
common to both confidence intervals. This process is depicted by 


Y APPROXIMATED PREDICTION 


pubs aA SECOND PREDICTION 
INTERVAL 


FIRST 
PREDICTION 
INTERVAL 


95% CONFIDENCE BAND 
ASSUMING X, KNOWN 


The Prediction Interval for a Conditional Forecast 


The results of this discussion help to accentuate some of the difficulties 
involved in the forecasting process. Even if the regression model has a good 
fit with statistically significant parameters, unconditional forecasts may not 
be very accurate. In addition, the process of forecasting or extrapolating 
the explanatory variables into the prediction period introduces additional 
forecast error which diminishes the reliability of the regression model. In 
forecasting, the problem is made even more difficult because an accurate 
initial conditional forecast may lead to a shift in policy and thus to an 
inaccurate forecast! 


Summary of Methodology 


A cross section data set across cities was selected because (1) demand 
predictions were required by utility service area and city data more closely 
approximated service areas than did state data, (2) deflators for cost-of- 
living are only available for cities and might be quite inaccurate if applied 
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to state data, (3) although time series data were collected the statistical 
difficulties inherent in time series analysis and the lack of consistency 
among previous time series studies justified an initial cross section 
analysis, (4) a large degree of variability in price and household consumption 
exists in the cross section data which allows forecasts without extrapolating 
outside the range of the data, and (5) cross section analysis incorporates the 
long run effects of price differences on both the rate of utilization of 
appliance stocks and the effects on the mix of electric versus alternative 
fuel appliance stocks without requiring measurement of stocks. Also, data on 
appliance stocks were suspected to contain errors of measurement. 


This study thus avoids the inappropriate application of city price 
deflators to state data, avoids the unresolved statistical problems of time 
series analysis, uses a sample which contains a highly desired wide range of 
prices and household consumption rates and avoids the necessity of using weak 
appliance stock data. The analysis was made possible by the acquisition of 
confidential, unpublished American Gas Association data by city. 


Other important unique attributes of the study include the use of a 
quadratic form for heating degree days, use of a Single variable to account 
for both heating and cooling requirements and the use of dummy variables to 
account for measurement error when missing city price indices were estimated 
by surrounding cities indices. 


A single equation log-linear model was formulated. While it is recognized 
that a supply-demand, two-equation model is more appropriate, the poor results 
of previous two-equation studies and the priority to forecast, rather than 
estimate, structural relationships led to the decision to use a Single 
equation model. Bias in the estimated coefficients is minimized if the error 
variance of the regression approaches zero. Unusually high coefficients of 
determination were attained for the demand estimates, leading to the 
conclusion that most of the remaining error variance was truly random. 

Neither the lack of a simultaneous model specification nor the exclusion of 
explanatory variables could cause bias in the estimated model under these 
circumstances. 


The success of the single-equation demand model is probably due to the 
more accurate deflation of monetary variables than exists in other studies. 
The use of "real" or constant-purchasing power monetary values across cities 
(or states, or time) is premised upon economic theory in the absence of either 
speculative purchasing or "money illusion." 


The use of the log-linear functional form circumvented the unresolved 
debate concerning whether marginal or average prices are appropriate in 
electricity demand and provided direct estimates of elasticity. The fit of 
the log-linear models was very good. However, the usual nontransference of 
statistical properties under nonlinear transformations applies to our demand 
estimates. 


CHAPTER IV. DATA DESCRIPTION 


The First Sample 


The estimating procedure was undertaken in two stages. The first stage 
involved the examination of a several alternative measures of the hypothesized 
causal variables to predict household electricity consumption. A list of the 
variables tested in the first stage is shown in table 6. The data set 
described by this list of variables is shown in the appendix. Observations 
over 102 cities were collected for the first stage investigation. The cities 
included in the first stage sample are listed on table 7. A number of 
alternative models were investigated using the 102 city data set. The results 
of the preliminary investigation led to two major conclusions. First, the 
specification for the demand equations presented later in this chapter was 
finalized and, second, a first approximation of a recursive model which could 
be estimated with recently acquired American Gas Association consumer price 
data was achieved. The latter research, which relied upon incomplete AGA 
data, supports further statistical study of a more complex system to explain 
and predict household energy consumption. A preliminary basis for this 
further research is presented in Chapter VI. 


Once a model appropriate for the prediction of household electricity 
consumption (given numerous trade-offs and compromises as discussed earlier) 
had been selected and tested using the 102 city sample, the second stage of 
the prediction effort could commence. The second stage consisted of the 
reestimation of the demand curve for household electricity consumption using a 
newly developed surrogate price index to adjust monetary variables in an 
expanded 138 city sample. The 138 city sample was more appropriate to the 
utility service areas for which household electricity consumption forecasts 
are required. The use of the surrogate deflator price index allowed the 
expansion of the deflated electricity demand equation estimates from a sample 
of 57 to the total 138 cities in the second stage sample. The use of deflated 
cross section data greatly improved the accuracy of the estimated demand 
equation (as expected from results obtained with the earlier experimentation 
with the 102 city sample). 


The Second Sample 


The sample data used to estimate the second stage household electricity 
demand function are shown in table 8. A short description of the variables 
follows. 

XPD -——deflated average household electricity price in 1975 

XINCD--deflated per family income in 1975 


AGAPD--deflated average household natural gas price in 1975 (data are 
confidential) 


HDD75--heating degree days in 1975 
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Table 6--List of variables for household electricity demand analysis 


Variable Explanation 

AGAP Gas price from AGA (household) 1975: $ per million Btu 

AGAPP Propane price from AGA (household) 1975: $ per million Btu 

AGAEP Electricity price from AGA (household) 1975: $ per million Btu 

AGAO2P #2 Fuel oil price from AGA (household) 1975: $ per million Btu 

AGAO1P #1 Fuel oil price from AGA (household) 1975: $ per million Btu 

AGACP Coal price from AGA (household) 1975: $ per million Btu 

POPSM Population/square mile (possible index of apartments with commercial 
rates) 

MEDAGE Median age: years 

MEDINC Median income: dollars 

HUM75 Average humidity at 1 AM for 1975: percent 

HDD75 Total heating degree days for 1975: days 

CDD75 Total cooling degree days for 1975: days 

AVHDD Average heating degree days for 1941-1970: days 

SSPC Percent of possible sunshine in 1975: percent 

AVDB Temperature exceeded 5% of the time (dry bulb thermometer): F 

AVWB Temperature exceeded 5% of the time (wet bulb thermometer): F 

ELEV Elevation: feet 

GP60 Gas price from AGA (household) 1960: cents per therm 

02P60 #2 fuel oil price from AGA (household) 1960: cents per gallon 

01P60 #1 fuel oil price from AGA (household) 1960: cents per gallon 

BCP60 Bituminous coal price from AGA 1960: $ per ton 

ACP60 Anthracite coal price from AGA 1960: $ per ton 

CKP60 Coke price from AGA 1960: $ per ton 

BP60 Butane price from AGA 1960: ¢ per gallon 

PP60 Propane price from AGA 1960: ¢ per gallon 

KP60 Kerosene price from AGA 1960: ¢ per gallon 

EP60 Electricity price from AGA 1960: mills per kWh 

NTOT60 Total occupied housing units 1960 

NGHF60 Number of houses with gas heat 1960 

NOHF'60 Number of houses with oil or kerosene heat 1960 

NCHF60 Number of houses with coal or coke heat 1960 

NEHF60 Number of houses with electric heat 1960 

NPHF60 Number of houses with propane heat 1960 

NOTHF60 Number of houses with wood or other 1960 

NONEHF60 Number of houses with no heat 1960 

NGWHF60 Number of houses with gas water heating fuel 1960 

NEWHF60 Number of houses with electricity water heating fuel 1960 

NCWHF60 Number of houses with coal or coke water heating fuel 1960 

NPWHF60 Number of houses with propane water heating fuel 1960 

NOWHF60 Number of houses with oil or kerosene water heating fuel 1960 

NOTWHF60 Number of houses with wood or other water heating fuel 1960 

NONEWHF60 Number of houses with no water heating fuel 1960 

NGCF60 Number of houses with gas cooking fuel 1960 

NECP60 Number of houses with electricity cooking fuel 1960 

NPCF60 Number of houses with propane heating fuel 1960 

NOCH60 Number of houses with oil or kerosene cooking fuel 1960 

NCCF60 Number of houses with coal or coke cooking fuel 1960 

NOTCF60 Number of houses with wood or other cooking fuel 1960 

NONECF60 Number of houses with no cooking fuel 1960 


Continued 


Variable Explanation 

NWRI60 Number of houses with wringer washer 1960 

NAUT60 Number of houses with automatic washer 1960 

NWD60 Number of houses with washer-dryer combination 1960 
NONEW60 Number of houses with no washer 1960 

NGD60 Number of houses with gas heated clothes dryer 1960 
NED60 Number of houses with electric heated clothes dryer 1960 
NONED60 Number of houses with no Clothes dryer 1960 

NFF60 Number of houses with one or more food freezers 1960 
NONEFF60 Number of houses with no food freezers 1960 

N1AC60 Number of houses with one air-conditioner 1960 

N2AC60 Number of houses with 2 or more air-conditioners 1960 
NCAC60 Number of houses with central air-conditioning 1960 
NONEAC60 Number of houses with no air-conditioning 1960 

N1TV60 Number of houses with ] television set 1960 

N2TV60 Number of houses with 2 or more television sets 1960 
NONETV60 Number of houses with no television sets 1960 

TRIND75 Revenue large industrial customers 1975 (FPC electricity) 
MUNI75 Dummy indicating municipal retailer of electricity 
TEB50075 Typical electric bill for 500 kWh 1975 

TEB75075 Typical electric bill for 750 kWh 1975 

TEB10007 Typical electric bill for 1000 kWh 1975 

NDTOT7Q #$Total number of dwellings 1970 

NAC170 Number of dwellings with air-conditioning in 1] room 1970 
NAC270 Number of dwellings with air-conditioning in 2 or more rooms 1970 
NACC70 Number of dwellings with central air-conditioning 

MEDR70 Median rent 

NRESI75 Number of residential customers 1975 FPC 

NCOM75 Number of commercial customers 1975 FPC 

NIND75 Number of industrial customers 1975 FPC 

KWHRES75 Residential electric consumption 1975 FPC: thousands of kWh 
KWHCOM75 Commercial electric consumption 1975 FPC: thousands of kWh 
KWHIND75 Industrial electric consumption 1975 FPC: thousands of kWh 
TRRESI75 Revenue residential customers 1975 (FPC electricity) 
TRCOM75 Revenue commercial customers 1975 (FPC electricity) 
NDWELTOT Total dwellings 1970 

B6970 Dwellings constructed 1969-—March 1970 

B6568 Dwellings constructed 1965-1968 

B6064 Dwellings constructed 1960-1964 

B5059 Dwellings constructed 1950-1959 

B4049 Dwellings constructed 1940-1949 

BEARLY Dwellings constructed 1930 or Earlier 

NTOT70 Total occupied housing units 1970 

NGHF70 Number of houses with gas heat 1970 

NOHF70 Number of houses with oil or kerosene heat 1970 

NCHF70 Number of houses with coal or coke 1970 

NWHF'’70 Number of houses with wood heat 1970 

NEHF'70 Number of houses with electric heat 1970 

NPHF70 Number of houses with propane heat 1970 
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Table 6—-List of variables for household electricity demand 


analysis (continued) 


Continued 


60 


Table 6--List of variables for household electricity demand 
analysis (continued) 


Variable 


NOTHF70 
NONEHF70 
NGWHF 70 
NOWHF'70 
NCWHE'70 
NWWHF'’70 
NEWHF70 
NPWHF70 
NOTWHF70 
NONEWHF70 


NGCF70 
NECF70 
NPCF’70 
NOCF70 
NCCF70 
NWCF70 


NOTCF70 
NONECF'’70 


NWR70 
NAUT70 


NONEW70 


NGD70 
NED70 


NONED70 


NDW70 


NONEDW70 


NFF70 


NONEFF70 


N1TV70 
N2TV70 


NONETV70 


Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 


S21 through S68 
48 dummy variables indicating state of Residence 


houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 


with 
with 


Explanation 


other heat 1970 
no heat 1970 


with gas water heating fuel 1970 


with 
with 


oil or kerosene water heating fuel 1970 
coal or coke water heating fuel 1970 


with wood water heating fuel 1970 


with 
with 
with 
with 
with 
with 


electricity water heating fuel 1970 
propane water heating fuel 1970 
other water heating fuel 1970 

no water heating fuel 1970 

gas cooking fuel 1970 

electricity cooking fuel 1970 


with propane cooking fuel 1970 


with 
with 
with 
with 
with 


oil or kerosene cooking fuel 1970 
coal or coke cooking fuel 1970 
wood cooking fuel 1970 

other cooking fuel 1970 

no cooking fuel 1970 


with wringer washer 1970 


with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 


automatic washer 1970 

no washer 1970 

gas heated clothes dryer 1970 
electricity heated clothes dryer 1970 
no clothes dryer 1970 
dishwasher 1970 

no dishwasher 1970 

food freezer 1970 

no food freezer 1970 

1 television set 1970 

2 or more television sets 1970 
no television sets 1970 


State Dummies 
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Table 7--Cities included in first stage sample 


WM~AAYHANO FWHM YE 
e e 


AT. 


49. 
oan 


Albuquerque, New Mexico 
Anaheim, California 
Appleton, Wisconsin 
Atlanta, Georgia 

Atlantic City, New Jersey 
Baltimore, Maryland 
Billings, Montana 
Binghamton, New York 
Birmingham, Alabama 
Boise, Idaho 

Boston, Massachusetts 
Bridgeport, Connecticut 
Buffalo, New York 
Charleston, West Virginia 
Charlotte, North Carolina 
Chicago, Illinois 
Cincinnati, Ohio 
Cleveland, Ohio 

Colorado Springs, Colorado 
Columbia, Missouri 
Columbia, South Carolina 
Columbus, Ohio 

Dallas, Texas 

Dayton, Ohio 

Decatur, Illinois 

Denver, Colorado 

Des Moines, Iowa 

Duluth, Minnesota 

El Paso, Texas 

Erie, Pennsylvania 
Evansville, Indiana 

Fall River, Massachusetts 
Fargo, North Dakota 
Fayetteville, North Carolina 
Flint, Michigan 

Fort Lauderdale, Florida 
Fort Wayne, Indiana 

Fort Worth, Texas 
Galveston, Texas 

Great Falls, Montana 
Green Bay, Wisconsin 
Greenville, South Carolina 
Harrisburg, Pennsylvania 
Hartford, Connecticut 
Houston, Texas 
Indianapolis, Indiana 
Jackson, Mississippi 
Kansas City, Missouri 
Knoxville, Tennessee 
Lafayette, Indiana 
Lafayette, Louisiana 


D2. 
53. 
54. 
55. 
56. 
Oi. 
58. 
3)5 
60. 
6l. 
62. 
63. 
64. 
65. 
66. 
67. 
68. 
69. 
70. 
ie 
a2. 
2c 
74. 
15. 
76. 
77. 
78. 
Le 
80. 
Sl. 
82. 
83. 
84. 
85. 
86. 
87. 
88. 
89. 
90. 
91. 
O27. 
93. 
94. 
ules 
96. 
as 
98. 
99... 
100. 
161. 
102. 


Lake Charles, Louisiana 
Lexington, Kentucky 
Lincoln, Nebraska 

Los Angeles, California 
Louisville, Kentucky 
Phoenix, Arizona 
Manchester, New Hampshire 
Memphis, Tennessee 
Milwaukee, Wisconsin 
Minneapolis, Minnesota 
Muncie, Indiana 

Nashua, New Hampshire 
Nashville, Tennessee 

New Orleans, Louisiana 
New York City, New York 
Newark, New Jersey 
Norfolk, Virginia 
Oklahoma City, Oklahoma 
Omaha, Nebraska 

Orlando, Florida 
Owensboro, Kentucky 
Peoria, Illinois 
Philadelphia, Pennsylvania 
Phoenix, Arizona 
Pittsburgh, Pennsylvania 
Pittsfield, Massachusetts 
Portland, Maine 

Portland, Oregon 
Providence, Rhode Island 
Racine, Wisconsin 
Raleigh, North Carolina 
Richmond, Virginia 
Roanoke, Virginia 
Rochester, Minnesota 
Rochester, New York 

St. Louis, Missouri 

Salt Lake City, Utah 

San Antonio, Texas 

San Diego, California 

San Francisco, California 
San Jose, California 
Scranton, Pennsylvania 
Seattle, Washington 

Sioux City, Iowa 

Sioux Falls, South Dakota 
Springfield, Massachusetts 
Topeka, Kansas 

Tucson, Arizona 
Washington, District of Columbia 
Waterbury, Connecticut 
Wilmington, Delaware 
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XPINA--a unit dummy variable which indicates those cities for which the 
Surrogate price deflator must be used 


XPI --the price deflator 
SPI -—the surrogate price deflator 
TRDN-~actual deflated spending per family for household electricity 


Although not used directly in the statistical analysis the rather large 
expenditure variation by city is revealing. Average real spending (by city), 
per family for household electricity in 1975, varied from $151 to $545. It is 
the task of the econometric model to explain this variation. 


Table 9 shows the means and standard deviations of the variables included 
in the stage two regression equations. A rather large variability can be 
noted in both the deflated average electricity price across cities and in the 
electricity consumption per family. In each case, the standard deviation is 
more than one quarter of the mean value. The large variation in the demand 
variables (deflated electricity price and consumption per family) implies that 
the cross section data reflects the desired long run comparison (in which the 
sample data include a wide range of price-quantity sets). The large variation 
in the sample data is desired because future forecasts will not require 
extrapolating as far beyond the range of the sample data as would be the case 
if the less variable time series sample data were used. It is also easily 
shown also that large variation in the dependent variable results in less 
sensitivity of the estimated equation to slight errors in the independent 


Table 9--Means and standard deviations of sample data used in 
regression estimation* 


Standard 

Variable Mean Deviation 
Deflated Electricity Price 
$/Thousand kWh (Avg. 1975) 38.6370 10.5864 
Electricity Consumption Per 
Family, Thousand kWh/year (1975) 8.6055 2.5328 
Deflated Income Per Household | 
$/Year (1975) 11,853.88 1 seen 
Deflated Natural Gas Price 
$/Million BTU (Avg. 1975) 1.9413 0.6177 
Heating Degree Days (1975) 2,915.93 1,041.22 
Heating Degree Days Squared 9,578,923.8 6,193,885.3 


*Sample data include 57 observations which are deflated using cross-city cost 
of living indices and 8] observations deflated by surrogate price indices 
calculated from nearby cities. State dummy variables are not shown. 
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variables. The increased robustness of the estimated model is a further 

desired outcome. Natural gas prices also show extreme variation among cities (the 
standard deviation is almost one-third of the mean). High variation in the sample 
natural gas price data is particularly desirable since DOE forecasts of future 
consumer energy prices show that natural gas prices are expected to grow very 
rapidly. 


CHAPTER V. FINAL DEMAND EQUATION AND PREDICTION METHODOLOGY 


Six equations are presented in the second stage of the empirical 
investigation. All of the final equations may be described as log-linear in 
the monetary variables with quadratic terms included to show the effect of 
heating degree days on household electricity consumption. Since part of the 
data were deflated using actual city price deflators and part were deflated 
using surrogate deflators, separate and pooled regressions are presented to 
allow some adjustment for any bias resulting from the inclusion of surrogate 
deflators. The general form of the estimated demand equations is shown in 
both the log--linear form (the form estimated by computer using the linear 
regression technique) and the equivalent product equation form. While the 
log-linear form is necessary to apply the linear regression technique to 
estimate the equation coefficients, the product form is commonly referred to 
in discussions. As is well known, the exponents Kj through Kg represent 
the elasticities of the variables to which they are attached. For example, 
K) indicates the price elasticity of demand for the variable electricity 
price (deflated by the price index). Similarly, the coefficient Kj 
indicates the price elasticity estimated in cities where the surrogate price 
deflator was required. Coefficients K3 and K4 show the respective income 
elasticities of demand and coefficients Ks and Kg show the respective 
cross elasticities of electricity sales to natural gas price. Since the 
heating degree days variable is of quadratic, rather than log-linear form, no 
direct elasticity interpretation of coefficients K7 and Kg is possible. 

The variables 9 and 1- P appear in the equations as a means of statistically 
identifying the elasticities estimated for cities with real price deflators 
and those with surrogate price deflators. The variable @ takes the value of 
unity when the real deflator log-linear version of the demand equation reveals 
that the @ and 1-9 variables are premultiplied by each of the independent 
monetary variables. Thus, when only a surrogate deflator is available, the 
first term observation (following K)) is zero. If the real deflator is 
available, the second term observation (following Kj) is zero. In this 
manner, separate estimates of price, income, and cross elasticity are provided 
for cities with each type of deflator. The use of a dummy variable indicating 
a real or surrogate deflator is not required for the HDD variable since it is 
a nonmonetary variable not requiring adjustment to a common measure. The 
equation in log-linear form is 


Whi pen years 
In SRovean a. Kral kK; (0) In Pe 


customer PT SPT 
+ Ky (0) Tn eee 
Pa SPI 
Pn 
+ Ks (D) In 509 + ke (1-9) In eng 
48 
+ Ky HOO Ke COD a) papi aMiediuns 
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The equation in product form is 


customer OPI SPI PI SPI 
a : te 
Tes K HDD + K HDD “ . 
Bee nga ie” a Ae 
PI SPI 


The Demand Models 

Six demand equations are presented in the following tables. Each of the 
models fits the data well (the coefficient of determination varies from .842 
to 0.978). A coefficient of determination in the 90 percent range is 
extremely good for cross section data which excludes most spurious correlation 
often caused by contemporaneous trends in economic time series data. In fact, 
a higher correlation may be suspect since measurement error is known to exist 
and part of the variation in consumer purchases may be truly random and not 
susceptible to prediction. 


Model I (table 10) is the subsample of 57 cities for which the actual 
price deflator data are available. Regressions on the data which are deflated 
by the actual price deflator tend to fit more precisely than do the 
regressions where the surrogate deflator is required. (When no deflation at 
all is used the fit of the estimated equations falls drastically.) Model II 
(table 11) is similar to Model I except that it is applied to the sample data 
of 81 cities where the surrogate deflator is required. The R2 falls from 
-95 to .842, partly, because of measurement error when the deflation is 
accomplished with the surrogate deflator. Models I and II incorporate dumny 
variables which shift the demand equation up or down to adjust for variation 
among states in the rate of household electricity consumption. The state 
dummy variables are included to adjust for geographic market or climatic 
factors for which more specific data are unavailable. The inclusion of dummy 
variables does not add to the forecasting difficulty usually encountered when 
the set of independent variables is enlarged. No prediction error exists for 
this type of variable. 


The earlier analysis of specification error suggests some of the serious 
difficulties involved with the model-building process. If we are unsure of 
which explanatory variables ought to appear in a model, we face several 
tradeoffs. The cost of excluding a variable which should appear in the model 
is bias and inconsistency. If the number of observations available is rather 
large, it seems reasonable to opt for the risk of adding irrelevant variables 
because the loss of degrees of freedom associated with the addition of several 
independent variables to the model is unlikely to be serious. If all 
potential model variables are known, the choice of model form must be made in 
terms of the bias-efficiency trade-off with the result dependent upon our 
objectives. If accurate forecasting is the goal, minimizing mean square error 
appears to be one reasonable objective since it accounts for both bias and 
efficiency (Mean Square Error = Variance + (bias)2). Thus, in Principle, we 
might estimate each of several alternative models and compare the mean square 
errors associated with each model. 
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Table 10--Demand Model I 
Oy peetgentaal Electricity Demand: Subsample with Real Price Deflator 
Dependent Variable: Log Transform of kWh per customer per year 
Sample Size: 57 SMSA's for the year 1975 


Functional Form: Log linear with quadratic term for weather 


Deflation of monetary variables: Prices and income deflated with a 
price index 


Coefficient of Determination: R2 = 0.950 (adjusted) 


Independent Variables Estimated Coefficient t-Value 
Electricity Price deflated 

by a price index (log transform) -0.843 - 17.05 
Income deflated by a price 

index (log transform) 0.634 28.89 
Natural Gas Price deflated by 

a price index (log transform) —0.144 - 2.18 
Heating Degree Days -0.000338 - 6.88 
Heating Degree Days Squared 0.0000000294 SVay7Ak 


State Dummy Variables: 


New York S28 —-0.141 ‘ie 2283 
Pennsylvania S29 0.074 1.91 
Tennessee S54 0.202 2.38 
Texas S58 -0.226 - 3.64 
Colorado S60 —0.401 - 5.80 
California S66 -0.835 -12.97 
Washington S68 -0.575 -— 7.49 


Constant Term Excluded Excluded 


yal 
Table 11--Demand Model II 
i} eaeedaanevar Electricity Demand: Subsample with Surrogate Price Deflator 
Dependent Variable: Log transform of kWh per customer per year 
Sample Size: 81 SMSA's for the year 1975 


Functional Form: Log linear with quadratic term for weather 


Deflation of monetary variables: Prices and income deflated with a 
Surrogate price index 


Coefficient of Determination: R2 = 0.842 (adjusted) 


Independent Variables Estimated Coefficient t-Value 


Electricity Price deflated by 
surrogate price index (log 
transform) —0.689 —11.69 


Income deflated by surrogate 
price index (log transform) 0.564 25 


Natural Gas Price deflated by 
surrogate price index (log 


transform) 0.175 4.76 
Heating Degree Days —0.000362 - 7.00 
Heating Degree Days Squared 0.0000000338 4.10 


State Dummy Variables: 


Massachusetts S22 —-0.342 - 3.42 
New Hampshire S24 0.130 1.80 
New York S28 —0.188 - 2.76 
Pennsylvania S29 0.149 2.47 
Iowa S35 0.204 3.97 
Minnesota Son 0.086 e235 
Deleware. S43 0.285 2.05 
Kentucky S52 =e - 5.18 
Tennessee S54 0.151 20 
Texas S58 -0.417 - 4.35 
New Mexico S63 —-0.269 - 2.81 


Constant Term Excluded Excluded 
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In terms of classical statistics, it is not difficult to test whether 
irrelevant variables are present. Since the coefficients of irrelevant 
variables have expected values of 0, we need simply apply standard tests if we 
wish to evaluate the relevancy of individual variables. Such testing is 
feasible if we know the set of variables which might conceivably appear in the 
model but it fails completely when we are ignorant to the underlying 
causation. 


Model III (table 12) combines the sample data included in Models I and II 
for a sample size of 138 cities. Dummy variables are used in the manner 
described earlier to allow pooling of the data while providing separate 
elasticity estimates for monetary variables deflated by real and surrogate 
price indices. The R2 of .903 falls between that of the first two models as 
might be expected. The elasticity estimates are very similar to those found 
in the two separate regressions in Model I and Model II. It is noteworthy 
that an anomaly appears in the demand functions for the cross elasticity of 
electricity demand to the price of natural gas deflated by the actual price 
index. In cities where a real price deflator is available the data indicate 
that high prices for natural gas are associated with a lower per family 
consumption of electricity. The sign of this variable is "wrong" if natural 
gaS is to be considered as a substitute, rather than a complement, to 
electricity. This theoretically unexpected finding has also been reported in 
several previous studies.(20) This serious error in estimation is removed in 
Model VI (table 13). Model VI does not differentiate among "real" or 
"surrogate" price deflators. Also, a second weather variable, humidity, is 
included. Model VI also is adjusted for heteroscedasticity by dividing all 
variables by the square root of real income per household. The resulting 
demand equation has an adjusted R2 of 90 percent. All variables in Model VI 
have the correct sign for their estimated coefficients. 


Models IV and V (tables 14 and 15, respectively) are Similar to Models I 
and II except that the stepwise regression program was calibrated to allow 
state price elasticity dummy variables to enter after (usually) the entry of 
the state dummy shift variables. In terms of the log-linear form of the 
model, this allows both the intercept and the slope coefficients to vary by 
state. Inclusion of both effects tends to lead to high correlation among 
independent variables and possible estimation errors. Price elasticity dummy 
variables with high intercorrelations with the shift dummy variables were 
excluded from the regressions. Further experimentation to extend the 
elasticity dummy variables to natural gas prices is warranted given the 
unexpected sign estimated for the cross elasticity discussed earlier. 


Although Models IV and V incorporate several added independent variables 
which adjust for cross state variation in price elasticity of demand, the net 
addition to the explanation is quite small. Model I has an R2 of 0.95 which 
rises to 0.978 in Model IV. Model II has an R2 of .842 which rises to .869 
in Model V. The added price elasticity variables are nominally statistically 
Significant but the high intercorrelations reduce the credibility of continued 
additions of dummy variables to the models. The forecasts will be little 
affected by the inclusion of the elasticity adjusting dummy variables. 


Model VI, which pools the two data sets contained in Models I and II, has 
been selected as the prediction equation best suited to provide consumption 
forecasts for 136 utility service areas. Although the R2 for Model VI is 
Slightly less than for Model I or Model IV, these models exclude 81 
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Table 12--Demand Model III 


Residential Electricity Demand: Total sample (model combining subsample 
with surrogate deflators and subsample with real deflators) 


Dependent Variable: Log-transform of kWh per customer per year 
Sample Size: 138 SMSA's for the year 1975 
Functional Form: Log-linear with quadratic term for weather 


Deflation of monetary variables: Prices and income deflated by a price 
index or a surrogate price index 


Coefficient of Determination: R2 = 0.903 (adjusted) 


Independent Variables Estimated Coefficient t-Value 
Electricity Price Deflated by 
Price Index (log transformation) -0.856 -16.67 
Electricity Price Deflated by 
Surrogate Price Index (log trans.) -0.655 -13.15 
Income Deflated by Price Index 
(log transformation) 0.644 30.85 
Income Deflated by Surrogate Price 
Index (log transformation) 0552 26.52 
Natural Gas Price Deflated by 
Price Index (log transformation) -0.170 - 2.43 
Natural Gas Price Deflated by 
Surrogate Price Index (log trans.) Ol 7 5435 
Heating Degree Days -0.000373 -11.08 
Heating Degree Days Squared 0.0000000368 6.86 


State Dummy Variables: 


Massachusetts S22 —0.360 - 3.89 
New York S28 -0.162 - 3.47 
Pennsylvania S29 0.103 3.00 
Towa S35 0.192 4.20 
Deleware $43 0.269 2.92 
Kentucky $52 -0.212 ohne yoe | 
Tennessee S54 0.179 3224 
Texas S58 0.286 5.39 
Colorado S60 —0.425 - 5.49 
New Mexico S63 -0.271 - 3.04 
California S66 0.842 -12.30 
Washington S68 -0.577 - 6.99 


Constant Term Excluded Excluded 
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Table 13--Demand Model VI 


Residential Electricity Demand: Total sample (no adjustment for 
variation between coefficients estimated on variables deflated by "real" 
versus "Surrogate" price indices) 

Dependent Variable: Log-transform of kWh per customer per year 

Sample Size: 138 SMSA's for the year 1975 


Functional Form: Log-linear with quadratic term for weather 


Deflation of monetary variables: Prices and income deflated by a price 
index or a Surrogate price index 


Coefficient of Determination: R2 = .906 (adjusted) 


Independent Variables Estimated Coefficient t-Value 
Electricity Price deflated by real 
Or surrogate price indices -0.718 =19.05 


Income deflated by real or 
Surrogate price indices 0.345 1 Ie5 


Natural Gas Price deflated by 


real or surrogate price indices 0.102 2.811 
Heating Degree Days -0.000209 - 0.610 
Heating Degree Days Squared 0.0000000172 3.16 

Humidity 0.440 © 8.296 


State Dummy Variables: 


Massachusetts S22 —3.564 - 4,01 
Rhode Island S25 —-2.073 - 2.30 
New York S28 —2./64 - 6.14 
Illinois S30 -0.839 - 2.57 
Michigan $32 -1.029 - 2.01 
Wisconsin S34 -1.153 - 2.82 
Towa S35 1.487 3u39 
Deleware S43 2.489 Pl S\D) 
Kentucky S52 —2.136 = D6 2 
Tennessee S54 1.746 Bos 
Arizona S59 6.435 8.72 
Idaho S6l 1.766 2.3 
California S66 -6.519 = 19553 
Washington S68 -5.00 = Wo)! 


Constant Term Excluded Excluded 
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Table 14--Demand Model IV 


Residential Electricity Demand: Subsample with real price deflators 
(model with price elasticity dummy variables) 


Dependent Variable: Log-transform of kWh per customer per year 
Sample Size: 57 SMSA's for the year 1975 
Functional Form: Log-linear with quadratic term for weather 


Deflation of monetary variables: Prices and income deflated with a 
Price index 


Coefficient of Determination: R2 = 0.978 (adjusted) 


Independent Variables Estimated Coefficient t-Value 

Electricity Price deflated by 

a price index (log transformation) ~0.782 -19.67 
Income deflated by a price index 

(log transformation) 0.614 35575 
Natural Gas Price deflated by 

a price index (log transformation) -0.177 = 2.97 
Heating Degree Days -0.000471 -12.06 
Heating Degree Days Squared 0.0000000615 7.82 


State Dummy Variables: 


New York S28 —0.072 - 1.69 
Pennsylvania S29 0.168 5.67 
Minnesota S87) -—0.096 - 1.62 
Tennessee S54 0.367 5287 
Texas S58 -0.115 - 2.38 
Colorado S60 -0.368 - 6.32 
California S66 -0.761 -17.00 
Washington S68 -0.375 -~ 5.78 


Price Elasticity Dummy Variables: 


Michigan PS32 0.018 2.04 
Ohio PS33 0.043 5.58 
Mississippi PS38 0.031 Sy 8 
District of 

Columbia PS42 0.074 4.93 
Maryland PS46 0.031 2, VE 
South Carolina PS47 0.081 55355 


Constant Term Excluded Excluded 
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Table 15-~Demand Model V 


Residential Electricity Demand: Subsample with surrogate price deflator 
(model with price elasticity dummy variables) 


Dependent Variable: Log transform of kWh per customer per year 
Sample Size: 81 SMSA's for the year 1975 
Functional Form: Log linear with quadratic term for weather 


Deflation of monetary variables: Prices and income deflated with a 
Surrogate price index 


Coefficient of Determination: R2 = 0.869 (adjusted) 


Independent Variables Estimated Coefficient t-Value 


Electricity Price deflated by 
surrogate price index 
(log transformation) —0.687 12.31 


Income deflated by surrogate 
price index (log transformation) 0.563 22.84 


Natural Gas Price deflated by 
Surrogate price index (log 


transformation) 0.164 4.51] 
Heating Degree Days -0.000365 - 8.12 
Heating Degree Days Squared 0.0000000370 5.40 


State Dummy Variables: 


Massachusetts S22 —0.350 - 3.78 
New Hampshire S24 0.096 1.47 
New York S28 —0.212 - 3.40 
Pennsylvania S29 On 32 2.39 
Towa S35 0.174 S15 (5y/ 
Deleware S43 0.280 3.05 
Kentucky S52 =0.221) = ale” 
Tennessee S54 0.151 2.29 
Texas S58 -—0.413 - 4.71 
New Mexico S63 —-0.279 - 3.19 


Price Elasticity Dummy Variables: 


Rhode Island PS25 -0.059 - 2.50 
Montana PS62 -0.049 - 2.83 


Constant Term Excluded Excluded 
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observations on cities where a "surrogate" deflator was required. Model VI 
has an R2 equal to that of Model III which included the total 138 

observation sample with adjustment for "real" or "Surrogate" deflators. Model 
VI is superior to Model III which estimated the wrong sign for natural gas 
prices deflated by the "real" price index. The use of weighted least squares 
to adjust for heteroscedasticity increased the efficiency of Model VI. A 
second weather variable, humidity, was also significant in Model VI. All 
signs of estimated coefficients in Model VI were of the theoretically correct 
Sign and reasonable in magnitude. Given the "noise" present in the data 
deflated by "surrogate" price indices, an R2 of 90 percent is extremely good 
for a cross section regression estimate. The required historical data and 
projected data inputs and the resulting econometric model and consumption 
forecasts are summarized in Figure 5. 


Comparison of Elasticity Estimates 


The cross section regression results presented in Tables 2 and 3 are the 
most appropriate models to compare with Model VI. The seven cross section 
models have price elasticities varying from -0.85 to -1.33. The simple 
average of the ten estimates (some studies provide a range of estimates) is 
-1.03. Our price elasticity estimate of -0.72 is slightly lower than any of 
the previous estimates. Cross elasticity estimates to natural gas price vary 
from -0.31 to +0.31. Excluding the estimate with the wrong sign, the average 
of the cross elasticity estimates is 0.207. Previous cross elasticity 
estimates are also higher than our estimate of 0.10. Previous income 
elasticity estimates vary from -0.46 to 1.13, with an average (excluding the 
estimate with the wrong sign) of 0.59. Our income elasticity estimate of 0.35 
is below the average of other pure cross section models. 


Although our elasticity estimates tend to be smaller in absolute value 
than previous cross section estimates, few of the cross section estimates 
differ markedly (excluding those with wrong signs). Several factors could 
account for our smaller elasticity estimates. Previous studies used data for 
1969 or 1970; or in one case 1959, 1965 and 1970; whereas our data are for 
1975. Some models measured electricity price differently although no 
systematic effect on the size of the elasticity estimates could be 
distinguished. Some previous models failed to deflate the monetary variables 
and the quality of most deflators is poor. Our natural gas price data are 
from a source not previously utilized (confidential AGA data). Our model 
contained several more independent variables than previous studies. Our Model 
VI may also be contrasted with our Models I-V. Generally our Models I-V have 
higher elasticity estimates more closely resembling the other cross section 
studies discussed above. The inclusion of an additional weather variable and 
the use of weighted least squares distinguish our Model VI from both Models 
I-V and earlier studies. Specification bias is a possible cause of the higher 
elasticity estimates found in models containing fewer significant explanatory 
variables. In any event, the range of cross section estimators is relatively 
small. 


The range of elasticity results when pure time series or pooled cross 
section and time series data sets were previously compared and shown 
previously in table 4. 
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Independent Variable Extrapolation Projection of Electricity Price 


The 1985 Series C electricity price (in 1978 dollars) is projected by the 
Department of Energy by region. The projections assume no natural gas 
regulation and an oil import price of only $15.00. It was necessary to divide 
all 1985 prices by 1.19 (Gross National Product implicit deflator) to adjust 
them to 1975 dollars (reversing the adjustment made by DOE). 


The DOE electricity price projections were made for 10 regions: New 
England, New York-New Jersey, Mid-Atlantic, South Atlantic, Midwest, 
Southwest, Central, North Central, West and Northwest. Table 16 shows the 
states contained in each of the above regions and Table 17 shows the sample 
cities used to estimate the electricity demand equation and the cities for 
which projections are required. 


Historical prices for the same regions were provided by the DOE for 1975. 
Key assumptions made by DOE are shown in Tables 18 and 19. These prices are 
in current or 1975 dollars. The percentage change (plus or minus) in 
electricity prices in 1975 dollars was calculated using the historical and 
projected prices provided by the DOE. The percentage change in electricity 
price for a region was applied to the sample price data used to estimate the 
electricity demand equation via a computer program designed to predict 
electricity sales by city. In this manner, DOE electricity price change 
projections were applied to our electricity price data in order to utilize our 
estimated demand equation to predict future electricity sales. The projected 
1985/1975 price ratios for residential electricity by DOE region are shown in 
Table 20. 


Projection of Natural Gas Price 


_ DOE historical and projected natural gas prices were used in the same 
Manner aS described above for the electricity price projections. The 
projected 1985/1975 price ratios for residential natural gas are shown in 
Table 20. 


The technique used by the Department of Energy to project household 
electricity and natural gas prices is briefly described in Energy Supply and 
Demand_in the Midterm: 1985, 1990, and 1995, prepared by the Office of 
Integrative Analysis, Midterm Analysis Division, April 1979. A mathematical 
computer model is used to make the midterm price projections. The computer 
model, Midterm Energy Forecasting System (MEFS), simulates the interactions of 
energy suppliers and consumers in the market place. Production and 
consumption responses to energy prices and explicit representation of the 
conversion of energy materials from one form to another for final consumption 
are contained in the model. The MEFS models price responses of producers, 
converters and consumers and projects price and ayant ety for each fuel at 
important points from extraction to final consumption.2// The medium 
scenario, "C", as described in Table 18, was selected for our use in 
projecting residential electricity and natural gas prices by DOE region. 
Important assumptions underlying all of the DOE scenarios are described on 


27/ Data used here are from special computer printouts provided by Terry 
Morlan, Director of the Demand Analysis Division, Office of Energy Use 
Analysis. 
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Table 16--List of States by DOE region* 


New England Southwest 
Connecticut Arkansas** 
Maine Louisiana 
Massachusetts New Mexico 
New Hampshire Oklahoma 
Rhode Island Texas 
Vermont** 

Central 
New York-New Jersey 

Iowa 
New Jersey Kansas 
New York Nebraska 

Missouri 


Mid-Atlantic 
North Central 


Delaware 
Maryland Colorado 
Pennsylvania Montana 
Virginia North Dakota 
Washington, D.C. South Dakota 
West Virginia Utah 
Wyoming** 
South Atlantic 
West 
Alabama 
Florida California 
Georgia Nevada** 
Kentucky Arizona 
Mississippi 
North Carolina Northwest 
South Carolina 
Tennessee Washington 
Oregon 
Midwest Idaho 
Illinois 
Indiana 
Michigan 
Minnesota 
Ohio 
Wisconsin 


*Also includes Washington, D.C. Alaska and Hawaii are not shown. 


**State not included in sample. 
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Table 17--List of sample cities by DOE region* 


New England 


Bridgeport 
Waterbury 
Hartford 
Portland 
Boston 
Fall River* 
Pittsfield 
Springfield 
Nashua* 
Manchester* 
Providence* 


New York-New Jersey 


Atlantic City* 
Newark* 

New York* 
Rochester* 
Binghamton* 
Buffalo* 


Mid-Atlantic 


Wilmington* 
Baltimore* 
Pittsburgh* 
Philadelphia* 
Scranton* 
Harrisburg* 
Erie* 
Richmond* 
Roanoke* 
Norfolk* 
Arlington 
Washington* 
Charleston* 


South Atlantic 
Birmingham* 


Fort Lauderdale 
Atlanta* 


Lexington* 
Louisville* 
Owensboro* 
Jackson* 
Charlotte* 
Raleigh* 
Fayetteville 
Greenville* 
Columbia* 
Knoxville* 
Memphis* 
Nashville* 


Midwest 


Peoria* 
Decatur* 
Chicago* 
Muncie* 
Indianapolis* 
Lafayette* 
Fort Wayne* 
Evansville* 
Flint* 
Duluth* 
Minneapolis* 
Rochester* 
Cincinnati* 
Dayton* 
Cleveland* 
Columbus* 
Milwaukee* 
Appleton* 
Green Bay* 
Racine* 


Southwest 


Lafayette 
Lake Charles 
New Orleans 
Albuquerque* 
Oklahoma City* 
Fort Worth* 


Houston 
Galveston 
Dallas* 

San Antonio 
El Paso* 
Lubbock 


Central 


Des Moines* 
Sioux City* 
Topeka* 
Lincoln* 
Omaha* 
Kansas City* 
St. Louis* 
Columbia* 


North Central 


Billings* 
Colorado Springs* 
Denver* 

Great Falls* 
Fargo* 

Sioux Falls* 

Salt Lake City 


West 


Anaheim* 

San Francisco 
San Diego 
Los Angeles* 
San Jose 
Phoenix* 
Tucson* 


Northwest 


Seattle* 
Portland* 
Boise* 


*Cities marked with an asterisk make up the set for which 1985 residential 


electricity sales are predicted. 


Table 18-—-Key assumptions by scenario 
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A B Cc D E 
I. Demand Curve 

Assumptions High High Medium Low Low 

Macroeconomic High- High- Trend- Low- Low- 

Forecast trend trend long trend trend 
II. Supply Curve 

Assumptions High Low Medium High Low 

Undiscovered Oil 

& Gas Resources High Low Medium High LOw 

New Technology 

Potential Penetration Medium High Medium Low Medium 

Nuclear Power Maximum 

Additions to Capacity High Low Medium High Low 
III. Cost Assumptions Low High Medium Low High 

Enhanced Oil 

Recovery Costs LOW High Medium Low High 

PACTEX Pipeline Tariff Low High Medium Low High 

Railroad Tariffs Low High Medium Low High 

Capital Cost for New 

Utility Plants LOW High Medium Low High 
source: Department of Energy computer printout. 
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Table 19--Assumptions underlying all of the scenariosl/ 


1. Crude oil prices are decontrolled by 1981. 


2. The alternative fuel cost for limiting the surcharge added to natural gas 
prices for low priority users is set at the wholesale price of distillate 
fuel oil. 


3. There is a 130 percent lifecycle cost test for new, large oil or gas 
boilers. 


4, Natural gas prices to producers are limited as specified by the Natural 
Gas Policy Act of 1978. 


5. Limited amounts of Canadian and Mexican natural gas and liquefied natural 
gas can be imported. 


6. The price of natural gas imports under new contracts is tied to the world 
oil price. 


7. The Alaskan Natural Gas Pipeline is completed by 1985. 


8. The Trans-Alaskan Pipeline can be expanded to a capacity of 1.6 million 
barrels of oil a day by 1985 and to 2.2 million barrels a day by 1990. 


9. The PACTEX oil pipeline is the only oil pipeline link connecting the West 
Coast to the midwestern and southwestern pipeline networks. 


10. No domestically produced crude oil can be exported. 


11. Petroleum product imports are constrained to encourage the domestic 
refining of oil. 


12. The only new, coal-fueled electric powerplants that will be available in 
1985 are currently in the planning stage. 


13. The rate reform provisions of the National Energy Act do result in the 
use of time-of-day rates that improve electric utilities’ efficiency of 
operation. 


14. Nuclear powerplant construction is constrained by noneconomic factors. 


15. Diesel cars achieve a 9.4 percent share of the market by 1985 and this 
share remains constant thereafter. 


16. Diesel light-duty trucks achieve a 7.8 percent share of the market by 
1985 and this share remains constant thereafter. 


Ly, One further assumption concerning the regional energy price forecasts 
from the DOE should be noted. The DOE forecasts are aggregated from regional 
projections. Conversation with the DOE researchers indicated, however, that a 
national price index was used to maintain constant dollar prices in the fore- 
casts. Thus, the regional price projections do not account for regional 
variation in rates of inflation when adjusting forecasted prices to constant 
dollars. The implicit assumption is that there is no differential in the rate 
of inflation among regions between 1975 and 1985. The alternative would be to 
predict the rates of inflation by region which is beyond the scope of this 


study. 
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Table 20--Projected growth rates of electricity and natural gas price 
for residential use, by DOE region, 1975 to 1985* 


Ratio of 1985 


natural 

Ratio of 1985 electricity gas price to 1975 

DOE region price to 1975 electricity price natural gas price 
New England 0.9195 152797 
New York-New Jersey 0.8935 1.4117 
Mid-Atlantic 1.0394 1.6285 
South Atlantic Wels 1.8346 
Midwest 1.1064 1.7541 
Southwest 1.4329 2.0275 
Central 1.1164 139535 
North Central oll AA 2.3276 
West Manley 22325 
Northwest 1.0891 1.7693 


*Constant dollars 
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pages 9-11 of the DOE report and are summarized in table 18. A more detailed 
description of the DOEs'assumptions are shown in the appendix of the DOE 
report. 


Income Per Family 


Specification of the electricity demand equation required family income 
since that measure is more appropriate as an indicator of ability to pay by 
electricity customers. The consuming unit is more accurately measured as the 
family than as the individual. The Regional Economics Information System of 
the Bureau of Economic Analysis provided computer printouts of historical 
personal income per capita for 1975 and projected personal income per capita 
for 1985. The 1975 income figure was in current or 1975 dollars and the 1985 
income data was in 1967 dollars. The Personal Consumption Deflator is used by 
the Bureau of Economic Analysis to adjust personal income to constant 
dollars. Therefore, adjustment of the 1985 data to 1975 dollars required 
multiplication of all 1985 income data by 1.5295. The calculation of the 
required adjustment to put all income data in 1975 dollars was accomplished by 
using the Personal Consumption Deflator for 1967 and 1975 from the Survey of © 
Current Business. A further adjustment was required because the base year for 
the Personal Consumption Deflator was changed from 1958 to 1972 in January of 
1976. The change of base required an adjustment of 1.382 times the indicated 
ratio of the deflator for 1975 divided by the deflator for 1967. The total 
adjustment was 


1.265 (1,382) = 1.5295. 
1.143 


In addition to the adjustment to constant dollars, the income data were 
weighted to reflect the populations of the service areas of the electric 
utility districts for which the demand forecasts were required. The old BEA 
economic area codes were used in conjunction with computer printouts provided 
by the Bureau of Economic Analysis showing population by economic area. For 
example, the demand Standard Metropolitan Area denoted by our study as 
"Birmingham" consists of BEA Economic Areas 40, 45, and 137 (Montgomery, 
Birmingham, and Mobile, Alabama). The relative share of population was used 
to compute a weighted average income for the utility service area denoted as 
Birmingham. Historical population data from the BEA were used to weight the 
BEA historical per capita income data and BEA 1985 projected populations were 
used to weight the 1985 BEA projected per capita incomes. As noted earlier, 
the correct measure to indicate ability to pay for household electricity is 
assumed to be family income. BEA forecasts of family income were not 
available so the percentage growth of per capita income from 1975 to 1985 was 
calculated from historical and projected per capita incomes provided by the 
Bureau of Economic Analysis. A computer program, designed to predict 
household electricity sales by city, applied the projected rate of change of 
per capita incomes to the family income sample data used in the estimated 
residential electricity demand equation. In this way, BEA data were applied 
to our family income data in conjunction with our estimated demand equation to 
predict future electricity sales. 
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Heating Degree Days. 


The long run average heating degree days was used as the predictor of 
heating degree days in 1985. In the absence of any recognized systematic 
predictor relation for heating degree days, the expected value serves as the 
maximum likelihood estimator. It is, of course, possible that existing or new 
theories of climatic change could allow meteorological experts to devise 
forecasts which could be incorporated into our electricity sales prediction 
model in place of the average heating degree days measure. 


Bias 


As discussed in the previous section, single point estimates are 
substituted for the various independent variables (for a given city) in the 
statistical demand equation. As is well known, the Jensen Inequality implies 
that greater accuracy can be obtained when the prediction equation is 
nonlinear if the distribution of values which will exist over the prediction 
period 1S incorporated into the forecasting process. Our demand equation is a 
double logarithmic transformation and also contains quadratic terms for 
heating degree days. The use of single point projections of the independent 
variables, such as expected price in 1985 or average heating degree days, may 
lead to some bias in the forecast. For example, consider the log-linear 
demand function depicted in Figure 6. Suppose that price in year 1985 takes 
the value P+ A exactly half the time and the value P - a the rest of the 
time. P is the expected price based upon DOE projections. As is shown in 
Figure 6, when the estimated demand equation is convex, the use of the single 
point projection of. price tends to understate the quantity demanded in 
comparison with the prediction which is obtained if the prediction is based 
upon a weighted average of the amounts consumed at the two prices which 
actually exist in 1985. Of course, the latter prediction requires knowledge 
of the probability distribution of price in 1985 which is more difficult to 
estimate than the single expected value. Given the difficulty and complexity 
of arriving at a point price forecast, determination of the probability 
distribution of prices is impractical. Nonetheless, it is important to note 
the direction of the bias inherent in the prediction technique. A downward 
btas of unknown magnitude will exist in the forecast of electricity 
consumption because of the lack of information on the expected variability of 
1985 electricity prices. A similar effect occurs for other variables which 
appear as arguments in the demand equation. It should be noted that it is the 
variability of the independent variables over the course of the year 1985 
within a given city, not across cities, that is relevant. In some instances, 
data to allow an estimate of this variability might be obtained. One example 
of this is in the heating degree days variable. The average variability of 
heating degree days for a given city for a year might be used to project the 
variability for 1985. The heating degree days relationship to quantity 
demanded is also convex as shown in Figure 6 and this again leads to an 
understatement of quantity demanded if the probability distribution of HDD is 
neglected. 


Future Electricity Consumption City by City 


The following tables present predictions of per family and total household 
electricity consumption by city in 1985. In addition to the base consumption 
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Figure 6--The log-linear demand function 
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projections the change in consumption from 1975 to 1985 is partitioned out for 
each of the causal variables. Thus, the separate influences of changes in a 
particular city's energy price structure and income level on electricity 
consumption are indicated. The effects of changes in prices and incomes 
depend both on the elasticity estimated by the regression model and on the 
amount of price and income changes projected by DOE and BEA for 1985. The 
effects of an error in the price or income forecasts are easily measured by 
observing the constant elasticities estimated by the regression model. 
Examination of the estimated coefficients in Model VI reveals that, for each 1 
percent increase in the projected electricity price, the per family 
consumption of electricity will decline by 0.72 percent. Likewise, for each 1 
percent increase in income, per family consumption will rise by 0.35 percent. 
The cross elasticity effect of changes in natural gas price on electricity 
consumption shows that, for each 1 percent rise in gas prices, electricity 
consumption will decline by 0.10 percent. A 1 percent increase in humidity 
will cause electricity consumption to rise by 0.44 percent. The effects of 
changes in the heating degree days projections on the electricity consumption 
forecasts are not easily seen from the regression equations. The quadratic 
equation used to model the effects of HDD on electricity consumption does not 
lend itself to easy interpretation via elasticities. For this reason, the 
computerized prediction program will provide estimates of the effects on per 
family electricity consumption for increments of 10 and 50 percent to Heating 
Degree Days for each city electricity consumption forecast. Although the 
constant elasticity coefficients provide easy interpolation of the effects of 
changes in prices and incomes, a given percentage error in the DOE price 
forecasts results in differing amounts of adjustment depending on the 
magnitude of the original forecasted change for the particular DOE region. It 
is useful, therefore, to present a sensitivity analysis which shows the 
resulting consumption prediction if alzZ DOE price projections are assumed to 
be understated by 10 percent. Likewise, the effects of an understatement of 
the BEA income projections by 10 percent is shown. It should also be noted , 
that total consumption is based upon the BEA projections of population by 
city. The forecast sensitivity to population forecast errors is proportional 
and need not be simulated on the computer, 1.e., a 10 percent understatement 
of projected population simply understates consumption by 10 percent. This 
possible source of error should not be ignored, however. 


Two sample tables, typical of the computer-generated output for 136 
utility service areas in the United States, are shown. The comolete set of 
136 tables is contained in Appendix III. In addition to the salient features 
of the tables discussed above the heating degree days variable requires 
further comment. As shown in Figure 7, the amount of electricity consumed 
declines as cooler weather (higher HDD) is encountered. This is because less 
air conditioning is required. When extremely cold weather occurs, the rela- 
tionship reverses and more electricity is consumed. Electric space heating 
and increased water heating requirements offset the savings in air 
conditioning energy use. Typically, kWh declines from level A to level B in 
Figure 7 when HDD increases by 10 percent. When HDD is increased by 50 
percent a decline in kWh from A to C typically occurs. For a few cities, 
where HDD is above average (cold climate), a 50 percent increase in HDD 
actually results in increased kWh. This is depicted by a movement from A to D 
in Figure 7. The second sample table for observation 53 is an example of the 
latter situation while the sample table for observation 1 is an example of the 
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Figure 7--Heating degree days variable relationship 


former. Cities with cold climates, where a 50 percent increase in HDD results 
in increased KWh, include Duluth, Minneapolis, Fargo, Sioux Falls, Great 
Falls, Manchester, Nashua, Rochester, Milwaukee, Racine, Appleton, and Green 
Bay. 


CHAPTER VI. THE NEED FOR FURTHER RESEARCH 


The primary goal of the research presented here was to project household 
electricity demand in 1985 in specific utility service areas. However, the 
process of model building and empirical testing serves the dual purpose of 
explaining and predicting through the development of cause and effect 
relationships. The broader contribution to theoretical and empirical analysis 
is to increase knowledge concerning the operation and behavior of markets, 
institutions and consumers. A more accurate perception of these elements 
allows decision units to promulgate policies which are appropriate to real 
world problems. The development of more complex models to describe the 
operation of markets for electricity may or may not allow more accurate 
prediction of sales in 1985. It is more likely that the extension to more 
realistic models would enable decisionmakers to better predict the 
ramifications of public policy actions which may impact the functioning of 
markets and consumer decisions. For example, more detail on the reaction time 
to price changes or other exogenous influences can be supplied through a model 
which describes both the time path of change in appliance utilization rates 
and the investment or disinvestment process in electric appliances. 28/ For 
prediction purposes, we have utilized a simplified model which combines 
changes in appliance stocks and changes in utilization rates into the price 
effect. For policy analysis, available data can be utilized to study more 
precisely how consumers react to changing prices and other variables. (1) 


Our research has concentrated on collecting and developing an improved 
Gata set with which to study residential electricity consumption. We have 
replicated previously developed demand models for our service areas and have 
developed improved models. A data set based on observations across Standard 
Metropolitan Areas has been assembled. Part of this data set is derived from 
information, provided by the American Gas Association, measuring prices of 
substitute fuels such as natural gas, propane, heating oil and wood. Most 
other studies have relied upon crude estimates to approximate these sub- 
stitute fuel price effects.29/ The AGA data provide us with superior, 
previously unused observations for residential electricity demand analysis. 


Previous studies have often failed to adjust for differences in the cost- 
of-living between regions. They have also used cross-state observations in 
their statistical analyses. Unfortunately, cost-of-living estimates are 
unavailable on a cross-state basis. Using the SMA observational unit and the 
AGA price data allowed us to develop cross-city cost-of-living data. Each of 
the models we test can be estimated both with nominal price and income 
variables and with monetary variables adjusted for cost-of-living differences 
between SMA's. This is extremely important in the EPA-USDA research because 
of its emphasis on small area interregional competition. The level of detail 
on the demand side must match the level on the supply side. 


28/ According to a study by the Center for Policy Alternatives at the 
Massachusetts Institute of Technology the electrical consumption for 
residences was broken down into the categories shown in Figure 8. 

29/ Previous electricity demand studies suffer from numerous defects, many 
of which are related to the paucity of accurate data. 


a3 


Figure 8--Distribution of residential electricity consumption, 1974 
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Existing studies of residential electricity demand are extremely 
heterogeneous in the structuring of their models and in their research 
Findings. All studies agree that prices are important in projecting the 
future sales of electricity. Due to differing methodologies employed by 
researchers, reflecting differing assumptions concerning the appropriate way 
to describe the behavior of purchasers and sellers of electricity, quite 
different estimates of the effects of price changes on quantity sold have been 
reported. In the process of replicating some of the existing studies using 
Our improved data set new insights suggest certain modifications in the 
structuring of model specifications. Also, combining the ideas of previous 
researchers has resulted in new models to be investigated. Preliminary work 
with some of these models has commenced but much of the investigative work on 
the new models cannot be included in this report. 


We have determined that long run electricity demand projections should 
consider changes in shares of appliances powered by electricity. However, our 
test of this hypothesis30/ depends on a single year (1960) of substitute 
fuel price data (provided by the AGA). A complete and accurate statistical 
analysis of the lagged and cumulative effects of historical fuel prices on 
current and future electric appliance stocks depends on the AGA data. 


Further research support would enable us to (1) continue the development 
of our data set, (2) more fully utilize the superior AGA substitute fuel price 
information, (3) continue the development and testing of several alternative 
models of the market for residential electricity, and (4) estimate the 
coal-fired electricity demand of those utilities not represented by SMA data. 
Continued research would permit the development of a residential electricity 
Gemand model which is more realistic and sophisticated than those now in 
existence. This model promises to provide a more complete understanding of 
the structure of residential electricity demand. It will permit forecasts 
which incorporate far more of the variable factors influencing future 
electricity sales. The understanding of these factors and their relation to 
household electricity consumption should provide useful information for policy 
analysis. 


Our research, to this point, suggests further examination of a two-tier 
model. The first equation relates quantity of electricity demanded per 
customer to the price of electricity, weather conditions, and the level of 
saturation of ownership of residential electrical appliance stocks. A second 
equation relates the saturation levels of ownership for each residential 
electrical appliance to the price of electricity, the prices of competing 
fuels, 31/ income, climate, and other regional characteristics. 


Equations to predict appliance saturation rates for five appliances have 
been estimated. In general, these five equations relate the percentage of 
households with an electrical appliance (automatic washer, air conditioning, 
dishwasher, food freezer, and electric water heater) to the saturation level 
of these appliance stocks ten years earlier, to prices for electricity and to 


30/ See Appendix II. 

31/ Historical fuel prices should and will be incorporated here using the 
AGA price data set. Currently a version of the "Short cut" distributed lag 
model is used to account for the missing historical fuel price data. 
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the prices of competing fuels. (Climatic conditions for the region are 
relevant for some appliances.) Thus far, only simple linear functions have 
been estimated for the appliance saturation prediction equations. More 
sophisticated equations can be developed with the yearly data on fuel prices 
provided by the American Gas Association. 


Exploratory Modeling 


1. Development and testing of a two-tier (recursive) model to forecast 
residential electricity demand in urban areas is described above. This model 
is complete except for historical fuel price data. Experimentation with 1960 
AGA fuel price data has shown that saturation levels of residential electrical 
appliances are partially determined by historic fuel price levels. The AGA 
has supplied yearly observations, by city, as required to complete the 
development of the appliance saturation prediction equations. 

These data have not been previously released for use in statistical modeling 
and no previous electrical demand models have utilized historical price data 
to predict electric appliance saturation levels. 


2. Development of a data set across cities (SMA's) and consisting of 177 
variables with over 300 more variables added by transformations. An 
additional set of variables can be calculated from DOE data to permit the 
testing of a simultaneous equation (supply and demand) model to predict 
residential electricity use. Our data set include price deflators for cities 
which will allow the testing of all models with and without price deflation. 
Previous studies have not usually deflated across Se Sy aa dere they 
have deflated over time. This data set will allow an analysis of the effects 
of interregional price differentials. The data set also includes electric 
bill data provided by the Department of Energy by city. These data will allow 
the testing of hypotheses which argue either for marginal pricing or average 
pricing and/or combinations of fixed charges plus marginal prices. As 
mentioned above, a controversy exists over the proper specification of 
consumer behavior in this regard. 


Summary of Work in Progress and/or Proposed 


1. The two-tier residential demand model has been validated but the 
electrical appliance saturation can be much more accurately specified using 
the historical AGA fuel price data.33 


2. Development of a two equation simultaneous model of supply and demand 
is in progress. Currently, data are being collected which will be used in the 
supply portion of the two equation model. This model will include electric 
appliance saturation rates predicted by the AGA fuel price data. This model 
structure will be a simultaneous supply and demand model with a recursive 
component to show how appliance stocks are affected by historic fuel prices 
and how appliance stocks affect current electricity demand. Several effects 
not previously tested in supply functions will be tested in the simultaneous 
model. Supply variables could include an average weighted fuel input price, 


32/ In fact, only cross-city deflators exist and these may be inappropriate 
for the cross-state data generally used in prior studies. 
33/ See Appendix II. 
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the weighted average age of boilers, a combination of variables which will 
adjust for the plant size mix for any given utility34/ and, of course, the 
output produced. 


3. The development and testing of an alternative demand equation similar 
to one specified by consultants for EPRI (95). This demand equation 
incorporates appliance stocks weighted by their average utilization rates for 
the total U.S. A utilization prediction equation attempts to forecast how 
regional variations in fuel prices, weather, and other regional 
characteristics affect the rate of use of the appliance stocks. Our model 
differs from the EPRI model in that we use city data and the EPRI model 
related utilization to the absolute capacity of appliances, if used 24 hours 
per day, rather than explaining deviations from normal use. 


4. Development and testing of commercial and industrial electricity 
demand prediction equations. Model development by researchers in these 
markets is rudimentary and the disaggregation incorporated by previously 
published models is at the state level. The data set which has been collected 
includes price and quantity data for commercial and industrial demand by SMA. 
We have recently added data to show commercial business activity and 
industrial production by SMA. Thus the data set for this portion of the 
electricity demand analysis is nearly complete. Sorting and merging of data 
sets is all that remains before the statistical estimation analysis can 
proceed. 


5. The two-tier residential electricity demand model can be replicated 
for years other than 1975, the base year for our current model. This will 
allow testing of the hypothesis that consumer behavior is changing over time 
aS awareness of energy shortages and government policies affect consumer 
tastes. 


34/ See Appendix I 
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APPENDIX I 


A TECHNIQUE TO ADJUST FOR PLANT SIZE MIX WHEN ESTIMATING 


LONG RUN COST FUNCTIONS 


LAC 


AC 


ee MEN OER 


OUTPUT 


If we assume that physical generating capacity is roughly equivalent to 
the short run average cost concept, as shown above, then the mix of plant 
sizes owned and operated by an electric utility company affects their overall 


efficiency. Other factors which affect cost and efficiency are ignored in 
this discussion. 


Given that scale affects efficiency, we can write 


SAC = f(c) (1) 
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where AC is the average cost per unit of electricity produced for a given mix 
of plant sizes. The mix of plant sizes can be demarcated as a probability 
distribution Py(c). The expected LAC is the sum of all possible (existing) 
SAC's weighted by the relative frequencies of occurrences of their output 
Capacities, 1.@e., 


E(LAC) = (SAC) (P-(c)) de (2) 
or, by substituting in equation (1), we find 
E(LAC) = /£(c) Pr(c) de. (3) 


We can express the expected value or average of the LAC'’s associated with 
a given mix of plant sizes in terms of the moments of the frequency 
distribution of plant sizes (capacities). To accomplish this, substitute a 
Taylor's expansion for f(c), letting c =C + c' where c' is any deviation from 
the mean plant size c. The Taylor's expansion about the point c is 
f(c) = #(€) + #(c)(c - &) + Af (Ec - &)? + af (E)(c - )° 


bs cus ae (4) 


Expected LAC is found by substituting equations (3) and (4). 


E(LAC) = S£(¢) Pr(c) de + SJ £ (c) (c-€) Py(c) de + 1/2 SE (c) 


(c-€)2 Pe(c) de + 1/6 SE" (&) (c-é)3 Py(c) de +... ad (5) 


The second term of equation (5) is always zero. The first term reduces to 
f(c), the third term becomes 1/2 £” (c) va and so forth. The generalized 
expression for expected LAC can be written: 


E(LAC) = £(c) + 1/2 £" (¢c) va + 1/6 £'" (€) So 


+ 1/24 £"' (c) Ko + -—- (6) 


where c is the average capacity, ve is the variance of plant capacities, 
Sc is the skew and Kc is the kurtosis of the plant capacity distribution 
for a utility. 


A simple way to estimate equation (6) is to regress LAC on the RHS of 
equation (6); this results in E(LAC f(c)). Addition of c, vé, So and 
Ke to a long run cost function should more completely and accurately specify 
the elements affecting costs. If the plant size distribution is normal, then 
only c and Ve are required to completely describe the distribution. 
Also, if the LAC is linear, then only c need be included since f' and higher 
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derivatives in (6) will all be zero. Likewise a quadratic LAC would only 
require that c and v@é be included in the cost curve equation. It seems 
likely that at least one of the above.conditions would hold, reducing the 
necessary scale variables to c and v,”. Whether or not higher order 
moments of the plant size distribution need be included can be determined 
empirically. 
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APPENDIX II 


EXPLORATORY ANALYSIS OF THE TWO-TIERED APPLIANCE STOCK MODEL 
WITHOUT COMPLETE HISTORICAL PRICE DATA 


In order to capture both the effect of electricity price on the saturation 
levels of electric appliances and to incorporate utilization effects of 
Soman pee price changes, we adopted the recursive model discussed 
below.i/. 


Our reasoning leads us to believe that the market is most accurately 
depicted by a two-tiered model. The first equation explains the rate of 
utilization of appliances. It relates the quantity of electricity demanded 
per capita to the price of electricity, weather conditions, and the level of 
ownership saturation. The level of ownership saturation is related to the 
price of electricity, the historical prices of competing fuels, income, 
climate, and other regional characteristics. 


Specifically, the first equation relates per capita electricity consump- 
tion to the percent of houses with electric hot water, electric clothes dryer, 
electric space heat, air conditioning, electric range, automatic washing 
machine or electric food freezer. In addition to the saturation levels of 
these appliances, the model includes current prices and income, several 
measures of weather and state dummy variables which allow for shifts in demand 
or for variation in demand elasticity across regions. Several different 
models have been developed, and these are presented in Tables E-] through E-8 
near the end of this appendix. 


The second equation set estimates the saturation levels of these electric- 
al appliances. This must precede the projections of the long run demand for 
electricity. Equations to predict appliance saturation rates for eight appli- 
ances have been estimated. In general, these eight equations relate the per- 
centage of households with an electrical appliance (automatic washer, air con- 
ditioning, food freezer, electric hot water, etc.) to the saturation level of 
these appliance stocks ten years earlier,2/ to prices for electricity and to 
prices of competing fuels. (Climatic conditions for the region are relevant 
for some appliances.) Thus far, only simple linear functions have been 
estimated for the appliance saturation prediction equations. (More so- 
phisticated equations can be developed with annual data on fuel prices 
supplied by the American Gas Association.) These equations are presented in 
Tables S-1 to S-8 at the end of this appendix. 


Development and initial testing of a two-tiered (recursive) model to fore- 
cast residential electricity demand in urban areas is complete except for 


i1/ Cross section regression is often suggested to measure long run 
equilibrium relationships. If appliance stocks are excluded from the demand 
specification, and if current electricity prices act as a proxy for historical 
Or expected prices, then the estimated price elasticity is an approximate 
measure of the long run sensitivity of sales to price. This kind of estimate 
provides little understanding about the mechanisms by which consumer behavior 
is manifest in the market place. 

2/ A version of the "short cut" distributed lag to substitute for missing 
historical fuel price data. 
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missing historical fuel price data. Experimentation with 1960 AGA fuel price 
data nas shown that saturation levels of residential electrical appliances are 
partially determined by prior fuel prices. The AGA has supplied yearly 
observations, by city, as required to complete the development of the 
appliance saturation prediction equations. 


The following two sets of tables show preliminary equation estimates which 
can be revised to incorporate the AGA fuel price data for 1960-1975. 
The first set of tables, Tables E-1 through E-8, show some alternative ap- 
proaches to the formulation of the short run electricity demand, the demand 
for electricity when stocks of electrical appliances are fixed. These short 
run models all predict the kWh/household for the year 1975. Most of the 
models contain dummy shift variables for each state. Previous research has 
noted a wide variation in price elasticity by region, in some cases our state 
dummy variables adjust the elasticity while in others they shift the entire 
demand function. The final model can be estimated both ways to obtain the 
best fit. 


Model E-1l 


Model E-1 is a common log linear model in which all variables, except the 
dummy variables, are transformed into logarithms. The estimated B coeffi- 
cients can be interpreted directly as elasticities. Price elasticity is -.747 
and income elasticity is +.835; the cross elasticity with fuel oil is -.168. 
The fuel oil price appears to have the wrong sign Since it is not normally 
considered to be a complement to electricity. This equation is one of the few 
regression models to contain income as a Significant determinant of quantity 
demanded. 


Model E-2 


Model E-2 is linear which may reduce the fit and the number of variables 
which entered. A new variable is introduced, LWA, which is a weighted average 
of the percentage of houses having each type of electrical appliance. The 
weights used to aggregate the appliance saturation levels are the national 
average rates of use for the respective appliances. 


Model E-3 


Note that both Model E-1 and E-2 use deflated monetary variables. Model 
E-3 uses nominal or undeflated money variables and this reduces the accuracy 
of the model. Model E-3 is expressed in logarithms so that the B coefficients 
can be interpreted directly as elasticities. The state dummy variables are 
multiplied by electricity price so that the estimated coefficients of the 
dummies are adjustments to the price elasticity for particular states. Note 
that, in place of the weighted average appliance saturation level, the 
saturation percentages for automatic washers and for electric heat are entered 
individually. The saturation rate for electric heat also enters as an 
interaction term with electricity price. It indicates that price elasticity 
is lower in regions where electric heat is more prevalent. Price elasticity 
in this short run model is -.597. 


iu 
Model E-4 


Table E-4 shows the correlations among the various measures of electric 
appliance saturation rates. As a preliminary trial, it was decided 
toaggregate all appliances except air conditioning. Air conditioning has a 
low correlation with most other appliance saturation rates. It is intuitively 
clear that the climatic variables in the models are related to air 
conditioning in different ways than to other appliances. Table E-4 shows very 
high correlations among some saturation variables. High correlations suggest 
aggregation since the separate influences of these variables cannot be 
statistically identified. Future research can investigate alternative schemes 
to aggregate the saturation measures. 


Models E-5 and E-6 


Model E-5 is an appliance utilization rate model. The variables in the 
model regulate the rate of use of WA, the weighted average of all appliance 
saturation rates. The variable LP, for example, is the coefficient of 
electricity price on the rate of use of all appliances. This model is 
disaggregated later into an air conditioning saturation rate versus all other 
appliances. This will be shown in Models E-7 and E-8. 


Model E-6 is similar to E-5 except the state dummy variables do not 
indicate differences in the coefficients of the price variable. Instead, 
these dummies shift the entire equation to account for regional variations in 
electricity demand. Note that, in this case, the demand shift regional 
adjustment appears to fit slightly better. The R¢ rises slightly from .93 
to .94 and the F values for the dummy variables are higher. 


Models E-7 and E-8 


Models E-7 and E-8 are similar to models E-5 and E-6 except that the rate 
of use of air conditioning is assumed to be affected differently by the price 
_ and weather variables than is the rate of use of other electric appliances. 
The variable (WB)(P), for example, shows the effect of changing electricity 
price on non-air conditioning appliance use while the coefficient of (WAC) (P) 
shows the effect of changes in electric prices on air conditioning use. Note 
that the electricity price effect on air conditioning appears to be much 
smaller. Model E-7 is for nominal money values while model E-8 is for 
deflated money values. The R* increases from .679 to .953 when the 
adjustment for cost-of-living is made. The effect of price on the rate of 
utilization of non-air conditioning appliances is lower with deflated data. 


Models S-1 through S-9 


Tables S-l1 through S-9 show preliminary estimates of the long run de- 
termination of electric appliance saturation. In each case, the dependent 
variable to be predicted is the percentage of houses having a particular elec- 
tric appliance.3/ Note that the prices of electricity and of substitute 
fuels enter in many of the equations. The fuel prices in 1960 are used as a 
Surrogate for expected fuel prices; in some cases, the Saturation rate of the 


3/ Weighted aggregates of these electric appliance saturation rates appear 
as independent variables in some of the short run models. In combination, 
they form a recursive system. 
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appliance ten years previously is a significant determinant also. This is a 
form of distributed lag effect which indicates that the fuel price data for 
1960 is not an adequate indicator of fuel price expectations. The detailed 
AGA fuel price data will allow us to construct Almon lag models to 
moreprecisely estimate expected fuel prices. These exploratory models simply 
show the appropriateness of the further development of these ideas. It is 
interesting to note that alternative fuel prices do not enter the short run 
models but, as would be expected, they do enter into the determination of the 
type of fuel to use in many of the appliances. Climatic variables, as opposed 
to the weather variables which entered the short run demand functions, are 
significant in several of the electric appliance saturation models. As 
previously discussed, these latter equations can be combined with the short 
run models to predict the total. long run sensitivity of electricity sales to 
changes in price and income. 
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TABLE E-1l 
DEPENDENT VARIABLE: Log transform of SAMPLE SIZE: 40 SMSA's 
~~ kWh/Customer in 1975 
INDEPENDENT VARIABLES; FUNCTIONAL FORM: 
S's are state shift dummy variables | 
P is deflated electricity price in 1975 LOG LINEAR MODEL - (Coefficients 
OP is deflated fuel oil price in 1975 can be interpreted directly as 
CDD is cooling degree days in 1975, a elasticities except for the state 
weather variable dummies) 


INC is deflated median income per capita 
HDD is heating degree days, a weather 


variable 

Multiple R 992 

R Square -985 

Adjusted R Square 966 
Sa a VARIABLES IN THE EQUATION-—-————-————————————_—————= 
VARIABLE B t-VALUE 

SIGNIFICANCE 

S38 -.96801144-01 - 1.56 
S37 67929893 8.52 
S56 - 26769338 2.11 
S42 —. 72449080 - 2,43 
S58 ~ 37544829 3.67 
S30 -.21719585 - 3.78 
S45 . - 19340005 2.31 
S27 | ~.37033114 = A.10 
S22 ~. 21710953 - 2.60 
S31 ~.21027925 - 2.64 
S23 —. 21468990 - 2.41 
S47 - 36847097 4,35 
S60 —. 41251729 - 6.70 
S44 -91365971 Soule) 
S28 -.20590192 - 2.95 


S34 —.13691083 - 2.77 
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TABLE E-1 (cont'd) 


~---~------------------------ VEMAUSEVSS) Th Wala) POU/UNT Oe 


VARIABLE B t-VALUE 
SIGNIFICANCE 
S54 - 38391326 3.86 
Log P -.74720132 =U2 0) 
Log OP =. 16793127 =) 3.60 
CDD - 23025337 6.24 
Log INC -83480305 2.81 
HDD - 25209084 6.43 


(Constant) -6.2864911 Si Acd sl 
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TABLE E-2 
DEPENDENT VARIABLE: kWh/Customer in 1975 SAMPLE SIZE: 40 SMSA's 
INDEPENDENT VARIABLES: FUNCTIONAL FORM: 
WA is the weighted average electric 
appliance saturation rate LINEAR MODEL 


S's are state shift dummy variables 

P is real (deflated) electricity price 
in 1975 

CDD is cooling degree days in 1975, a 
weather variable 


Multiple R 0.972 
R Square 0.941 
Adjusted R Square 0.926 


~------------------~--------- VARIABLES IN THE EQUATION----------------------—~- 


VARIABLE B t-VALUE 
SIGNIFICANCE 
WA 12.997 4.60 
$37 11.462 7 
S58 5.587 5.96 
P -0.976-01 M0) + S297 
$33 2.134 3.37 
CDD 0.153-02 Di 2ece 
S54 aL? ) 2.87 
S66 -1.087 - 1.92 


(Constant) 8.462 9.20 
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TABLE E-3 
DEPENDENT VARIABLE: Log transform of SAMPLE SIZE: 95 SMSA's 
kWh/Customer in 1975 
INDEPENDENT VARIABLES : FUNCTIONAL FORM: 
% AUT, % EHF are actual appliance saturation 
rates for individual appliances (automatic LINEAR MODEL 


washer and electric heating fuel) 

P is nominal electricity price in 1975 

EPEHF is an interaction between nominal electricity 
price in 1975 and percent houses with electric 
heating fuel 

S's are state shift dummies times the log of 1975 
electricity price. They can be interpreted 
directly as adjustments to price elasticity by 


state 
Multiple R 0.904 
R Square 0.818 
Adjusted R Square 0.796 
wen enn VARIABLES IN THE EQUATION------------------------- 
VARIABLE B t-VALUE 
SIGNIFICANCE 
Log P I BS) 7) - 8.21 
CDD75 0.453-03 9.48 
(HDD) (EH) 0.118-02 4.86 
HUM 75 0.488-02 4.03 
(S37) (Log P) 0.103 4.15 
SAUT70 0.6714 2.81 
MUNI75 -0.485-01 - 2.16 
(S56) (Log P) -0.630-01 - 2.09 
SEHF70 -5.449 - 2.41 
(Log P) (%EHF) 1355 2.10 


(Constant) 2.988 8.60 
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TABLE E-4 


Correlation Matrix of Electric Appliance Ownership Saturation in 1970 


SEWHF ECF SEHF %ED AUT %AC SDW $FF 
SEWHF il 
%ECF 0.848 1] 
SEHF Oo713 0.570 } 
%ED 0.572 0.765 02322 1 
SAUT 0.232 0.491 0.022 0.548 il 
SAC -0.047 —-0.314 -0.060 -0.352 ey “213 a 
2DW -0.005 0.139 0.245 0.167 0.294 -0.107 iI 
SFF 0.347 0.410 0.176 OL 572 0.382 -0.165 0.014 1 
SEWHF = percent houses with electric water heating fuel 
%ECF = percent houses with electric cooking fuel 
SEHF = percent houses with electric heating fuel 
SED = percent houses with electric dryer 
AUT = percent houses with automatic washer 
%AC = percent houses with one or more air conditioners or central air con- 
ditioning 
DW = percent of houses with dishwasher 
FF = percent houses with food freezer 
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TABLE E-5 
DEPENDENT VARIABLE: kWh/Customer in 1975 SAMPLE SIZE: 40 SMSA's 
INDEPENDENT VARIABLES : FUNCTIONAL FORM: 
(HUM) (WA) is Humidity times WA (defined below) LINEAR MODEL 


(P)(WA) is real (deflated) electricity price 
in 1970 times WA 

(CDD) (WA) is cooling degree days in 1975 times 
WA 

(S)(P) are state shift dummy variables times 
electricity price 

WA is a weighted average of all electric appli- 
ance saturation levels for each SMSA...the weights 
used are the national utilization rates of each 


appliance. 

Multiple R O27 

R Square 0.944 

Adjusted R Square 0.929 
een nee eee VARIABLES IN THE EQUATION-----------------------—- 
VARIABLE B t-VALUE 

SIGNIFICANCE 

(HUM) (WA) 0.312 rr Ee Bi 
(S) (537) (P) 0.294 11.65 
(S) (538) (P) 0.157 5.43 
(P) (WA) -0.577 - 6.68 
(S)(533)(P) 0.495-01 3.04 
(CDD) (WA) 0.123-01 3.84 
(S) (544) (P) -0.102 - 2.67 
(S) (531) (P) -0.671-01 ) - 1.89 


CONSTANT 6.143 11.02 


119 


TABLE E-6 

DEPENDENT VARIABLE: kWh/Customer in 1975 SAMPLE SIZE: 40 SMSA's 
INDEPENDENT VARIABLES: FUNCTIONAL FORM: 
Same as previous table except S is a state 
shift dummy, it is not multiplied times LINEAR MODEL 
the electricity price 

Multiple R 0.976 

R Square 02953 


Adjusted R Square 0.940 


----------------------------- VARIABLES IN THE EQUATION--------~-------------—-- 


VARIABLE B t-VALUE 
. SIGNIFICANCE 
(HUM) (WA) | 0.309 8.49 
$37 ti TAS i262 
S58 5.046 , 6.00 
(P) (WA) - 0.544 eon 
$33 2.364 4.29 
(CDD) (WA) 0).122-01 4.18 
S44 aes = a506 


$31 ape Ret oT = 2004 


(CONSTANT ) 3-905 1132 
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TABLE E-7 


DEPENDENT VARIABLE: kWh/Customer in 1975 SAMPLE SIZE: 102 SMSA's 


INDEPENDENT VARIABLES : _ FUNCTIONAL FORM: 
(HUM) (WB), Humidity in 1975 times weighted 
average of electric appliance saturation LINEAR MODEL 


rates excluding air conditioning. 

P nominal electricity price in 1975 

S's are state shift dummy variables 

(WAC) (CDD) Cooling degree days in 1975 times 
the percentage of houses with air condi- 
tioning, a weather variable. 

(WAC)(P) nominal electricity price in 1975 
times percentage of houses with air 


conditioning 

Multiple R 0.839 

R Square 0.704 

Adjusted R Square 0.679 
SS SS SE NSO sUdaS) ION) Dee) MNO] OP UO) 
VARIABLE B t-VALUE 

SIGNIFICANCE 

(HUM) (WB) 0.393 8.56 
P -0.652 = 85h) 
(WAC) (CDD) 0.556-02 7.83 
S37 3.482 8.52 
S44 -6.125 - 3.45 
(WAC) (P) -0.789-01 = ANS 
S61 4.053 2.41 
S38 1S 97 1.86 
(CONSTANT) 6.597 13.43 
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TABLE E-8 
DEPENDENT VARIABLE: kWh/Customer in 1975 SAMPLE SIZE: 40 SMSA‘s 


INDEPENDENT VARIABLES; 
Same aS previous table except all prices are 
in real terms (deflated) LINEAR MODEL 


FUNCTIONAL FORM: 


Multiple R 0.981 
R Square 0.962 
Adjusted R Square 0.952 


-—--------~~----------------- VARIABLES IN THE EQUATION------------------------- 


VARIABLE B t-VALUE 
“SIGNIFICANCE 
(HUM) (WB) 0.384 11.53 
S37 11.964 14.66 
S58 4.502 6.32 
S33 2.663 5.44 
(WAC) (CDD) 0.726-02 8.38 
P - 0.390 =) 621 
(WAC) (P) - 0.821-01 = 4.71 
S44 - 5.148 ~ 3.82 
(CONSTANT) 5.491 13.06 
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TABLE S-1l 
DEPENDENT VARIABLE: % of houses with SAMPLE SIZE: 102 SMSA's 
electric water heating fuel in 1970 
INDEPENDENT VARIABLES : FUNCTIONAL FORM: 
S's are state shift dummy variables 
SEWHF60 is the percentage of houses with LINEAR MODEL 


electric heating fuel 10 years earlier 
(1960, a distributed lag effect) 
P60 is the nominal price of electricity 


in 1960 
GP60 is the nominal price of natural gas 
in 1960 

Multiple R 0.944 

R Square 0.801 

Adjusted R Square 0.873 
a VARTABLES@ IN THE EQUATION-————————— ~~ 
VARIABLE B t-VALUE 

SIGNIFICANCE 

GP60 0.789-02 2.89 
S61 0.696 a7 
SEWHF60 0.204-01 4.14 
S40 -1.161 =3.07 
S54 0.446 8.31 
S67 0.568 6.93 
S47 0.456 8.99 
S31 0.164 4,46 
S48 0.424 | Ge2i. 
S68 0.487 5.44 
S44 0.406 S55 
S49 O223 4.68 
S34 0.113 2016 
P60 -0.379-07 =o, 


(CONSTANT) 0.658-01 2.64 
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TABLE S-2 
DEPENDENT VARIABLE: % of houses with an SAMPLE SIZE: 102 SMSA's 
electric dryer in 1970 
INDEPENDENT VARIABLES : FUNCTIONAL FORM: 
S's are state shift dummy variables 
AVHDD the average heating degree days, LINEAR MODEL 


P60 is nominal electricity price in 1960 

BEARLY is the percentage of houses built 
before 1939 

AVDB iS average temperature uSing a dry 
bulb thermometer 


Multiple R 0.878 
R Square 0.771 
Adjusted R Square 0.731 


----------------------------- VARIABLES IN THE EQUATION-----------~----------~--- 


VARIABLE B t-VALUE 
SIGNIFICANCE 

AVHDD 0.478-04 6.05 
P60 -0.612-07 5.34 
$25 0.357 6.47 
S68 i Qi255 4.09 
S67 0.275 4.97 
BEARLY -0.341 =a 
S61 0.212 3.87 
S50 O95 3.56 
AVDB -0.5608-02 —3.70 
S33 0.101 3.66 
S54 0.104 Sizo 
S62 0.123 S207 
S39 0.114 2290 
S31 0.702-01 2.80 
S29 0.614-01 Zoer 


(CONSTANT) 0.804 5.54 
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TABLE S-3 
DEPENDENT VARIABLE: % of houses with SAMPLE SIZE: 102 SMSA's 
electric heating fuel in 1970 
INDEPENDENT VARIABLES : FUNCTIONAL FORM: 
S's are state shift dummy variables 
B6568 is the percentage of houses built LINEAR MODEL 


between 1965 and 1968 
AVHDD is the average heating degree days 
GP60 is the nominal price of utility gas in 1960 


Multiple R 0.865 
R Square 0.749 
Adjusted R Square 0.731 


—--------~---—--—--—-------—- VARIABLES IN THE EQUATION-~---------------------—- 


VARIABLE B t—VALUE 
SIGNIFICANCE 
S54 0.311 9.88 
S44 0.242 5.02 
B6568 0.578 3.48 
S67 0.199 3.73 
S68 0.188 3.47 
AVHDD -0.180-04 -3.24 
GP60 0.423-02 725839) 


(CONSTANT ) 0.655-02 0.21 
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TABLE S-4 


DEPENDENT VARIABLE: % of houses with SAMPLE SIZE: 102 SMSA‘'s 
air conditioning in 1970 

INDEPENDENT VARIABLES: 

AVDB is average dry bulb temperature 

AVWB iS average wet bulb temperature 

P60 is the nominal electricity price in 1960 

S's are the state shift dummy variables 

B6568 is the percentage of houses built 
between 1965 and 1968 


FUNCTIONAL FORM: 


LINEAR MODEL 


Multiple R 0.774 
R Square 0.594 
Adjusted R Square 0.573 


~-------~-----—-------------- VARIABLES IN THE EQUATION-----------------------~- 


VARIABLE B t-VALUE 
SIGNIFICANCE 
AVDB 0225-01 6.35 
AVWB 0.164-01 4,42 
P60 0.115-06 4.03 
S54 -0.230 2004 
B6568 0.828 “rel! 
S31 0.136 2.08 


( CONSTANT ) 


Deg a0) 
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TABLE S-5 
‘DEPENDENT VARIABLE: % of houses with electric SAMPLE SIZE: 102 SMSA's 
cooking fuel in 1970 
INDEPENDENT VARIABLES : FUNCTIONAL FORM: 
P60 is nominal electricity price in 1960 
B6064 is the percentage of houses built LINEAR MODEL 


between 1960 and 1964 
%ECF60 is the percentage of houses with 
electric cooking 10 years earlier (1960, 
a distributed lag effect) 
GP is the nominal price of utility gas in 1960 
S's are state shift dummy variables 


Multiple R 0.891 
R Square 0.794 
Adjusted R Square 0.749 


-------------------——------- VARIABLES IN THE EQUATION------------------—----——- 


VARIABLE B t-VALUE 
SIGNIFICANCE 
P60 -0.127-06 -6.38 
B6064 1.067 4.62 
SECF60 0.122-01 Sia iky 
S40 | -0.794 -2.28 
S61 0.414 4 
S58 -0.129 -3.23 
S56 -0.202 -3.44 
S25 -0.371 3.77 
S67 0.378 B83 
S54 O27 3.69 
S68 0.333 3.25 
S27 -0.213 ~2.99 
GP60 0.145-01 fo) 
S42 -0.235 | -2.35 


S62 0.201 2.82 
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TABLE S-5 (cont'd) 


VARIABLE B t-VALUE 
SIGNIFICANCE 
S47 0.171 2.80 
S48 0.199 2.62 
S64 0.241 2.43 


(CONSTANT ) 0.184 3.82 
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TABLE S-6 
DEPENDENT VARIABLE: % of households with SAMPLE SIZE: 102 SMSA's 
automatic washing machines in 1970 

INDEPENDENT VARIABLES: FUNCTIONAL FORM: 
S's are state shift dummy variables 
P60 is the nominal electricity price in 1960 LINEAR MODEL 

Multiple R 0.793 

R Square 0.568 


Adjusted R Square 0.516 


~-----~-----~~---~---—------- VARIABLES IN THE EQUATION-———————_-- === 


VARIABLE B t-VALUE 
SIGNIFICANCE 

S25 lls SAIL —6.67 
P60 -0.398-07 -4.13 
S21 0.732-01 Zoo 
S56 0.703-01 2.49 
S40 -0.143 21196 
$3] =O 16-0 8) 5 fas 
S30 -0.791-01 a 713) 
S58 -0.489-01 -2.58 
538 -0.687-01 -2.44 
S41 -0.108 21126 
$35 -0.675-01 TOO 


(CONSTANT ) 0.658 | 105,92 
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TABLE S-7 
DEPENDENT VARIABLE: % of houses with food SAMPLE SIZE: 102 SMSA's 
reezer in Go ae ee 
INDEPENDENT VARIABLES : FUNCTIONAL FORM: 
S's are state shift-dummy variables 
P60 is the nominal price of electricity in 1960 LINEAR MODEL 


%FF60 is the percentage of houses with a 
food freezer 10 years earlier (1960, a 
distributed lag effect) 

BEARLY is the percentage of houses built 


before 1939 

Multiple R 0.942 

R Square 0.888 

Adjusted R Square 0.855 
oe VARIABLES IN THE EQUATION------------------------- 
VARIABLE B t-VALUE 

SIGNIFICANCE 

S25 0.654 15369 
P60 -0.631-07 ="7 60 
$22 -0.104 - 4,58 
S24 2535 - 4,54 
S62 0.158 5.47 
S61 0.190 4.69 
S56 0.102 4,24 
SFF60 0.356-02 4,37 
S44 -0.146 - 4.86 
S27 eee cone 
S21 -0.768-01 lps Ful, 
S34 0.883-01 4.16 
53) 05992-01 4.11 
S67 0.132 Suet 


S23 -0.876-01 =e 


———— VARIABLES IN THE EQUATION 


VARIABLE 


S53 
S59 
BEARLY 
$52 
S49 
S41 
S30 


S68 


(CONSTANT ) 


TABLE S-7 (cont'd) 


B 


0.877-01 
-0.955-01 
-0.111 
=) Syl eh 
asi / 0h 

0.917-01 

0.515-01 

0.837-01 


0.325 
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t-VALUE 


SIGNIFICANCE 


Za 
= 2220 
= eideis: 


- 2.48 
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TABLE S-8 
DEPENDENT VARIABLE: % of houses with a SAMPLE SIZE: 102 SMSA's 
Gishwasher in 1970 
INDEPENDENT VARIABLES: FUNCTIONAL FORM: 
S's are state shift dummy variables 
B6064 is the percentage of houses built between LINEAR MODEL 


1960 and 1964 
BEARLY is the percentage of houses built 


before 1939 

Multiple R 0.916 

R Square 0.840 

Adjusted R Square 0.819 
= SSS 55 5 55355 SS Se yea StS) TES) AREY 0 
VARIABLE B t-VALUE 

SIGNIFICANCE 

S25 0.688 17.06 
B6064 | 1.055 5.18 
S67 0.196 4.90 
S68 0.192 4.79 
S60 0.118 4.08 
S66 0.797-01 $:95 
S42 0.146 3262 
S61 0.130 3.23 
S58 0.617-01 3.40 
S45 0.113 2279 
S21 0.661-01 2579 
BEARLY 0.149 2.48 


(CONSTANT) eo -5 tla 2 


24. 


2B. 
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APPENDIX III 
DEFINITION OF VARIABLES IN EXPLORATORY DATA SET 

1975 Gas Price (AGA Confidential) 
1975 Propane Price (AGA Confidential) 
1975 Electricity Price (AGA Confidential) 
1975 No. 2 Oil Price (AGA Confidential) 
1975 No. 1 Oil Price (AGA Confidential) 
1975 Coal Price (AGA Confidential) 
Population per Square Mile 
Median Age 
Median Income 
1975 Average Humidity 
1975 Heating Degree Days 
1975 Cooling Degree Days 
Average Heating Degree Days 
Percent Sunshine 
Average Dry Bulb Humidity 
Average Wet Bulb Humidity 


Elevation 


1960 Gas Price (AGA Confidential) 


1960 No. 2 Oil Price (AGA Confidential) 
1960 No. 1 Oil Price (AGA Confidential) 
1960 Coal Price (AGA Confidential) 
1960 Coal Price (AGA Confidential) 
1960 Coal Price (AGA Confidential) 
1960 Coal Price (AGA Confidential) 


1960 Propane Price (AGA Confidential) 


26. 
Zi 
20. 
Zo. 
30. 
aie 
32. 
33. 
34. 
35. 
36. 
Sie 
38). 
398 
40. 
4l. 
42. 
43. 
44, 
A5. 
46. 
47. 
48. 
49. 
50. 
Sn Oe 


O2s 
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1960 Coal Price (AGA Confidential) 


1960 Electricity Price (AGA Confidential) 


Total Houses in 1960 


Number of houses with 


Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 


Number 


of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 


of 


houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 


houses 


with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 


with 


gas heat in 1960 

oil heat in 1960 

coal heat in 1960 

electric heat in 1960 
propane heat in 1960 

other heat fuel in 1960 

no heat source in 1960 

gas water heat in 1960 
electric water heat in 1960 
coal water heat in 1960 
propane water heat in 1960 
oil water heat in 1960 

other water heat in 1960 

no water heat in 1960 

gas cooking fuel in 1960 
electric cooking fuel in 1960 
propane cooking fuel in 1960 
oil cooking fuel in 1960 
coal cooking fuel in 1960 
other cooking fuel in 1960 
no cooking fuel in 1960 
wringer type washer in 1960 
automatic type washer in 1960 


washer-dryer in 1960 


Dee 


54. 


Dey 


56. 
Bho 
8p 
aos 
60. 
6l. 
62. 
6S. 
64. 
Da 
66. 
67. 
68. 
69. 
70. 
Wile 


W2o 


US. 


74. 


Ue 
76. 
vue 
18. 


We 


Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 


Number 


1975 Total Revenue from Electricity Sales to the Industrial Sector 


of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 


of 


houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 


houses 
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with no washer in 1960 
with gas dryer in 1960 
with electric dryer in 1960 
without a dryer in 1960 
with a food freezer in 1960 
without a food freezer in 1960 
with a single air conditioner in 1960 
with 2 air conditioners in 1960 
with central air conditioning in 1960 
without air conditioning in 1960 
with one television set in 1960 
with 2 television sets in 1960 


without a television set in 1960 


Dummy Variable indicating Municipal Ownership 


Total Electric Bills 500-750 block for 1975 


Total Electric Bills 750-1000 block for 1975 


Total Electric Bills 1000 + block for 1975 


Total Houses in 1970 


Number of houses with one air conditioner in 1970 


Number 
Number 
Median 
Number 
Number 


Number 


of houses with 2 air conditioners in 1970 


of houses with central air conditioning in 1970 


rent in 1970 


of residential electric customers in 1975 


of commercial electric customers in 1975 


of industrial electric customers in 1975 


kWh consumed by residential customers in 1975 
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80. kWh consumed by commercial customers in 1975 

81. kWh consumed by industrial customers in 1975 

82. 1975 Total Revenue from Electricity Sales to Residential Sector 
83. 1975 Total Revenue from Electricity Sales to Commercial Sector 
84. Total number of dwellings in 1970 

85. Dwellings built between 69-70 

86. Dwellings built between 65 and 68 

87. Dwellings built between 60 and 64 

88. Dwellings built between 50 an 59 

89. Dwellings built between 40 and 49 

90. Dwellings built prior to 40 

91. Total houses in 1970 

92 Number of houses with gas heating fuel in 1970 

93. Number of houses with oil heating fuel in 1970 

4. Number of houses with coal heating fuel in 1970 

95. Number of houses with wood heating fuel in 1970 

96. Number of houses with electric heating fuel in 1970 

97. Number of houses with propane heating fuel in 1970 

98. Number of houses with other heating fuel in 1970 

99. Number of houses with no heating fuel in 1970 

100. Number of houses with gas water heating fuel in 1970 
101. Number of houses with oil water heating fuel in 1970 
102. Number of houses with coal water heating fuel in 1970 
103. Number of houses with wood water heating fuel in 1970 
104. Number of houses with electric water heating fuel in 1970 
105. Number of houses with propane water heating fuel in 1970 


106. Number of houses with other water heating fuel in 1970 


107. 
108. 
109. 
110. 
IILke 
A Ze 
ILS) 
114, 
1Ol)e 
116. 
JULI 
118. 
AUG 
120. 
21. 
WAZ. 
23: 
124, 
25%. 
126. 
U27 2 


2G; 


Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 
Number 


Number 


of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 


of 


houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 
houses 


houses 


with 
with 


with 


136 
no water heating fuel in 1970 
gas cooking fuel in 1970 


electric cooking fuel in 1970 


propane cooking fuel in 1970 


with 
with 
with 
with 
with 
with 
with 
with 
with 
with 
with 


with 


oil cooking fuel in 1970 
coal cooking fuel in 1970 
wood cooking fuel in 1970 
other cooking fuel in 1970 
no cooking fuel in 1970 
wringer type washer in 1970 
automatic type washer in 1970 
no washer in 1970 

gas dryer in 1970 

electric dryer in 1970 

no dryer in 1970 


a dishwasher in 1970 


without a dishwasher in 1970 


of houses with a food freezer in 1970 


ofhouses without a food freezer in 1970 


of houses with one television set in 1970 


of houses with 2 television sets in 1970 


of houses without a television set in 1970 


137 


APPENDIX IV 


ELECTRICITY CONSUMPTION FORECASTS, INDEPENDENT VARIABLE PROJECTIONS, 
AND SENSITIVITY ANALYSIS FOR 136 UTILITY SERVICE AREAS 


NAMES ALABAMA POWER CO. 
CITY NAME’ BIRMINGHAM 


UTILITY 


1 


UPSceRVATION NUMBeR IS! 


4) 


(Fed REGION 


TABLE l=-INDEPENDENT VARIABLE PROVECTIONS AND THEik IMPACT 


Sem owre 2 2Oe e228 OVO S2S SST 2 MOB BTF@S22oo 22S eeewoeos Boece moose owes ee aooas 


N 
N 
1580. 


STATE DUMMY VARIABLE IS NOT RELEVANT 


MONETARY VARIABLES ARE DEFLATED BY A SURROGATE PRICE INDEX 


DEMAND COMPONENT ANALYSIS 


RICE 
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BS Chance IN AAbORsL CaS 


521905 CHANGE IN HEATING DEGRE 


Cele CHANGE IN ELECTRICITY PRICE 
a 

43218 

7 


POrnm 
mt ty oo 


Sees es Sooo Ss Ses Soe Seer eeeoeoes a= e@eeoooowo ene ses Soa =a 
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TABLE 2@=eLECTRICITY CONSUMPTION FORECASTS 


SSS SSS OS SS SSS 2 SSS 28S SCT S 2S V@SVSVC SES STO O22 eB TSS eoroosasoersoewea mesos eao one oem omeoees ee oe aoe eo oan 


TASLE 3—-=SENSI TIVITY OF DEMAND TO VARIATION IN THE INDEPENDENT VARIABLE PROJECTIONS 
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CALCULATED ON THE 1975 BASE. 
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ARIZONA PUBLIC SERVICE CO. 


PHOENIX 
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2 
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TABLE L==INDEPENDENT VARIABLE PROJECTIONS AND THEIR IMPACT 
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TABLE 3==SENSITIVITY OF DEMAND TO VARIATION IN THe INDE PEND 
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TABLE 2==tLECTRICITY CONSUMPTION FORECASTS 


Perse SS lS SSS Sol @ O22 22 SPB V2 SSS SBT SSD OTB OS S28 SSeS @BVeeeroeoeocwro see Genrer esenoesonsmws seo eamoe ee ce 


TABLE 3-—SENSI TIVITY OF DEMAND TO VARIATION IN THE INODEFENDENT VARIABLE PROJECTIONS 
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TABLE 2=-ELECTRICITY CONSUMPTION FORECASTS 
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eS Se Oe OS ae NS Oe ee EP EP Se SS ES SS SE SD GE Ge GE Ee Oe SD SF SS Ow CS GS OE SS SS ES OR ES SS SS SS SSS OSTEO SOS SSS OSES SS OSE SST SOS Swe ee ewe =e 


TABLE 3==SENSITIVITY OF DEMAND TO VARIATION IN THE 
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